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Explicit calculations are made to examine whether the 
tetrahedral model of CH, is really the most stable in the 
Hund-Mulliken and Slater-Pauling theories. The secular 
determinant involved in the H-M procedure is too compli- 
cated (degree 8) to solve generally, but three methods of 
approximation are given applicable in three different 
limiting cases, in ali of which the tetrahedral model proves 
to be that of least ener:xy quite irrespective of repulsions 
between the H atoms. The fact that two of the methods 
do not assume the carbon s-p separation to be small 
shows that s-p hybridization is not a necessary condition 
for tetrahedra! valence bonds. Calculations are also given 
which show that in the Heitler-London-Pauling-Slater 
method the regular tetrahedron is superior to other 
models of somewhat lower symmetry. It is further calcu- 
lated according to both the H-M and H-L-P-S schemes 
that in compounds of the form CH2X2, CHX;, and CH;X, 
the most stable models are tetrahedra of less symmetry 
than the regular tetrahedron; i.e., models with the valence 


angles somewhat distorted from 109.5° unless the C—H 
and C—X bonds are of equal intensity. The predicted 
directions of the deviations from 109.5° are in agreement 
with x-ray diffraction data for CH2Cl. and CHCl;. The 
conclusion is rather prevalent in the literature that with 
s-p hybridization and electron pairing two bond axes tend 
to set themselves at 109.5°. This is shown incorrect; 
instead the angle can be anything between 90° and 180° 
depending on the relative intensities of the s and p bonds. 
If the s bonding power is not negligible, the angle between 
an NH axis and the pyramidal axis in NH; should be 
somewhat greater than the value 54.7° which is obtained 
if the three NH axes are orthogonal and which is char- 
acteristic of pure p-bonding. Actually Dennison and 
Uhlenbeck find 68°, and Lueg and Hedfeld 73°, from 
band spectra. It is shown that CH,* should be a flattened 
rather than regular tetrahedron, conceivably even being 
plane. Also CH; should be a flatter pyramid than NHs3. 





N Part I' a general comparison was given of the 

Hund-Mulliken and Slater-Pauling? theories 
of valence in carbon compounds. It was shown 
that the H-M procedure gives excessive accumu- 
lation of electrons on one atom, while the S-P 
viewpoint loses rather too much of the symmetry 
properties characteristic of an individual electron 
in fields of tetrahedral symmetry. We do not 
intend these remarks as disparaging criticisms, 
for insuperable mathematical difficulties prevent 
any really accurate solution, and both methods 
are to be admired for the headway they have 
made in the understanding of carbon valencies. 


‘J. H. Van Vleck, Part I, J. Chem. Phys. 1, 177 (1933). 
*See Part I for references to the literature. 
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However, since none of the procedures are 
strictly correct, unless modified by exceedingly 
difficult higher approximations, it appears de- 
sirable as one aim of the present paper to 
prove that the tetrahedral model of CH, is 
the most stable by as many methods as 
possible, rather than to accept such a model 
solely because of the tetrahedral geometry of 
the sp* wave functions of Pauling and Slater. 
If, as it turns out, the tetrahedral configuration 
invariably has the least energy in the various 
feasible limiting cases represented in the different 
theories, we can feel reasonably sure that this 
configuration is also the best in the intractable 
intermediate case actually realized. Although we 


220 J. H. VAN VLECK 


shall give most of our space to the H-M scheme 
of attack (Sections 1-4), we do not wish to 
convey the impression that it is superior to the 
H-L-P-S one (Section 5). The H-M procedure 
is adaptable formally to a greater variety of 
models, but the H-L-P-S is probably more 
satisfactory if numerical success is desired. For 
numerical rather than formal results, everything 
in the H-M method would depend on a happy 
choice of a self-consistent field. 

Before proceeding to the details of the calcu- 
lations we shall state certain simplifications 
which we shall make throughout without further 
mention. We shall neglect all types of forces 
between H atoms, notably the Coulomb repul- 
sions of the H nuclei in polar models, and terms 
due to mutual overlapping of the H wave 


of the tetrahedral model since in the tetrahedron 
the H atoms are evidently in the aggregate as 
far separated as possible. It is, however, of 
interest to show that the tetrahedral model 
evolves itself in virtue of the directional valence 
of the central C atom without invoking any 
H-H repulsions. To avoid excessive complexity, 
it will be supposed that all four H atoms are 
equidistant from the C atom. We shall always 
make the usual neglect of excited states of the 
C and H atoms. That is to say, we assume that 
in expanding the complete wave function in 
terms of unperturbed C and H wave functions, 
the coefficients of the 3,4..quantum states of C 
and 2,3..quantum states of H are comparatively 
small. Consideration of these higher states 
would be required if we were to include cor- 











rections for polarization. The very small bonding 
effect of the K shell of carbon will be disregarded. 


functions. These H-H forces are doubtless 
repulsive,* and so only accentuate the stability 


1. First STAGES OF CALCULATION ACCORDING TO GENERALIZED HuND METHOD 


Following Hund, we shall assume that the molecular wave functions of CH, can be expressed to a 
sufficient approximation as linear combinations of the eight wave functions 


¥(1s; Hi), @=a, 6, c, d); ¥(2s; C), ¥(2po,; C), ¥(2po,; C), ¥(2po,,C). (1) 


Here ¥(2p0,; C) means a carbon 2 state which has m,;=0 when viewed along the x axis. The func- 
tions (1) are to be taken as solutions of the equations 


[v?+ (822m /h?)(W.° — Vc) JW (2k; C) =0 (R=s or fp), 


where Vy and V¢ are the potential energies of an electron in the field of an H and C atom, re- 
spectively. These potential energies need not be considered those due to a bare H nucleus and a 
C nucleus surrounded by — 6e, but rather these nuclei surrounded by an appropriately chosen 
centro-symmetric electron density clouds, so that ¥(1s; H) is not necessarily an ordinary hydrogenic 
wave function. The question of the best choice of the cloud need not concern us, as our aim is a 
proof of the tetrahedron’s stability rather than numerical calculation of the binding energy. The 
cloud screening the C nucleus will be greater than — 6e if it attracts much charge from the H atoms, 
making the molecule appreciably polar. The choice of the screening cloud must, for tractability, 
be taken independent of the spatial arrangement of the H atoms. It may be urged that the ‘“‘best 
cloud’”’ does actually depend on the arrangement; but this does not vitiate our proof, since with 
any given cloud we will prove the tetrahedral model that of least energy. The embodiment of electron 
repulsions through a statistical charge density or cloud is essentially the Hartree method of the 
self-consistent field except that, as noted to the writer by Professor Kramers, our procedure does 
not represent the Hartree method in its full generality. Hartree uses different clouds for obtaining 
different solutions corresponding to various classes of electron orbits, whereas we must use the 
same screening cloud for all eight solutions of our forthcoming secular problem. 

We shall use the following notation for integrals which will enter in the construction of the secular 


determinant: 


[v?+ (822m /h?)(W* — Vu,) W(1s; Hi) =0, (2) 


3 Cf. Heitler and Rumer, Zeits. f. Physik 68, 12 (1931) and discussion to be given in Part III. 
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=f f fvass H,)?V cde, C= ~ ff foros Cade 


D=—{ f {¥@siC¥QeuC)Vude, e=—f ff vs: Hobs C)Vnde, 


v 


(3) 


— f { [vas:Hov2ks C)Vods, Tr= f f U(1s; H)y(2ke; Cd, 





with k=7, o, or s. Note that 0,=R,=T,=0. 


The notation 7;, ¢; means that the spatia! quantization is relative to the axis C — Hi, and the indices 
pr, po are here and henceforth abbreviated to z,o. A zm state has m?=1; a o state m,;=0. The 
vanishing of Q, and R, follows from the axial symmetry of Vz,, Vc on integrating over the azimuthal 
coordinate ¢g involved in dv=r’ sin @drdédy. The integrals are independent of 7 because all four 
H atoms are equidistant from C. 

The expressions C, D, Q, T, defined in (3) may all be considered positive. Here and elsewhere 
we omit subscripts when our remarks apply collectively to the 7, o, and s types. It is true that 
some of the integrals are really arbitrary in sign because of arbitrary phase factors in the wave 
functions, but the final formulas which we shall obtain for the energy will naturally involve these 
expressions only through products (e.g., Q,7.) which are invariant of the choice of phases. It is 
convenient, however, to take these expressions individually positive. To see that this is possible, 
one notes that ¥(1s;H;) is nodeless and so may be taken everywhere positive, while ¥(2s;(C), 
¥(20;;C) can be taken positive where they overlap the atom H,, which will be well beyond the 
node of ¥(2s;C). The potential energy Vy, is, of course, everywhere negative, since the charge 
screening the H nucleus is surely less than that of a proton. The signs of B, R are immaterial for 
our purposes. 

The model will not be assumed tetrahedral. Instead it will be supposed that an axis C—H; (t=a, 
b,c, d) has arbitrary direction cosines /;, m;, n;. The secular equation is then 








E st Sm Sn @ @ @ @ 
St Ft fim tin Qi tte 
Sm Tim F,, ’mn Om Dm Cm dm 
Sn Tin mn Pa an b, Cn dy, 
=0, 4 
q ai Am An G 0 0 0 ( ) 
q bh bn bn O GO O 
qd Ci Cm Cn 0 0 G 0 
qd di dm dy, 0 0 0 G 
with the abbreviations 
a =l,[—-Q,+(W,°-—W)T-], g=—-Q0,+(W.°-—W)T,, 
E=W,°—4C,— W, F,= W> P —4C,+(C.—-—C,) > il?-—W, (5) 
G=W*+B-W, rim = —(C,—C,) >-1hmi, sp= —D>I;. 
It is to be understood that bn, --- d, denote expressions similar to that in (5) for a; except that 


mM, +++ Nq replace 1;. Analogous remarks apply to the rest of the notation (5). 

That the direction cosines enter in the way that they do in (4) follows from the fact ¥(2s; C) 
is centro- -symmetric, and that the wave functions ¥(20,; C), ¥(20,; C), ¥(20.;C) are of the form 
*f(r), f(r), 2f(r), and hence transform from one coordinate system to another like the components 
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of a vector.’ For instance 
V(202; C) =ly(20;; C) + (1-12) Y(2 73; C). (6) 


Eqs. (5) show that tems —WT,, — WT, involving W appear off the diagonal in (4). Such terms 
have arisen because the C and H wave functions are not orthogonal. In place of the formula for 
gq in (5) we could equally well have set g= —R,+(Wi”—W)T., as Eqs. (2) and the self-adjoint 
property of the Hamiltonian operator H give us the relation 


weT.—Q= ff [¥(is:HoBYQs: Cdo= ff [¥(2s; OHY(ts; H)do=W"T,—R, (7) 


Analogous remarks apply to the definitions of a, b, c, d. We are using only real wave functions, 
which involves no loss of generality in the absence of magnetic forces. Hence (4) is symmetric as 
well as Hermitean. 

Eq. (4) differs from Eq. (2) of Part I only in the system of representation, as (2) of I was based 
on sf® functions for carbon instead of the usual nonhybridized functions (1). The zero elements 
occurring in (4) are due to our agreement to neglect H-H forces. 

Twice the sum of the four lowest roots of (4) is the energy of the CH, molecule. The sum is to 
be doubled because the Pauli principle allows two electrons to each solution of the orbital wave 
equation, so that only four roots in all are needed to absorb the eight electrons available from the 
H atoms and the L shell of C. 

Eq. (4) is still too complicated to solve except with simplifying assumptions. 


2. THE ‘““NEAR-IONIC” CAsE W¥—W°SDQ 


This case would be realized if CH, were nearly of the ionic structure C*-(H,)**, for it requires 
that the difference between the electron affinities of the C and H attracting centers (appropriately 
screened) be large compared to the various integrals (3), notably the resonance terms Q, R which 
are the main cause of the molecular binding. It further requires W” —W°>(Q so that electrons are 
in the mean attracted from H to C rather than vice versa. Actually this case will be a rather poor 
approximation in CH,, for the ionization potentials of C and H are nearly equal, and their electron 
affinities small. Nevertheless it is instructive to carry through a calculation for this case, because 
it is desirable to solve the secular determinant (4) for as many different limiting cases as possible, 
and especially because it here proves possible to obtain simple and illuminating formulas for the 
energy with arbitrary arrangement of the H atoms. 

- The mathematical simplicity of the present case rests in our ability to use a series development 
in the ratios of the integrals (3) to W7— W¢. It is convenient to introduce a parameter \ into Eq. 
(4) by the following scheme: multiply each Q or T term in (5) by X, and each C or D term by d*. The 
series development may then be regarded as a power series in \. The convergence of the development 
is, of course, due to the smallness of the expressions (3) rather than of X, since finally one must 
set \=1. The réle of \ is thus only a formal one, and is primarily to describe the degree of approxi- 
mation which first yields the tetrahedral model. The calculation will yield the following interesting 
result: the molecule’s energy is a series in even powers of ) and is independent of the way the H 
atoms are arranged if the development is carried only to terms in )? inclusive, but is a minimum 
for the tetrahedral model when terms of order 4 are included. These \* terms are of two types. 
One kind arises from fourth powers of Q, T; the other involves the product of a first power of C, D 


‘In using (6) it must be mentioned that to specify a | elements 7 in (4), use must be made of the fact that if 
real + wave function completely, it is necessary to specify | yx, Wx’ denote « wave functions of C with different 
the phase ¢ in its azimuthal factor cos (g—e). In (6) the | phases, 


value of ¢ is zero if the apse line for measuring ¢ is taken . : 
in the plane x—C-—zi. In calculating the values of the SS [vere Vyz,du= —C, cos (e—e’). 
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and a squared term in Q, T. (Squared terms in C, D prove to be absent.) Both types need not be of 
coordinate importance, but fortunately both yield the same dependence of energy on angular 
arrangement, favoring the tetrahedron. It appears impossible to say which type is the more im- 
portant, as an answer to this question would require knowledge of the magnitude of C, D compared 
to Q. One knows, to be sure, that Q is large compared to C, D, but this is insufficient. 

We now proceed with the details of the calculation. The first step is to apply a transformation 
which will for our purposes reduce the eighth order secular determinant to a fourth order one. 
The transformation in question is 


S(ij) = — S(jt) = — K (ij) /hv(ij), S(ii’) = 6,"4, S(ij’) = 6)" (8) 


where 6 is the usual Kronecker symbol and the K(mn) are defined by writing the matrix elements 
of (4) as K(mn) — W6,”. Note that because of the T terms in (5) the K(mn) invoive W. The indices 
i, i’ refer to any of the C levels and j, 7’ to any of the H ones (i.e., the ground states of any of the 
four H atoms). One can verify that in the transformed matrix S-'KS the elements 7j (which replace 
a, b,c, d, q in (4)) are of the order )* and so affect the four lowest, i.e., perturbed carbon, roots of 
(4) only in the approximation \°. As we shall stop with \‘, we therefore need retain in this trans- 
formed matrix only the fourth order piece (ii’) which embraces the four C roots. On calculating® 
this piece by means of (8) one finds that it gives the secular equation 











st | L(it’) —6*W| =0 (9) 
wit 
Kai)? KGnK Gi) KG KiKi) Kaj K Gb) 
ie — type 
hv(2j) hv(tj)hv(2’7) hv(1j)hv(0'j)hv(ij’) j 
KGa) K Gi 
Li) =K@) 45). (11) 
hv(1'j) 


The terms in braces are of degree \*; the others of degree 2 or lower. We have carried the expressions 
for the L’s far enough to give the expansion of (9) to \*. This means that it has only been necessary 
to carry the off-diagonal elements (11) to the approximation )*. If we were interested in individual 
roots of (9), instead of merely the sum, we would need to carry these elements to the approximation 
\‘, since to determine individual roots of (9) to \4 one would need to know the coefficients of W?, 
W, W® to A8, A”, 16, respectively, whereas we are interested purely in the coefficients of W, W*. 
In virtue of (10), (11), (5), to our approximation the explicit form of (9) is 











mt+L&OiDi-j)-W MC Mm, MLN 
-M,yd I; {QO-W JdTlm J>dIln: i 12) 
M,>m; Jdlm; f(m)-W JXmn; 
with MpUn: JZTIln; JXmn; f(n)-W 
f(Q) = W,°-4V C+ I Dig? + Didi gallet+ Nip] 

H, = W.° —4¥°C, — 4740." +4A47,90,"", J=X2(Ce—Ce) — rn” _ 

T= M147 pDQu'Qo!’ +167 pT Qe"Qs", = Mu = —ND—2"r10,'Qe' 

N=)*r,?Q."(Co—Cz) +047 p°Qo”! : 


*The matrix elements of S~! needed for computing SKS are found by noting that (8) gives nearly a unit matrix, 
So that 


S-1=1—(S—1)+(S—1)?—(S—1)4-++. 








224 J. H. VAN VLECK 


where the i and j summations each embrace a, b, c, d and where 
1/7, = W"+B-—W,°>0, (k=s or p); O,.’=Q.+Ti(Wi°-— W), (k=s or oc) (14) 
(i-j) =1d;+mun;+nn;=cos (C—Hi, C—H,), (z, 7 =a, b, c, d). (15) 


We have treated B as part of the H atoms’ unperturbed energy in defining the 7, or 1/hv(ij). This 
avoids the necessity of treating the B terms as perturbations. The réle of the B terms is thus a 
trivial one, without appreciable influence on the properties of directional valence. When one drops 
terms beyond X, the expansion of (12) is of the form 


(W—W,°)?(W— W,°) +2 (a5 W>+a,W1+ «+ +a0) +A4(b;W7 + + + + +09) =0. (16) 


The coefficients a of the \* terms are easily seen to be independent of the molecular arrangement 
since the sum of the squares of three direction cosines is unity. The terms proportional to W*®, W®, W’ 
exist because the second relation of (14) makes the elements (13) nonlinear in W. The fact that 
(16) is a seventh degree expression in W gives the superficial appearance of a seventh degree secular 
problem, but this is only an illusion, for it is only legitimate to take the four roots which reduce 
to C levels at infinite C—H separations. The other three roots are devoid of any significance, as 
nearness to C levels was assumed in reducing the determinant to the fourth order. The occurrence 
of higher powers than W* in (16) may be obviated by depressing the higher powers of W by the 
approximate relation 


W1=6u® W?—8u®W+3wt+[(6a;w? +03) W?+ (—8a;w? + 6ayw*+a2)W?+--- J 


obtained by solving (16) first with 4\=0 and then correcting to \* inclusive. Here for simplicity 
we have assumed the origin for the energy so chosen that the mean root is zero when \=0, ie., 
that W.°= —3w, W,° =w where 4w is the carbon s-p separation. This assumption involves no loss 
of generality and simplifies the calculations materially. Thus (16) is transformed into a fourth 
degree equation AW!+ BW?+---=0. The sum X of the four relevant roots is thus, 


B d[a3+6a;w? ]+A4[b3+a304+0205;+ w?(6b;+ 12a4a5) — 8w3 (be +a;?) +39b,wt ]+--- a7 
catntiah is . (17 


A 1+Naj+--- 


a = 





It will be noted that the absence of orthogonality makes the procedure somewhat more complicated 
than it would otherwise be. If we were willing to assume perfect orthogonality and so set 7=0 
in (5) there would be no terms higher than W* in (9) or (16) and one could use the usual spur relation 
X = —>-L(ii) which actually is invalidated by the non-orthogonality. It is clearly to be understood 
that the spur in question is that of (9) or (16) and not that of (4), since we are interested in the 
sum of only half the roots of (4). It does not appear a legitimate approximation to assume orthogo- 
nality as (7) shows that T is of the same order of magnitude as Q/hv(C; H) unless fortuitously 
Q and R are very nearly equal. Although the non-orthogonality invalidates the ordinary spur 
relation, it is still very much easier to calculate the sum than individual roots, as the latter would 
require solution of a biquadratic equation, which we have avoided. 

When the explicit values of the coefficients in (16) are calculated by means of (14-15) it is found 
that (17) becomes 

X=athprt+yDidili-j’, (18) 


where a, 8, y are independent of the arrangement of the H atoms, (i-j) is defined as in (15), and 
p is a vector of components }-/;, 5>m;, >-n;. This vector is clearly proportional to the electric moment 
of the molecule. It is particularly noteworthy that to our degree of approximation \‘, the dependence 
of X on the molecular arrangement, i.e., on 1, , «++ ms, is only through the two combinations 
>~>X(i-j)? and p?. The coefficients in (18) have the following values in terms of the integrals (3) 
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B/M! = 2747 pDLQ.0.+ 2(T:Q. = Q; T,)w]+ THT Qs?(Qz ws 4T,w)?+ 7 TeQ0(Qs +4 Tw)? 


+D(t.+ Tp) (Q.7.+ T Qc) + tT pL2Q07 Qe? + 20.7 Q.+4(Q.°T?? ae: Q/T.”)w], 
7/M= $r,°(C,.— Cx) Qo? + 27 p'Qo'+2rp(C, —C,)T,Q.+21,°0,*T, 


(19) 


(20) 


with 7 as in (14). The corresponding formula for a could be given, but is of no interest. 


The all important question is now the signs 
of 8 and y. These are both positive because of 
the following considerations. We have already 
mentioned after Eq. (3) that the C, D, Q, T may 
all be considered positive, and 7 is positive in 
virtue of (14). Furthermore C,—C, is positive 
since the o; wave function overlaps the hydrogen 
atom H; much more than do the 7; functions. 
Hence every term in (20) is positive. Eq. (19) 
contains two negative terms, but these are of 
subordinate importance for the following reasons. 
In the first place, the carbon s-p separation 4w 
may be small compared to the separation 7 
between the H and C levels. If this is the case 
Tw will be small compared to Q, since (7) shows 
that Tis of the order of magnitude Q/7 inasmuch 
as Q and R are usually very roughly of the 
same order. Furthermore, even if the s-p interval 
should be so large that 7w and Q are of the 
same order of magnitude, there are so many 
positive terms that it would be only with a 
very curious and unreasonable collection of 
values for the integrals (3) that the negative 
terms could predominate. In particular, if the 
s and p bonds are nearly equal, the value of 
T.Q, will not differ greatly from that of 7.Q,, 
so that the factors (7,0,—Q.T.) and (Q,?7.? 
—Q,T,?) are comparatively small. 

The most stable molecular arrangement is that 
which minimizes (18), since 2X is the binding 
energy. A convenient alternative way of writing 
(18) is 


X=a+(16/3) y+(6—3y)p? 
+2yd0i, parol: +3. (21) 


Eqs. (18) and (21) are identical inasmuch as 
P= i, i-i=1, and as the sum in (18) included 
jSi in addition to j>i. Now the tetrahedral 
model alone possesses both of the two following 
Properties: (a) absence of an electric moment p; 
(b) a common value —3 for every term of the 
form i-j(ij). (The easiest way to see that the 
tetrahedron fulfills (b) is to note that the 


projections on C—H, of unit vectors directed 
along C— Hs, C—H,., C—Ha4 are all equal, and 
of magnitude —4 because of the absence of an 
electric moment). Eq. (21) shows immediately 
that the tetrahedral model possesses the least 
energy of all models if B=#y. It is not unlikely 
that the latter inequality is fulfilled. For instance, 
if we suppose that the Q,.T are more important 
than the C, D terms, that 7w/Q is small, that 
T,=T, and that the s and p bonds are of 
equal intensity (ie, Q.=Q,, T,=T;.), Eqs. 
(19-20) show that 37 is between one-sixth and 
one-third as large as 6; actually the s bonds 
are weaker than the ~, but perhaps not enough 
to make 37>. 

One can also see beyond all reasonable doubt 
from (18) or (21) that the tetrahedral model is 
superior even if B<?y. It is an immediate 
consequence of (21) that the tetrahedron has 
the least energy of all models devoid of an 
electric moment, since with p=0 the question 
of the sign of B—%y does not enter. Among 
the models with an electric moment, the only 
arrangement that seems at all probable is a 
pyramidal model with the C atom at its apex 
and the four H atoms at the corners of a square 
base. One easily finds that in this model the 
sum )-(i-j)? has its minimum value when each 
C—H axis is inclined at cos~! 3-!=54.7° to the 
axis of the pyramid. Then some of the (i-j) 
(tj) have the value +3, the rest the value —}. 
The 54.7° pyramid thus gives the same value 
for the third term of (18) as the tetrahedron, 
but the pyramid has an electric moment and so 
gives a larger value of the second term. Of course 
this still does not disprove the possibility of some 
freak, unsymmetrical model with an electric 
moment, but we must not forget the physical 
fact that dielectric constant data unquestionably 
reveal that CH, has no moment. Irrespective 
of this fact, one can progress a step further 
towards complete mathematical rigor by noting 
that the tetrahedral model, regardless of the 
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sign of B—?2y4, fulfills the analytical criterion 
DV iLixj](i-j) =0 (22) 


that >-;;(i-j)* be an extremum. That (22) is the 
condition for the minimum of this sum can be 
verified by writing the sum explicitly in terms 
of the /’s, m’s, and n’s, and then differentiating 
with respect to, say, the m’s after the n’s have 
been eliminated by the relations /;?+m?+n7=1. 
Satisfaction of (22) by the tetrahedron shows 
that it is stable with respect to adjacent arrange- 
ments, but this does not say that there are not 
radically different arrangements also satisfying 
(22) and conceivably giving still lower energy if 
there is an electric moment. The 54.7° pyramid, 
the square coplanar model, and collinear arrange- 
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ments i-j = +1 are all solutions of (22), but each 
gives a higher value of (18) than the tetrahedron. 

Case W©—W*">@Q. This is another limiting 
ionic case, differing from the preceding in that 
the drift of electronic charge is from C to H 
rather than from H to C, making the model 
nearly C*+(H,)*-. This case can be treated by 
the same kind of perturbation development as 
the preceding one, but we shall omit the calcu- 
lations. Without appreciable computation, one 
can show that the tetrahedral model has the 
least energy if one neglects the corrections for 
non-orthogonality, so that one can use the spur 
theorem.® Of course the corrections are im- 
portant, but in the preceding case, at least, they 
reinforced rather than counteracted the other 
terms. 





3. ANOTHER METHOD OF APPROXIMATION WITH THE HUND-MULLIKEN METHOD 


When Q is of the same order of magnitude as, or larger than, the separation W” — W° between 
the C and H levels, the molecule is no longer usually highly polar, and it is naturally no longer 
legitimate to use the perturbation developments employed in Section 2. Because of this difficulty 
it appears no longer possible to find the binding energy for arbitrary arrangements. Instead it is 
necessary to restrict our attention to models with some degree of symmetry so as to permit factori- 
zation of the secular determinant (4). Fortunately the symmetrical models are, of course, the most 
likely. Three types of symmetrical models will be described and calculated below under (a), (0), (¢), 
and it will be shown that of these the tetrahedron is the most stable. Even with the symmetrical 
types the calculation is too complicated to be tractable’? without the following three simplifications: 
(I) neglect of the s-p separation in carbon; (II) neglect of the C and D terms in (4-5); (III) neglect 
of terms arising from want of orthogonality (the 7 terms in 4-5). 

Assumption (I), viz., that W,°— W,° is very small compared to Q, means conditions such that 
s-p hybridization would be a good approximation were the calculation to be made by the Heitler- 
London method. It is not necessary to make the more stringent restriction that |W” —W°|<Q, 
as the latter proves to permit little additional simplification beyond that arising from W,°= W.° 
= W°. We can clearly suppose that W” — W° is at most of the order Q, since if W“— W°>Q we 
would merely revert to the case already covered by perturbation theory. In consequence of this 
fact, we shall show in the next paragraph that (II) and (III) are good approximations if (I) is granted. 

Assumption (II) means that the molecular bonds are supposed to arise mainly from the Hund 
resonance (our Q) terms rather than from penetration (our C, D) terms. In the preceding cases 


‘If we set T=0, the root sum Y which interests us 
when W°—W“>>0Q is Z—X where Z is the spur of the 
complete matrix (4) and X is our previous root sum. 
The value of Z is independent of the angles, and hence Y 
has an extremum when X does. When W°—W#>>(Q the 
tetrahedral arrangements gives a minimum of Z—X and 
a maximum of X rather than the reverse as previously. 
This is because sign changes in 7 and in certain of the 
integrals (3), occasioned by the electronic drift being from 
C to H rather than from H to C, make the expressions 
B, y defined in (19, 20) become negative. 





7 For the regular tetrahedron, which is a special case of 
the models of class (a) or (c), and which has the most 
symmetry of all, none of the approximations (I), (II), 
(III) are necessary, and an exact expression can be derived 
for the energy, which was given in Eq. (4), Part I. The 
constants employed in Part I have the following signifi- 
cance in terms of the integrals (3), etc. 


a+3b=4C,4+ W.°, a—b=4C,+3C,+W;,°, 
6+3e=2(-Q.+T.W.°), c—e=—(4)! (Q.—TeWs°), 
d+3f=—2T., d—f=—(4)'Ts, g~0,h~0, R=W"+8. 
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we stated that it was not allowable to forget the C, D terms, even though C, DQ, inasmuch as 
they affected the dependence of energy on arrangement in the approximation CQ?, DQ?, whereas Q 
alone appeared only in the fourth order (cf. Eqs. (19-20)). In the present case, however, the develop- 
ment in ascending powers of Q/hvy(H;C) is no longer used and in consequence the Q terms will 
influence the energy of spatial arrangement in the first rather than fourth approximation, thus 
surely becoming more important than the C terms. This can be seen concretely from later specific 
energy formulas; the situation is roughly similar to that of an expression like (a?+0? cos? @)! 
+(a?+0? sin? 0)', which depends on angle in the approximation 5‘/a* if b°<a?, but in terms of the 
order 6 otherwise. Likewise neglect of orthogonality is a much better approximation than when 
W" —W°>@Q. This can be seen from dimensional considerations. The expressions 7,, 7’, are dimen- 
sionless and small compared to unity, while Q has the dimensions of energy. Terms of the form 
Q*(W"—W°)"T" can thus appear along with those of the form Q*+¥ and could be coordinate in 
importance with the latter in the previous cases where |W” —W°| >Q but not in the present one, 
where W” — W° is at most comparable with Q. 

(a) One type of model which is readily treated and which we shall call class (a) is obtained by 
the following scheme. Let the coordinates of the C atom be 0, 0, 0 and those of the four other atoms 
be x, 0,2; —x, 0,2; 0, y’, 2’; 0, —y’, 2’; Let w=tan-!2/x and w’=tan— 2’/y’. This class includes the 
regular tetrahedron as the special case w = — w’ = tan~! 2—} = 35.3° and the Guillemin-Henri pyramidal 
model (H atoms at corners of square base, C at vertex) as the special case w=w’, z=2’. In order 
to make our calculations applicable to molecules of the form CH2Xz as well as just CH,, we shall 
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Fic. 1 (class (a)). For CH, or CH2X2. Here Oa=Ob; Fic. 2 (class (c)). For CHy, CH3X, or NH;. Here Ob = Oc 
Oc=Od, with H atoms at c, d and a C atom at O. The =Od and ab=ac=ad with H atoms at b, c, d; a C (or N) 
atoms at a, b are either both H or both X. The regular atom at O and an H or X (or no) atom at a. The g axis is 
tetrahedron is the special case w= —w’ =tan-! 2-$=35.3°. directed along Oa and the xz plane contains a0b. The 
n the text the notation 9 or P is used for integrals in- regular tetrahedron is the special case @ =cos~!— } = 109.5°. 
volving either a or 6; and N or Q for those involving cord. The integrals N or P apply to a; and N or Q to d, ¢, d. 
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not assume that the atoms in the xz plane are the same as those in the yz plane.® It will not be 
necessary to suppose? x?-+2?=y"+2”, i.e., that the two kinds of atoms are equidistant from the 
C atom, but unfortunately it must be assumed for tractability that they all have the same energy 
W*. The resonance integrals between the C atom and the atoms in the xz plane will be denoted 
by the letter P rather than Q in order to exhibit the fact that they are held by bonds of different 
intensity than those in the yz plane. It is convenient to choose the eight wave functions as follows 


[y(2s), 2(VatyotvetYa), ¥(202), 3(WatYe—ve—Wa) ], 
[y(20.), 2-3(Ya— Yo) J, [y(20,), 2-4(We—wa) J. 


with notation as in Eq. (1) except for obvious contractions. The secular determinant, originally 
of degree 8 (cf. Eq. (4)) then factors into a biquadratic and two quadratics. The biquadratic is a 
quadratic in (W— W*")(W—W‘) and so easily soluble. Wave functions bracketed together belong 
to the same factor. The factorization is readily verified and explicit form of the factorized determi- 
nants constructed from the transformation property (6) of the carbon p wave functions under 
rotation. Four of the roots are as follows: 


(23) 


W,=F(P.? cos? w), W2= F(Q,? cos? w’), (24) 
Wa. =Ga(P,? sin? w+ Q,? sin? w’+P2+0,?; 2P.Q0, sin w—2P,Q, sin w’), (25) 
where F(a) =3(W#-+W*) —4y[(W¥—W?)+ 8a], as 


Ga(b; c) =} (W727 + W*) —3V0(W? — W°)?+4044(0? —c?)*). 


The radical sign is used throughout to denote the positive square root. The other four roots are 
similar except that the factor before the radical sign is +} rather than —}. These other roots are 
“anti-bonding,”’ i.e., give higher energy than the lesser of the two quantities W”, W°, and hence 
are of no interest to us. The Pauli principle allows two electrons to each bonding root, and so the 
molecule’s energy is 

We first note that for given absolute values of w, w’, and nonvanishing P, Q, the expression (27) 
has a lower value when w and w’ have opposite rather than similar signs. This follows from (24-26) 
and the inequality 


G,(b; c’) + G_(b; c’) <G,(b;c)+G_(b;c) (b2c’>c>0). 


Hence for any molecule of the type CH, or CH2X_ the pyramidal model is never that of least energy. 

Now let us make the specialization Q=P appropriate to CH, if all four H atoms are equidistant 
from C, and plot (27) as a function of the variables w, w’. The value of (27) is then symmetric with 
respect to the line w= —w’ in the w, w’ plane, as well as with respect to w=w’. Hence if we consider 
the behavior of (27) along any curve in the w, w’ plane, (27) will have extremums where the curve 
crosses either of these lines. Furthermore, an extremum of (27) along the line w= —w’ occurs at 
w=tan-! 2-3, as can be verified by differentiation of (27) with respect to w after first setting w’ = —o. 
Hence the regular tetrahedron gives an extremum of the energy, which closer examination shows 


7 
6 


8 Most commonly the atoms X will be halides, in which | (3) obviously are to be calculated for an mpo rather than 
case they will be in ?P rather than in 2S states like hydro- | 1s state of X. Also no harm is done if the letter H similarly 





gen. Despite this difference, the analysis is still appli- 
cable if we suppose that in the halide’s incomplete shell 
ms? mpx* mpo, the axis of mpo is directed towards the 
central atom C and if we neglect the bonding power of ms 
and mpzx. Such a choice of the axis of mpo is clearly 
demanded for overlapping favorable to bonding. With 
these assumptions, only the 2p0 electron of the halide is 
involved in the secular problem (4), and the integrals 





symbolize some atom other than hydrogen and not 
necessarily in an S state. Thus the conclusions of out 
calculations also apply to compounds such as CCls 
CCI,Br2, CCl;Br, as well as CH4, CH2X2, CH3X. 

The assumption x?+2?=x"+2” is, however, always 
needed for the simplification P=Q which we make in 
connection with CH,. 
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to be a minimum rather than maximum. The tetrahedral model is thus stable with respect to adjacent 
models. This is not quite the same as demonstrating rigorously that it has the least energy of all 
models of the class, as conceivably there might be a lower minimum at some other values of wo, w’. 
Because of algebraic difficulties, it would be difficult to solve the equations 0X /dw=dX/dw’ =0 
explicitly for w, w’, but rough inspection of the behavior of (27) for various values of the arguments 
shows that the occurrence of a lower minimum is exceedingly unlikely. In fact the only other mini- 
mum appears to be the 54.7° pyramid (w=w’ =35.3°), which is a minimum along the line w=’, 
but is a higher minimum because of the considerations in the preceding paragraph. It is interesting 
to note that this pyramid and the tetrahedron are minima for all values of the relative intensity 
Q./Q. of the s and p bonds, including Q,/Q,=0 and excluding only the trivial case Q,/Q,= © where 
there is no directional valence. The existence of a nonvanishing Q, is necessary only to make the 
two minima unequal, as the signs of w, w’ are immaterial in (24-25-26) if Q,=P,=0. Hence just 
as in Section 2 the tetrahedron takes precedence over the 54.7° pyramid only because of the interplay 
between the s and ~ bonds. 

(b) Plane models. Consider a plane model formed by placing the C atom at the center of a 
parallelogram, an atom X at each of two opposite corners, and an atom H at each of the two other 
corners. In this type of model the energy obviously has the same value when the angle a between 
the diagonals of the parallelogram has the value a as when it has the value r—a. Hence a= 7/2 
gives an extremum, easily seen to be a minimum. With such a value of a, this plane model is merely 
the special case w = w’ = 0 of the models of class (a) and hence has higher energy than the tetrahedron. 

Another readily soluble model in which all atoms are in the same plane is that in which the X 
atoms are at the upper corners, the H atoms at the lower corners, of an isosceles trapezoid, with C 
on the axis of symmetry. This model is obtainable from class (a) by rotating the plane labelled dOc 
in Fig. 1 by 90° around the z axis. We can thus take the coordinates of the atoms to be the same 
as in class (a) except that those of the last two atoms become x’, 0, 2’; —x’, 0, 2’ and w’ =tan- 2’ /x’. 
The selection of wave functions and factorization are the same as in (23) except that ¥(2c,) gives 
an isolated root W°, while the other three wave functions in the second line of (23) give a cubic of 
which one root is W”. Two of the roots are still given by (25) while (24) is replaced by 


W,' = F(0), W,! = F(P,? cos? w+ Q,? cos? w’). (28) 


Using (26) one finds W,’+ W2’ > W,+ Wz, and so the plane model is less firmly bound than the spatial 
model of the same w, w’. 

(c) Another tractable case, which we shall call class (c) and which is illustrated in Fig. 2, is that 
in which the C atom is at the center of a polar coordinate system 7, 6, ¢ while the coordinates of 
the other four atoms are r, 0, 0; 7’, 6, 0; r’, 0, 37; 7’, 0, $x. In order to make our calculations appli- 
cable to CH;X, we shall not necessarily assume that the atom a on the polar (z) axis is identical 
with the other three atoms held by the C atom, which we may take to be identical H atoms. As 
before, resonance integrals connected with the C—X bond are denoted by the letter P. Only in 
the particular case CH, shall we suppose r=r’. The regular tetrahedron is included as the special 
case 9=cos—! (— $) =109.5°. The choice of wave functions and factorization are indicated by the 


following scheme: 
[¥(2s), ¥(202), Ya, SF4'(Yot+ye+a) J, 
[¥(2.2), (3) bh 2Ve— 2a) 1, [¥(20,), 2-*(Ye—Ya) ] 
and the bonding roots are 


W, = W.= F(3Q,? sin? 6), W.. =G..(3P.22+30,? cos? 0+ 3P,2+30,; 3! cos 0P,0,—3'P.Q,) (30) 


(29) 


with notation as in (26). One verifies by differentiation with respect to 0, that with (30) and the 


specialization P=(Q appropriate to CH,, the root sum (27) has a minimum at the tetrahedral angle 
9=109.5°, 
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4. Case W"—W,° But not W” — W,° LARGE COMPARED TO Q IN HUND-MULLIKEN METHOD 


This case is a mixture of those in Sections 2 
and 3, and is perhaps a better approximation 
than either to reality since the s-p separation in 
C is not small, and is probably greater than 
W*"—W,°. Here the 2s state of C is assumed 
to lie much deeper than either W” or the 2p 
states of C. One twice-filled root of the secular 
determinant is then substantially a carbon 2s 
level; in other words by far the largest coefficient 
in the corresponding linear combination of wave 
functions is that of the carbon 2s one. If W” 
= W,°, and if we neglect the interaction between 
s and p bonds, which is small in the present 
case, the six electrons not absorbed by 2s will 
in the mean be distributed equally between the 
totality of the four H atoms and the three p 
levels of C. Thus the 2s but not the 2 part of 
the problem will be amenable to the near-ionic 
approximation and the molecule will be sub- 
stantially of the structure C+(H,4)~, except inso- 
far as W¥#W,°. As a first approximation, the 
bonding effect of the 2s wave functions of C can 
be neglected, since this could be treated by 
proper adaptation of the ‘‘near-ionic’’ perturba- 
tion method of Section 2, and so will give rise 
only to terms of order Q?/hv(H; C,) and higher. 
Such terms are small compared to the p bonding 
energies, which are of the order Q. A further, 
less important effect enhancing the importance 
of the » as compared to the s bonds is the fact 
that the 2s wave functions represent firmer 
binding in the C atom and so less “‘overlapping”’ 


of the H atoms than do the 2 functions, thus 
making Q,>Q,. As a first approximation the 
bonding action of the 2p states may be treated 
by the method of Section 3 with Q,=P,=0. 
There is now, to be sure, the difference that 
there are six rather than eight electrons in the 
secular problem, since two electrons have been 
absorbed and isolated by the deep 2s level. 
This difference, however, is immaterial inasmuch 
as in each example of Section 3 one root W_ 
becomes nonbonding when Q,=P,=0, and we 
have merely to eject 2W_, which is now inde- 
pendent of angular arrangements, from the sum 
(27). The analysis of Section 3 can thus be used 
and shows that the two most stable models are 
the regular tetrahedron and the 54.7° pyramid. 
We saw that with 0,=P,=0 these two models 
were on a par. However, correction for the fact 
that the s-p interplay (not to mention H-H 
repulsions) is not completely negligible will make 
the tetrahedron have the lower energy of the 
two. This could be demonstrated by introducing 
the interaction between the s and p bonds as a 
perturbation of order Q?/hvy(H;C,) and above. 
We shall, however, omit the proof, since most 
readers will probably consider the result suffi- 
ciently obvious in view of the fact that the 
present case is intermediate between the cases 
of Sections 2 and 3, and in both these cases we 
showed that the s-p interplay favored the regular 
tetrahedron over the 54.7° pyramid. 


5. CALCULATION BY HEITLER-LONDON-PAULING-SLATER METHOD OF ELECTRON PAIRS 


It unfortunately appears impossible without undue labor to make calculations by the Heitler- 
London method for anything like as great a variety of arrangements as in the preceding Hund- 
Mulliken procedure. The difficulty with the H-L method is that the secular determinant is of degree 
14 and hence is too complicated for explicit solution without either the assumption of electron pairing, 
wherein the members of a pair have their spins anti-parallel, or else the Heitler-Rumer assumptions 
to be discussed in Part III. A pair, of course, represents a saturated chemical bond, and the total 
wave function is invariant when the purely orbital coordinates of the members of a pair are per- 
muted. The qualitative idea of electron pairs as units of chemical structure is due originally to 
Lewis,'° who proposed it before the days of quantum mechanics and hence without the wave- 


mechanical interpretation. 


Since the carbon p wave functions are proportional to x, y, z, and so transform like vectors, 
the most general real linear combination of the carbon s and p wave functions is of the form 


0G. N. Lewis, Valence and the Structure of Atoms and Molecules, Chem. Cat. Co., 1923. 
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¥i° =aw(2s) +B (203). (31) 


Clearly the H atom to be paired with this wave function should be located on the axis of quantization 
C—i of the p term in (31), as then the overlapping is the largest possible. This concept is the basic 
idea of the Slater-Pauling theory of directed valence. Unfortunately it is impossible to construct 
four wave functions of the form (31) which are normalized and mutually orthogonal and which 
project respectively in four arbitrary directions i=a, b,c, d. This is true since the requirements of 
normalization and orthogonality yield 10 relations while there are only eight coefficients a, ---, Bs 
at our disposal. The need of orthogonality is vital, since two electrons with non-orthogonal wave 
functions have a constraint on the angle between their spins and so are not free to make valence 
bonds with electrons in other atoms. (The easiest way of seeing this is to pass to the limiting case 
in which the non-orthogonal functions are identical: the spins associated with the two functions 
are then anti-parallel and the valences of the two electrons are hence saturated intra-atomically. 
Abandonment of normalization would not increase the number of allowable directions, as it merely 
would introduce a trivial constant factor.) 

The only models, other than highly unsymmetrical ones, for which it appears possible to construct 
wave functions with the necessary directive and orthogonality properties are those of class (c), 
Seciion 3, for which cot? @=3 and those of class (a) for which 


cos 2w cos 2w’=sin w* sin w” with |w’|=47, |w|=3r, w’/w<0. (32) 


The highly unsymmetrical models presumably would be needed for say CWXYZ, but are quite 
unlikely for CH,, CH3X, or CH2X2, so that we shall not consider them. For class (c) the proper 
set of orthogonal functions is 


Wa° = (1—2 cot? 6)4¥(20.) —2! cot 6 (2s), 


33 
¥i° = (3 cosec? 0)4¥(20;)+(3—§ cot? @)4y(2s) (j=b, c,d), 
while for (a) 
¥;° =(1+cos 2w)—![y(20;) + (cos 2w)4¥(2s)] (j=a, b), 


(34) 
¥i° =(1+cos 2’) ![¥(20;) + (cos 2’)!y(2s)] (j=c, d) 


with notation as in Section 3. The orthogonality is readily established with the aid of the relation 
SS S'¥(2¢;)(20;)dv=cos (i, j) after the cosine is expressed in terms of @ or of w, w’ with the re- 
striction (32). When cos @= — 3 or tan w= —tan w’ = 2-! the models reduce to the regular tetrahedron, 
and then (33) or (34) become the usual Slater-Pauling sp* functions. 

Because of the inability to find many symmetrical models which permit natural C—H paired 
bonds, it is in a certain sense impossible to demonstrate in completely quantitative fashion that 
the tetrahedral model is the most stable in the Heitler-London method. Pauling’s criteria of ‘‘maxi- 
mum overlapping,” to be sure, make this model very reasonable indeed, as he has shown." Quite 
irrespective of this fact, the mere inability to “localize” the bonds with symmetrical models other 
than those cited must be regarded as strong evidence that other models do not fit naturally into 
the scheme of things and so yield less bonding energy.” 


cunniiaestninensiiemes 


"L. Pauling, J. Am. Chem. Soc. 53, 1367 (1931). required is the same as the number of different kinds of 

“Essentially this point has also been noted in an | atoms attached by C. If this is the case then all the C 
interesting paper by R. Hultgren, Phys. Rev. 40, 891 | wave functions of CH, have the same ratio and the 
(1932). His wave functions are more general than ours in | tetrahedral model follows uniquely simply from the 
that d states are admitted to the linear combinations, but requirement of orthogonality. Nevertheless it seems of 
less general in that all the functions are supposed to have | some interest for us to verify that introduction of two 
the same hybridization ratio. Eqs. (33) and (34) each | unequal ratios in place of one raises the value of (35) 
involve two such ratios. Four would be required for | in CHy. 
CWXYZ, Presumably the number of different ratios 
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There is no difficulty in showing that for CH, the regular tetrahedron has the least energy of all 
the models of class (c) with cot? 6=3 or of class (a) obeying (32). Under the assumption of perfect 
pairing and neglect of everything but exchange terms the molecular energy is 


W=Dd:Ju- 4). 5, ii) I ij +constant (35) 


where J;; is the exchange integral for the combination y;”, ¥;°. Eq. (35) is the embodiment of the 
well-known results that the exchange energy associated with a saturated bond a, b, is +J.a and 
that this energy is —}J,- between a member 0 of a pair and any other electron c not its mate in 
the pair. As usual, we shall neglect all H-H forces or H-H overlapping. By means of the transfor- 
mation relation (6) it is possible to express (35) in terms of 6 or w, w’ and the essentially diatomic 
integrals 


Niu= —f fff ff vmcrrvscayetytyyn(2rdordvs (k, l=2s, 2po, or 2pm of C). (36) 


The minus sign has been inserted in order to make most of the N,: positive. Except for sign, the 
expressions (36) are exchange integrals when k=/. The value of N,, is zero. It is readily shown, from 
band spectrum data or by approximate calculation, that N,, and N,, are both positive, while N,, is 
negative. 

In the case of (33) we find 


W = —$Neet4Nart3 Nes t23(1—2 cot? 6)}(—3 cosec 6+3 cot 6)Nos— (Neo — Noe) (1—3 cot? 6) 
+ (Nes — New) (3 —3 cot? @) +N e— Noe t23(3Ns —3Nos) cot a(1 —2 cot? 6)}. (37) 


In order to make our calculations adaptable to CH;X or CHX;, we have allowed atom a in Fig. 2, 
Section 3, to be different from the other three atoms }, c,d held by C, and so yield values of the 
integrals (36) which may be different than those for the other bonds and which are distinguished 
by being printed in German type. Thus in CH;X, 9? denotes (36) with subscript X substituted for 
H, while {= N in dealing with CH, When one sets Jt=N, (37) depends on @ only in the terms 
proportional to N,., which are peculiar to s-p hybridization. One verifies by differentiation with 
respect to @ that then if V,,>0, Jt=N, Eq. (37) has a minimum at cos 6= —}, giving us the regular 
tetrahedron. The assumption N,,>0 is correct, for the charge cloud of a 2p0 wave function is 
dumb-bell shaped with opposite signs for ¥(2c) on the two ends, so that the sign of N,, is dependent 
on which end of the dumb-bell one places the H atom. To obtain favorable bonding one must choose 
the end for which N,, is positive, as is obvious from ‘‘overlapping”’ considerations or from the fact 
that the minimum of (37) is higher with N,,<0 than with N,,>0. When cot? 6=3 the three axes 
C—b, C—c, C—d are mutually orthogonal; there is then no hybridization and atom a is held only 
by an s bond; under these conditions it is not necessary to restrict atom a to the axis @=0, as it 
is easily shown that the energy is independent of the alignment of C—a relative to the three or- 
thogonal axes C—b, C—c, C—d. 
In the case of models of class (a) with the restriction (32) one finds that (35) becomes 


W = —$Nect4Nart4NestNosg(w’, w) + Neosg(w, wo”) + (Nes— Nos) [1 —2f(w’) +f(w) J 
4- (Nee — Nee) [ —2f(w) +cos? 2uf(w) +2 sin? w sin? w’f(w’)] 
+(Nr2—N-rx)[sin? 2wf(w) +2f(w’)(1—sin? w sin? w’)] (38) 
f(w) =(1+cos 2w)—, g(w’, w) = —2(cos 2w)§[1+f(w) ]+4(cos 2w’)'f(w’) sin w sin w. 


with 


In order to make the calculations applicable to CH2X2, we have allowed atoms a and b to be different 


Cf, for instance J. C. Slater, Phys. Rev. 34, 1307 13 See, for example the table given by Coolidge, Phys. 
(1929). His K is our J. Rev. 42, 198 (1932). 
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from c, d, and so give rise to integrals Jt rather than N. As before, when we make the specialization 
t= N appropriate to CH,, a dependence of W on angle exists only in virtue of the terms involving 
N,s. That (38) then has an extremum at the tetrahedral configuration is most easily demonstrated 
by a symmetry argument: viz., when J(=N the expression (38), regarded as a function of w, w’, 
is symmetric about the line w= —w’ in the w, w’ plane, and so will have an extremum along the 
curve (32) at the intersection of this curve with this line. The coordinates of the intersection are 
w= —w’=tan-! 2-}, the tetrahedral values. 

Slater“ has remarked that the Heitler-London and Hund-Mulliken methods are equivalent in 
the last analysis as they reach the same goal when pushed to final accuracy by inclusion of higher 
approximations (notably addition of polar terms in the H-L and antisymmetrization and inclusion 
of electron repulsion in the H-M). Of course what Slater says is ultimately true, but when one 
stops at the first approximation the resemblance between the two methods is very remote in poly- 
atomic molecules. The superiority of the tetrahedron over other molecules which could be calculated 
was shown above to owe its existence in the H-L-P-S method entirely to the integral N,, which 
represents the interplay between the s and p bonds and epitomizes the s-p hybridization. On the 
other hand in the H-M method the bonds alone were sufficient to insure the superiority of the 
tetrahedron over all other models with the single exception of the 54.7° pyramid. The tetrahedron 
was superior to this pyramid because, unlike the latter, the tetrahedron involved no s-p hybridization 
at all! 


6. CONFIGURATION OF SOME OTHER MOLECULES, NotaBLy CHX;, CH2Xs, CWXYZ, BeX2, OX, 


The ‘“‘near-ionic’’ perturbation calculation 
given in Section 2 can be extended without 
particular difficulty to the case in which the 
central atom is not necessarily carbon, and 
attaches any number of other atoms, provided 
only that in all there are precisely eight electrons 
available for the bonding. There are, of course, 
in general, plenty of electrons of lower principal 
quantum number, but these are treated as 
nonbonding. In place of (18) one finds 


X=a'+ Vid Guli- J+ViVivii-j)? (39) 
where i, j are unit vectors along the axes C—i, 
C-—j and where £;;, y:; are defined nearly as in 


(19), (20) except that they are symmetrical in 
i,j. For instance 
D*r43Q4'Tp'Qo? +D'75'Qs'tp'Qo* 

replaces 27,7,DQ,Q, in the first term of (19). 
Superscripts are, of course, now added to 
indicate that the integrals (3) have different 
values for a bond involving atom 7 than one 
involving atom j, and instead of (14) one has 
1/74= Wi+B;—W,° with W°, W’ denoting re- 
spectively the unperturbed energy of the central 
atom and that of atom j. Our discussion does 
not apply to double or multiple bonds, and so 


Ginemmeideamiet 


“J. C. Slater, Phys. Rev. 41, 255 (1932). 


we must suppose that atoms i, j, --+ each have 
only one solution of Schrédinger’s equation 
which is important for bonding. 

Eq. (39) throws considerable light on molecular 
arrangement. When we remember that the £;; 
and y;; are positive and that the terms y;; are 
due purely to the bonding by the carbon p wave 
functions, while the 8;; are due to the interplay 
between the s and p bonding actions (cf. Eqs. 
(19, 20)), we see that the axes C—7i, C—j joining 
any two atoms 7,7 to the central atom seek to 
align themselves at right angles to each other if 
the » wave functions alone are important for 
bonding. However, if the s-p interplay is ap- 
preciable, the angle i—C —j will tend to increase 
towards 180°, as i-j has its minimum value at 
180° as compared with 90° for (i-j)?(i#j). Of 
course it is impossible for all bonds to be mutually 
at right angles if more than three atoms are 
attached to the central atom, and in CH, the 
tetrahedral arrangement represents, so to speak, 
the nearest the bonds can come to such angles. 
In molecules in which the attached atoms are 
not all alike, the values of 8;;, y:; will be greater 
for certain pairs of bonds than for others. These 
pairs will tend to set themselves at the optimum 
angles even though this involves making some 
of the other axes inclined to each other at 
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unfavorable angles. For instance, consider a 
hypothetical ideal case in which y12, 13, Y23, Y14; 
‘Ye4 are infinitely large compared with ys4--- 
and with any 8. Then atoms 1, 2, 3 would be 
located along three mutually orthogonal axes 
with origin at the central atom, and atom 4 
would be located on the axis C—3 (probably on 
the opposite side of C from 3) though this leads 
to the worst possible value of the term involving 
Bi2, Bis, B23, Bis, Bos, ss In short, the most 
important pairs of bonds are greedy and tend to 
assume the most favorable angles at the expense 
of their fellows. Hence the correct structure 
when the carbon atom attaches four atoms not 
all of the same type may be described as a 
distorted tetrahedron. (Strictly speaking the 
model is still a tetrahedron, but not a regular 
one. We insert the adjective ‘‘distorted”’ because 
we use the word tetrahedron usually only in 
connection with the regular type.) One can 
easily see that pyramidal models, where all four 
attached atoms are on one side of a plane 
containing the central atom, and plane models 
will yield higher values of (39) than a properly 
chosen distorted tetrahedron. Some particular 
molecules are the following: 

CH2X»2. Here the model will be one of type 
(a), shown in Fig. 1, Section 3, but with in 
general — w’~w+tan-! 2-! unless the C—H and 
C—X bonds are very nearly equal. If P, is large 
compared with P,, Q,, and Q,, the angle w will 
be nearly 45° and w’ nearly 0. 

CH;X or CHX;. The appropriate model is one 
of type (c), shown in Fig. 2, with the tetrahedral 
value 109.5° for @ in CH;X only if Buy =Bun, 
Yux = Yun: If Yuu? Yux and y>8 then @ will be 
nearly tan-! 2-}, and the three C—H axes will 
be nearly at right angles. If yyx>vyy and y>8 
the value of @ will be nearly 90°, and the plane 
of the three H atoms will almost contain the 
central atom. 

NH. For mathematical purposes the ammonia 
molecule is equivalent to a special case of CH;X 
in which atom X is held by a bond of zero 
intensity (i.e., Bax = Yux =0) and so is effectively 
non-existent. As already stated, the fact that 
the central atom is N rather than C is immaterial 
for formal considerations so long as there are 
precisely eight electrons. The value of @ is found 


by minimization of (39) to be cos [(y—8)/3y]} 
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if B=y and 90° if B=y. Thus @ is bounded 
by 54.7° and 90°. If B=0 is a good approxi- 
mation, 6 will be 54.7° and the three NH axes 
will be inclined at 90°. If the s-p interplay 6 
were the most important thing, which is surely 
not so, the angle @ would be 90°, and the model 
would be coplanar. Actually, Lueg and Hedfeld 
determine @ from band spectra to be 73°, while 
Dennison and Uhlenbeck find @=68° by a 
different method but likewise based in part on 
band spectra.” Either of these values is within 
the critical range predicted by the theory. 

CWXYZ. Here, naturally, none of the angles 
should in general be equal. 

All the above conclusions are seemingly based 
on (39) and hence on a “near-ionic’”’ approxi- 
mation which is not usually an especially close 
one. To remove this objection we have therefore 
given another method of approximation for 
CHeX. and CH;X (or CHX;) in Section 3, 
which applies under different conditions, viz., 
items (I), (II), (III) of Section 3, and the 
assumption W”=W*. These other conditions 
are more likely to be met when the molecule is 
primarily nonpolar, and agree with (39) in 
predicting that the angles cease to be exactly 
tetrahedral in CH2X»2 and CH;X or CHX;, as 
(27) ceases to have a minimum at —w’=w 
=tan 2-} (with Eqs. (24-25)) or at cos @= —} 
(with Eqs. (30)) when P#Q. The distorted 
tetrahedra were again found in Section 3 to be 
better than the corresponding pyramidal and 
plane models. Similar results can also be obtained 
with the intermediate method of Section 4, 
which uses the near-ionic approximation only 
for the 2s part of the problem. The same con- 
clusions thus hold in our three methods of 
approximation based on the Hund-Mulliken 
viewpoint. Additional confidence is supplied by 
the fact that in the Heitler-London-Pauling- 
Slater method the tetrahedral angles disappear 
as soon as the calculations are adapted to 
CH2X», CH;X or CHX;, since (37) and (38) 
cease to have minima at exactly the above 
values of 29 or w, w’ if NAN. 

16 Lueg and Hedfeld, Zeits. f. Physik 75, 559 (1932); 
D. M. Dennison and G. E. Uhlenbeck, Phys. Rev. 41, 
313 (1932) or N. Rosen and P. M. Morse, ibid. 42, 210 
(1932); these determinations are presumably more accu- 
rate than older spectroscopic work which gave values of 4 
in the neighborhood of 60°. 





STRUCTURE OF CH, 


That the tetrahedral angles are not preserved 
in CH2X» and CHX; is in agreement with x-ray 
data’ which reveal the following values for the 
angles between the two C—Cl axes: CHCl, 
116.4°; CHeCle, 123.8°, whereas the value dis- 
closed for CCl, is 109.5°, exactly the tetrahedral 
value. The deviations from 109.5° in CHCl; and 
CH2Cle are in the direction to be expected from 
our theory. Namely, it is known experimentally 
that C—H bonds (94 kg cal./g mol.) are firmer 
than C—Cl ones (72 kg cal./g mol.). Hence if 
the p-bonding effects overshadow the s ones, as 
they doubtless do, both the H-M and H-L-P-S 
methods predict that the angles H—C—H and 
Cl—C—H will become nearer 90° than in the 
regular tetrahedron, but this entails making the 
angle Cl—C—Cl greater than 109.5°. It is clear 
that the deviation should be smaller in CHCl; 
than in CHCl. since in the latter there are more 
H—C-—H or H—C—CIl angles, causing more 
overpowering of Cl—C—Cl directional valency. 
If the latter and s-bonding were negligible, the 
angle Cl—C—Cl would be 120° in CHCl; and 
180° in CH2Cle. 


The same general conclusions on NH; as before with 
(39) are also obtained if one uses the other methods. 
With the method of Section 3 one minimizes (27) with 
P=0 in (30). The resulting value of 0 is given by 


cos? = 3 — 3(Q.?/Qo*) 


and 6=z otherwise. The method of Section 4 predicts 6 
to be slightly greater than 54.7°=cos~(3)4. With the 
Slater-Pauling method one uses the wave functions (33) 
but assigns two electrons to vc, Since nitrogen has five 
electrons. The resulting expression for the energy is (37) 
with N=0 but augmented by a contribution 


if Q?=3Q.? 


(40) 


:(1—2 cot? @) cos? 0Nog ; 
—3x24(1—2 cot? 6)! cot 6 cos 0Nos 


+2 cot? 6N..+(1—2 cot? 6)(W,—W,) (41) 


due to the fifth electron. The minimum can only be found 
numerically, but is located at a somewhat greater angle 
than 54.7°, The term proportional to W,—W, has been 
inserted in (41) to make some allowance for the influence 
of the finite s-p separation of the central atom upon the 
variation of the energy with angle. That the s-p separation 
does give an angular term is obvious from the fact that 
with the wave functions now being used, the N atom is 
in the configuration sp at @=54.7° and in the configuration 
sp‘ at @=90°, so that the internal energy of the N atom 


° L. Bewilogua, Phys. Zeits. 32, 265 (1931); P. Debye, 
Zeits. f. Electrochemie 36, 612 (1930). 
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is raised in passing from 54.7° to 90°. When we use notation 
such as sp* here we do not mean that individual electrons 
are of the 2s or 2p types but that the sum of the squares 
of the 2s and 2 parts of the individual functions add to 
1 and 4, respectively. The specific form of the term 
proportional to W, — W, in (41) is obtained by examination 
of the sum of the squares of the 2s and 2 parts for various 
angles. The corresponding term was omitted in our 
formulas (39), (38) for CH4, CH2X»2, CH;X; its inclusion 
would merely introduce an additive constant W,—W, 
independent of angle, inasmuch as our models of these 
molecules always involved the aggregate configuration sp* 
whereas the normal state of C is s*p?. It must be cautioned 
that even with such additions we make no allowance for 
the dependence of the energy of the central atom on L 
and S in distinction from on the individual /’s. The correc- 
tion for the L, S structure will be studied in Part III. 

In the cases of the molecules CH»X»2, CH3X, or CHXs3, 
attempts could be made to calculate quantitatively the 
deviations of the angles from the tetrahedral values by 
finding what values of w,w’ or 6 minimize (27), (37), 
(38), and (39), and so determine the values of the constants 
(3) or (36) which will yield agreement with experiment. 
However, this procedure involves algebraic difficulties not 
encountered in NH;, where simple formulas for @ were 
derived, and is promising only in a rough way, since, as 
in NHs, the precise minimizing angles will depend on the 
type of approximation used. For instance, the minimum 
of (27) with (24-25) does not yield values of w, w’ con- 
forming to (32). If in particular, as in the method of 
Section 4, one neglects all s bonding effects by setting 
P,=Q,=0 in (24-25) the minimum of (27) is easily found 
and comes at 


P,? cos? w=(Q," cos? w’ = 3(Ps?+Q,") if 
on’, dat if 
w’ =45° if 


$=P,?/Qe? — > f 
P7/Q°=2; 
Po? /Q.?S}-. 


w=0, 
This is for CH2X2; for CH;X (Eqs. (30) in (27)) the 
corresponding values are 


90," sin? 0=60,?+2P,* if Po?/Q.2=3; 


6=90° if P?/Qf=}. 
(If instead, (39) is used, the neglect of s bonding means 
B=0, and the minimizing angles for CH;X, for instance, 
are the same as those just given except that yyx, YHH 
replace P,?, Q,”. If the approximations (I), (II), (III) of 
Section 3 are made, yyHu/yvux reduces to Q,?/P,?, as should 
be.) The fact that t—2w=6=90° if P?>>(? is, of course, 
an illustration of the tendency of the most important bond 
angles to become 90° if the p bonding effects predominate. 
At this point we may mention a modification which we 
believe should be made in some of the statements in 
Pauling’s important paper.!! It conveys the impression 
that even when there are only two atoms attached the 
tetrahedral angle 109.5° is the best possible one, whereas 
we have already intimated that this angle is characteristic 





only of four equivalent attached atoms, and that two sole 
(or the two firmest out of many) attached atoms tend to 
set their bonding axes at an angle between 90° and 180° 
whose precise magnitude depends on the relative im- 
portance of pure p-bonding and s-p interplay. This can 
be seen not merely from the near-ionic approximation 
(39) but also with the two other methods of approximation 
for the Hund-Mulliken theory given in Sections 3-4, or 
with the Heitler-London-Pauling-Slater method. The 
assumptions of Section 4 largely suppress s bonding and 
so make the valence angle in CX2 nearly 90°. To calculate 
this angle with the procedure of Section 3 one has only 
to set Q=0 in (24-25) since CX.» is formally equivalent 
to CH2:X2 with the C—H bonds made infinitely small. 
The minimum of (27) then comes at P,? cos 2w=P, if 
P.=P, and at 2w=180° if P,=P,; whether there are 
four, eight or some intermediate number of electrons is 
immaterial in this particular case, since two of the roots 
(24-25) are nonbonding when Q=0, making any electrons 
beyond four superfluous. Similarly in the Heitler-London- 
Pauling-Slater method, the behavior of CX2H2 in the 
limiting case that the C—H bonds are very weak is given 
by Eq. (38) if one sets N=0 and remembers that w’ is a 
function of w given by (36). An explicit formula for the 
value of the angle r—2w between the two C—X axes 
which then minimizes (38) cannot be given because of 
algebraic difficulties, but the terms of (38) which involve 
Noe have their minimum at 2w=90° and those which 
involve Nz; or Mss, their minimum at w=0, making it 
clear that the minimum will in general come at a value 
of +—2w intermediate between 90° and 180° whose precise 
magnitude depends on the relative sizes of Voc, Ves, and 
N.ss Terms containing M,, are unimportant because Ire 
is small. The reason why Pauling appeared to obtain 
109.5° quite generally is that he used semiquantitative 
criteria of ‘overlapping’ rather than the more exact Eq. 
(35) and assumed a particular value of the relative intensity 
of the s and p bonds. Of course 109.5° is in the critical 
range 90°-180°, and may not be far removed from the 
actual best angle since the p bonding effects favoring 90° 
are probably the most important factor. In fact, according 
to Mecke,!? the actual angle in HO is between 103° 
and 106°. 

CX:, BeX2, OX2. We must mention that when the H 
atoms are completely removed from CH2X2, the wave 
functions appropriate to the third and fourth electrons of 
C are probably not correctly given by y,°, y, in (34). 
When the H atoms are missing, these two electrons have 
no attached atoms with which to be paired, and so need 
not have equal hybridization ratios as assumed in (34). 
Also they may even be in the same state, with their spins 
mutually paired. Another possible choice for the third 
and fourth wave functions in place of those in (34) is'8 


17R, Mecke and W. Baumann, Phys. Zeits. 33, 833ff. 
(1932). 

1% It seems quite probable that the functions (42) and 
first two of (34) are the proper wave functions for C 
double bonds, rather than the tetrahedral sp* functions 
assumed by Pauling and Slater. At any rate it is certain 
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Vel _ ¥(2ey), 
va’ =[k(w) J#[ —cosec w(cos 2w)4y(20:)+y(2s)] (42) 


with k(w) =1/(cosec? w cos 2w+1). Here the y axis has been 
taken perpendicular to the plane of X -C—X, and the z 
axis bisects the angle X—C—X. The exchange energy 
between the third and fourth electrons and the X atoms 
is minimized by placing both these electrons in the state 
vc’. The total energy is"® then 


W= W’ —2Nss+2Nee+[1 +k(w)][Wp- Ws], 


where W’ denotes the value of the expression (38) when 
one sets N=O and strikes out all terms involving f(w’) 
(including those through g(w’, w)). More likely, however, 
the third and fourth electrons both decide to settle in ya’, 
as this makes the carbon configuration nearly s*p? rather 
than nearly #* if 2w is nearly 90°. The total bonding energy 
is then 


W=W’'—2%..+[1—k(w) J[W,— W. ] 
+2k(w)[cos 2w(Nos+cot? wrx) 
—2(cos 2w) Nos +Nee ]. (43) 


If the third and fourth electrons are entirely missing, as 
in BeXe, this energy is W’—2%,s+[1+k(w) ][W,— W-]. 
If there are six electrons in the central atom, as in H,O 
or OXse, there are probably two electrons each to the 
states y.’, va’, and the energy is the expression (43) 
augmented by 29,7. In all cases the first two electrons 
are assigned one each to the first two states of (34). 
None of these modifications change the critical region in 
which the best value of the angle X —C —X is located. 
The foregoing discussion of CX2, BeX2 or OX: pre- 
supposes the two attached atoms to be identical. Calcu- 
lations can also be made for the case that they are not 
equal, but we shall not give them, as they yield nothing 
qualitatively new regarding the favored angles. When 
the two attached atoms are not identical and hence not 
held by bonds of equal strength, one can no longer assume 
that they are held by pairing with electrons of the central 


that the s-p hybridization ratios for double bonds will not 
be the same as for CHy. We suggest this combination of 
(42) with (34) because it gives three directed coplanar 
valences, viz., those corresponding to vee, We’, wa’. The 
C—C bond is here to be located along the z axis, so as to 
bisect X—C—X. This subject of double bonds will be 
discussed more fully in a paper by W. G. Penney. 

1” This and shortly following formulas for the bonding 
energy are obtained by application of (35) and (6), with 
allowance for the finite s-p separation in the same fashion 
as in connection with (41). The addition of —29., and 
the omission of f(w’) are consequences of the fact that in 
(38) those terms containing Nos, Nos, Mex Which arise from 
interaction of the third and fourth electrons of C with X:, 
are proportional to f(w’), while the corresponding terms 
containing Mt,, are proportional to 1—f(w’). Such terms 
must be dropped and new terms added because of the 
substitution of new wave functions in place of the third 
and fourth of (34). 
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atom which are represented by the first and second wave 
functions of Eq. (34). Instead the s-p hybridization 
ratios appropriate to the two bonds may be unequal, 
whereas (34) assumes equality. 


7. STRUCTURE oF CH; AND CH,? 


All of our conclusions concerning the superi- 
ority of the tetrahedral or distorted tetrahedral 
arrangement are contingent upon there being 
eight electrons available for bonding, so that 
each of the four lowest roots of the secular 
Eq. (4) are filled twice. When there are only 
seven electrons available, as in CH; and CH,*, 
the results are sometimes different. In the 
Hund-Mulliken method they depend upon the 
particular type of approximation which can be 
employed, and are as follows. 

(i) The near-ionic approximation of Section 2, 
if allowable, predicts that the most stable models 
of both CH; and CH," are plane provided the s 
state of C is beneath the p states by an energy 
difference large compared to the bonding energy. 
This can be seen from the following simple 
argument. In the near-ionic approximation, by 
far the largest part of the linear combination 
involved in the wave function is the carbon 
portion, so that it has a meaning to label the 
electrons 2s, 2p, etc., as in carbon. We showed 
in Section 2 that the energy of an eight electron 
system is invariant of the spatial arrangement 
to the approximation )?. The bonding power of 
any individual electron, however, is of the order 
\*. Hence the arrangements of the H atoms in 
CH; or CH,*+ which give the lowest energy will 
be the same as the arrangements of CH;~ or 
CH, for which one 2 electron is removed most 
easily. The plane models possess this character- 
istic, since when all the H atoms are in the xy 
plane, the wave function ¥(2c.) is nonbonding, 
except for the small term C, in (3). Another 
way of saying the same thing is that with the 
seven electron configuration 2s?2¢,°20,?2¢., in- 
creased overlapping is secured by placing all the 
H atoms in the xy plane, whereas this is clearly 
not true of the more symmetrical, eight-electron 
configuration 2s?2¢,220,220,7. 

(ii) The method of approximation of Section 3 
is the same for CH; as NH. The determination 
of the best value of @ is unaffected by the fact 
that the energy is now 2W,+2W2+2W,+W. 


instead of the sum (27) appropriate to an eight 
electron system, for W_ in (30) is a nonbonding 
root independent of 6,Q when P=0. Hence @ 
is given by (40) and the model is pyramidal if 
the s bonds are not too intense. On the other 
hand one finds with this method” that CH,+ 
will be plane even if the s bonds are weak unless, 
also, Q, is sufficiently large compared to W” — W°, 

(iii) With the method of approximation of 
Section 4, which is near ionic only for the 2s 
electrons, one concludes that the best models of 
CH; and CH,* are respectively models of type 
(c) with atom a missing and type (a) shown in 
Figs. 1, 2, Section 3, with the following values of 
6 or w, w’ ; 


8 1 
sin? @=-+——— (W"” — W,,°)? 
9 180,2 


if (W" —W,°)?<2Q,?; 





cos? w=cos?w’ =—+ (W¥—W,°)? 
9 240, 
if (W* —W,°)?<4Q,?; 
§6=7 or w=w’=0 otherwise. 


These results are obtained by minimizing”! 
2Wi+2W2+ W. with P,=P,=Q,=0 in (30) or 
with P,=Q,, P,=Q,=0 in (24-25). The values 
6=7, w=w’=0 represent plane models, and so 
we conclude that the spatial model is obtained 
by the procedure of Section 4 only if Q is suffi- 
ciently large compared with W” — W,°. 

To study CH; by the H-L-P-S method one 
uses the wave functions (33) and sets J{=0 in 
(37). The minimum of (37) then falls at” 


(Noe _ Nas) 


[4 Nee + (Neo — N,s)* }* 
if Neo= Nes t2-'Nes (44) 


20 For instance when P=Q, w= +w’, one finds W,+2W 
<W,'+2W%' only if Q,? cos? w>3(W" — W°)?. Here Wi, W2 
and W,’, W2' are levels for certain spatial and plane 
models as given in (24), (28). The remaining roots (25) 
are common to both models. 

21 Only five electrons are involved in this sum because 
two are absorbed by the deep 2s states. 

22 Eq. (44) would also give the best value of @ for NH; 
in the H-L-P-S method were it not for the fact that the 
energy of NH; involves the additional terms (41). The 
term proportional to W,—W, in (41) will tend to make @ 
smaller in NH; than in CHs. 





cot? @= —j+ 


wloo 
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and at @=7 otherwise. Hence this method 
predicts a pyramidal model only if N,, be 
sufficiently large compared with N,,, Nos. The 
structure of CH,+ cannot be simply studied by 
the H-L-P-S procedure since no natural scheme 
of electron pairing is possible. 

In view of all the above it is likely that CH; 
has a larger value of @ than does NH3;, and may 
conceivably even be plane. One thing is particu- 
larly clear: there is no reason for CH,* to be a 
regular tetrahedron like CH, Instead CH,* 
should be a flattened tetrahedron, i.e., a model 
of type (a), Section 3, with w= —w’, | w|=35.3°, 
and possibly w= w’ =0. 


H. VAN VLECK 


With all three variants (i-ii-iii) of the H-M 
method, one electron of CH; is more loosely 
bound than the rest. In fact, in (ii) one root is 
nonbonding, and in (i) nearly so. The H-L-P-S 
method also predicts one loose electron in CHs, 
viz., that electron of C which is unpaired to an 
H partner. Experimental confirmation is pro- 
vided by the fact that the ionization potential 
10.3 volts of CH; is lower than that, 14.5, of 
CHa, an observation due to Mulliken.” From 
qualitative consideration of the wave functions, 
he has, in fact, suggested a pyramidal model of 
CHs with one loosely held electron. 


23R. S. Mulliken, Phys. Rev. 40, 56 (1932). 
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The eigenfunctions corresponding to the fourteen singlet 
states of a system of seven or eight electrons are given 
and the corresponding matrix elements tabulated in a 
form convenient for computation, neglecting multiple 
interchange integrals. Proof is given that the general 


method described in earlier work for writing bond eigen- 
functions suffices to separate the multiplets. A method 
for calculating matrix components between bond eigen- 
functions directly from the bond systems is also described. 





N a series of papers from this laboratory! the 

problem of five or six spherically symmetric 
atoms, each with one valence electron, has been 
solved to the first order of the perturbation 
theory of quantum mechanics and the results 
applied to the calculation of the activation ener- 
gies of some typical chemical reactions. In the 
third paper of this series it was found to be 
necessary to extend the results to the case of 
eight atoms. In this paper we shall give the de- 
tails of these calculations and at. the same time 
supply proofs of certain theorems which were 
assumed in the earlier papers of this series. 

The method of approach is an extension of a 
method used by Slater? for the problems of three 
and four atoms. Since the seven-electron problem 
is immediately obtained from the eight-electron 
problem by putting one electron at infinity, it 
need not be considered separately. Let the eigen- 
functions (neglecting spin) of the eight atoms be 
a(x, y, 2), b(x, y, 2), «++, h(x, y, 2). Let the co- 
ordinates of the eight electrons be x, 1, 21, %2, 
‘++, 23. We shall abbreviate such an expression as 
a(x, 1, 21) to a(1). Then an approximate solution 
of the problem of eight atoms is the eigenfunc- 
tion a(1)b(2)c(3)---h(8). Any function which 
may be obtained from this by permuting the 
numbers 1 to 8 is an equally good solution. 

If we now introduce the spin of the electron, 
we find that instead of the 8! solutions we have 
just obtained, we have 2°8! solutions since with 


‘Sherman and Eyring, J. Am. Chem. Soc. 54, 2661 
(1932); Kimball and Eyring, J. Am. Chem. Soc. 54, 3876 
(1932); Taylor, Eyring and Sherman, J. Chem. Phys. 1, 
68 (1933), 

*Slater, Phys. Rev. 38, 1109 (1931). 


each space eigenfunction we may either associate 
the eigenfunction a representing the electron spin 
of +1/2 h/2r or the eigenfunction 8 representing 
a spin of —1/2 h/2z. The Pauli exclusion prin- 
ciple, however, limits us to the 2° linear combina- 
tions of these eigenfunctions which are antisym- 
metric in the electrons. Following Slater, we shall 
take as these antisymmetric combinations the 
determinant 


aa(1) ha(1) 
1 |aa(2) ba(2) ha(2) 


as Ves , 
(8!)4] ; ; ; : (1) 


ba(1) ++ 








aa(8) ba(8) ha(8) 


and the 2°—1 determinants which may be ob- 
tained from this by replacing any column or 
columns of a’s by f’s. We shall denote such an 
eigenfunction by a symbol such as 


pe Se Bae 


aaPBpaaBpeaea 


this denoting that determinant in which the c 
and f columns contain #’s and the rest a’s. 

Each of these eigenfunctions is an exact eigen- 
function of S,, the operator giving the total 
z-component of spin, since each term of the ex- 
pansion of the determinant is an eigenfunction of 
S, with the same eigenvalue. We may, therefore, 
classify these eigenfunctions according to their 
values of S., which are easily seen to be (%.— 1g) 
<h/4mr where n, is the number of columns of a’s 
and mg is the number of columns of §’s. The 
number of eigenfunctions with each value of S, 
is given in Table I. By a well-known theorem, if 
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spin interactions are neglected the matrix ele- 
ments of the Hamiltonian operator H between 
eigenfunctions with different values of S, all 
vanish, so that in the perturbation theory we 
can consider each group separately, instead of 
using all 256 eigenfunctions at once. 


TABLE I, 


S,(in unitsh/2x7) 43 2 1 0-1 -—2 -3 -4 
No. eigenfunctions 1 8 28 56 70 56 28 8 1 


Although these eigenfunctions are eigenfunc- 
tions of S,, they are not in general eigenfunctions 
of S*, the square of the total angular momentum. 
It is possible, however, to combine these eigen- 
functions in such a way that the new eigenfunc- 
tions are eigenfunctions of S? as well. Consider, 
for example, the 70 eigenfunctions for which 
S,=0. We shall show that from this set we can 
form 70 linear combinations which are, of course, 
still eigenfunctions of S, and in addition are now 
eigenfunctions of S*. Of these new eigenfunctions, 
1 will have the eigenvalue 20(4/27)?, 7 the eigen- 
value 12(4/27)*, 20 the eigenvalue 6(4/27)*, 28 
the eigenvalue 2(h/27)? and 14 the eigenvalue 0. 
Moreover, these eigenfunctions will have no 
matrix elements of H between the various groups. 

In order to find these combinations we define 
for each pair of atoms ab, ac, etc., a function 
5(k). This function 6(k) has the value 1 if in 
the k** eigenfunction in the set of 70 the spin of 
ais a and the spin of b is 8; it has the value —1 if 
the spin of a is B and the spin of } is a; and it has 
the value 0 if @ and b have the same spin. By 
means of these functions we may construct an 
eigenfunction corresponding to any way of draw- 
ing bonds between the eight atoms. Suppose, for 
example, that we wish the eigenfunction corre- 
sponding to the bonds a-b c-d e-f g-h. The re- 
quired eigenfunction is then 


Dnbav(k) Seah) 5e¢(k) 5on(k) Wes 
where y; is the k** eigenfunction in the original 
set of 70. Although there are 105 ways of drawing 


fo ee 


B Baauwciaaeaea 





(S.+48,)( 


§ Pauli, Zeits. f. Physik 43, 601 (1927). 
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4 bonds between 8 atoms, only 14 of the resulting 
eigenfunctions are linearly independent. 
We shall now show that these bond eigenfunc- 
tior are eigenfunctions of S*. By definition 
S= $7+S,7+S. (2) 
Remembering the commutation rules 
S,Sy—SySz:=th/2x S; 
SyS2—S,Sy=th/2r Sz 
S,S2:—S,S,=th/2r Sy 
we may write (2) in the form 
S?= (S,+4S,)(S2—-tSy) —h/2aS2+S2 (3a) 
S?= (S,—4Sy)(S2+iSy) +h/27S.4+S2. (3b) 


It follows that if a function is an eigenfunction of 
S, and S? it must also be an eigenfunction of 
(S2+i5,)(S;—-aiS,) and of (S,—i5S,)(S,+7S,). 
Now the operator 


S2t1Sy= Sri ttSyit S2ottSyo+ th oy + SzstiSys 


where S,; and S,; are the operators for the x and y 
components of the spin of electron 1. Since® 
Sra(1)=(h/4r)B(1); = S218(1) = (h/477) (1) 
Syra(1) = (th/47)B(1); — SysB(1) = (—th/4 7) (1) 
we have (Sz1+7S,:)a(1)=0, 

(Szi+7¢Sy1)B(1) = (A/27)a(1), 
and also (Sz1—7Sy1)a(1)= (h/27)B(1), 

(S21 —4Sy1)B(1) = 0. 


Let us now consider the effect of the operator 
S,+i5S, on a typical Slater eigenfunction. Acting 


= eee rt ss 


a @2 a@ a@ @ @ A@ 04 


the operator S,+iS, obviously gives 0. Acting on 
any other eigenfunction of this type it is easily 
seen that the result is 4/27 times the sum of all 
the determinants which may be obtained by re- 
placing one column of 6’s by a’s. Thus 


el pee ee, 


Baaeaaeaae 


2r 











































of 


or 
ng 


on 
ily 
all 











Similarly the operator S,—72S, gives 0 when it 
acts on 

(| bcdef g 3) 

B BBBRB B B B 


and acting on any other Slater eigenfunction 
gives h/2m times the sum of all the determinants 
which may be obtained by replacing one column 
of a’s by B’s. 

It is now easily verified that the operator 
S,+iS, acting on any bond eigenfunction cor- 
responding to four bonds gives 0 as a result. 
Hence, the operator (.S,;—75S,)(.S,;+7S,) also gives 
0, and since S, operating on any of the set of 70 
original eigenfunctions gives 0, it follows that the 
operator S®? acting on this eigenfunction also 
gives 0. The bond eigenfunctions corresponding 
to 4 bonds are, therefore, eigenfunctions of S? 
with the eigenvalue 0. A similar argument shows 
that those bond eigenfunctions corresponding to 
3 bonds are eigenfunctions of S* with the eigen- 
value 2(h/27)*, those corresponding to 2 bonds 
have the eigenvalue 6(//27)?, those corresponding 
to 1 bond have the eigenvalue 12(h/27)?, and 
that the single eigenfunction corresponding to no 
bonds has the eigenvalue 20(h/27)*. It follows 
that the matrix components of H between bond 





Hi,1—-d1,1E 
ITe,1—d2,1E 


IT1,2-—d1,2E 
He,2—d2,2E 





where 
Hy= { vetlivade, d= [ vtvde. 


We must, therefore, evaluate these integrals. 

We shall suppose that the field due to all the 
electrons in each atom except the valence elec- 
tron may be represented by a function V(x, y, 2), 
distinguishing the atoms by subscripts as 
V.(x, y, ). Then the Hamiltonian function takes 
the form 


H= D8 a1 —h?/822mvy2 
of b, eee), tem Vo Xi, Vis i) + DL i>i(e?/rij). 
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Hy4,1—-di4,1£ TH3,2—dy4,2E 
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eigenfunctions corresponding to different num- 
bers of bonds all vanish. 

Since in chemical problems it is almost always 
the state of lowest multiplicity which is of inter- 
est, we may from now on confine our attention to 
a set of fourteen linearly independent bond eigen- 
functions with S,=0 and S?=0. The actual set 
of fourteen used were those corresponding to the 
ways of drawing the bonds shown in Table II. 








TABLE II. 
yi a-b c-d e-f g-h vs a-c b-d e-g f-h 
v2 a-b c-e d-f g-h yo a-c b-e d-g f-h 
vs a-c b-d e-f g-h Yio a-d b-e c-g f-h 
vs a-c b-e d-f g-h - Wu a-h b-g c-d e-f 
vs a-d b-e c-f g-h Yi. a-h b-g c-e d-f 
ve a-b c-d e-g f-h Yis a-h b-d c-g e-f 
¥z a-b c-e d-g f-h via a-b c-g d-h e-f 








These eigenfunctions are independent, and since 
all other four-bond eigenfunctions can be ex- 
pressed in terms of them they form a complete 
set. 

By using these fourteen eigenfunctions as a 
basis for the usual perturbation theory, the 
energy is given to the first order by the secular 
equation 


IT, u—d, ue 


He, 14—d2, 14 E 
=0, (4) 





i, u—di4, ui 





Since the eigenfunction a(z) is by hypothesis a 
solution of the equation 


(—h?/82?my 2+ Va(xi, yi, 2:))a(t) = E,a(t), 


where E, is the lowest energy level of atom a, we 
may replace the operator —(h?/8x’m)V?2 by 
E.— V.(i) whenever it operates on a(i) and by 
similar expressions when it operates on the other 
atomic eigenfunctions. The matrix components 
of H between the Slater eigenfunctions have been 
discussed by Slater.? If we neglect the lack of 
orthogonality of the eigenfunctions a(i), 6(7), 
c(i): ++, the multiple interchange integrals are all 
small and can be neglected. To this degree of ap- 
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proximation the matrix components are then 
given by the following rules: 

(1) The matrix component between any Slater 
eigenfunction and itself is the coulombic integral 


C= face): - «h(8)I1a(1)b(2)+ + +h(8)dr 


minus the sum of all the interchange integrals 
between atoms having the same spin. The.-inter- 
change integral between two atoms, say a and 3, 
is the integral 
(ab)= | a(1)b(2)-++h(8)Z7b(1)a(2)- + + h(8)dr. 
(2) The matrix element between two different 
Slater eigenfunctions is zero unless the two 
eigenfunctions differ only by a single interchange 
of spins on two atoms, in which case the matrix 
component is the negative of the interchange in- 


tegral between these two atoms. For example, 
the diagonal matrix element for the eigenfunction 


Ce. eS 
aaaaBp Bp BB B 


is C— (ab) — (ac) — (ad) — (bc) — (bd) — (cd) — (ef) 
— (eg) — (eh) — (fg) — (fh) — (gh), while the matrix 
element between this and 
fg ;) 
B B B 
is — (ae). 


abede 
“ aaaea 

These rules enable us to find the matrix ele- 
ments between the various bond eigenfunctions. 
Suppose, for example, that we wished to find the 
matrix element between the eigenfunction repre- 
senting the bond system a-b c-d e-f g-h and 
that representing the bonds a-c b-d e-g f-h. 
We could, of course, actually write out these 
eigenfunctions in terms of the Slater eigenfunc- 
tions and, knowing the matrix components of 1 
between these, find the required matrix element. 
It is, however, possible to set up rules for finding 
the matrix element directly. It is evident that the 
matrix element we seek will be the sum of some 
multiple of C and certain multiples of the various 
interchange integrals. We shall now show how the 
coefficients of the various integrals in the matrix 
element may be simply determined. Let us first 
consider C. The only way C can occur in the re- 
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sult is from the presence of the same Slater 
eigenfunction as a part of both bond eigenfunc- 
tions. The coefficient of C is, therefore, equal to 
the number of Slater eigenfunctions which occur 
once in each of the two bond eigenfunctions, 
taken with a+ or — sign according to whether 
they occur with the same or opposite signs in the 
two bond eigenfunctions. In order to find this 
number we write the two bond systems above one 
another, thus: 


a-b cd ef gh 
a-c b-d e-g f-h. 


We now start assigning spins to the various 
atoms. Taking a arbitrarily to have the spin a, 
we must have the spin 8 on bd in order to have the 
bond a—b. We, therefore, write a on the a’s and 
6 on the 0’s. In the same way we see that the spin 
of c must be 6 and that of d must be a. We now 
have the following picture: 


aBp Ba 
a-b c-d e-f g-h 
a-c bd e-g f-h 
ap Ba. 


We may now choose the spin of e arbitrarily, but 
then all the spins are determined, thus: 


aB Ba aB Ba 
a—b c-d ef gh 
b-d e-g f-h 
Ba aB Ba. 


a-c 


a B 


Since we could have chosen a to be either a or 8, 
and likewise e to be either a or , there are four 
Slater eigenfunctions common to the two bond 
eigenfunctions. It is easily seen that these eigen- 
functions all enter both bond eigenfunctions with 
the same sign, so that the coefficient of C in the 
matrix element is +4. 

In order to find the coefficient of an interchange 
integral such as (ab) we proceed in a similar way. 
We note first that such an integral can arise in 
two ways: it may arise from the presence of a 
Slater eigenfunction common to both eigenfunc- 
tions in which a and b have the same spin; or it 
may arise from the presence of two Slater eigen- 
functions differing only in the interchange of the 
spins of a and b. The two cases are considered 
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separately. To obtain the contribution of the first 
case we again put down our two bond systems 
and now write a’s on the a’s and also on the b’s. 
We then proceed to fix the other spins by the re- 
quirement that wherever there is a bond between 
two atoms the spins must be opposite. The 
process is continued until every spin is fixed or 
an inconsistency is reached. For example, taking 
the same bond system as before the result for 
(ab) is immediately inconsistent because of the 
bond a-d in the first eigenfunction. If we take the 
integral (ae), however, we find 


aB Ba aB Ba 
a-b c-d e-f g-h 


a-c bd e-g f-h 
aB Ba aBp Ba. 


There are, therefore, only two Slater eigenfunc- 
tions giving rise to the integral (ae) (one with 
a’s, the other with 6’s on a and e). Each has a 
positive sign, so that the contribution to the co- 
efficient of (ae) is —2 (the negative sign arising 
from the original rules for matrix elements be- 
tween Slater eigenfunctions). Taking the integral 
(ad) as a further example we find 


aB Ba aB Ba 
a-b c-d ef g-h 


a-c bd e-g f-h 
aBp Ba aB Ba. 


In this case, when the spins of a, b, c and d were 
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assigned we were left free to choose either spin 
for e, so that there are four eigenfunctions giving 
the integral (ad). The contribution to the coeffi- 
cient of (ad) is, therefore, —4. 

In order to find the other part of the coefficient 
of (ab) we write a on a and £ on b in the first bond 
system and 6 on a and a on B in the second, con- 
tinuing the process as before. Thus we find 


aB aB aB Ba 
a-b c-d e-f g-h 


a-c bd e-g f-h 
Ba aB aB Ba. 


The sign of the first is negative and the second 
positive, and we could have chosen a@ or 8 on a 
and also on e. The contribution is, therefore, +4. 
The total coefficient of (ab) is, therefore, +4. 
The diagrams for (ae) and (ad) and their contri- 
butions are 


ap Ba B 
a—b c-d ef g-h 
(ae) (inconsistent) 
a-c bd e-g f-h 
BB aa 


aB aB aB Ba 

a—b cd ef gh 
(ad) (—4) 

a-c bd e-g f-h 

Ba Ba aB Ba. 
The total coefficient of (ae) is, therefore, —2 and 
that of (ad), —8. 


TABLE III. 
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Having found the matrix of Hin this way it isa 
simple matter to find the matrix of unity, for if we 
neglect the integrals arising from the lack of 
orthogonality of the eigenfunctions, d;; is just 
the coefficient of C in H;;. If the various inter- 
change integrals and the coulombic integral can 
be found, we can, therefore, set up the secular 
equation and complete the first order perturba- 
tion theory. 

Since the bond eigenfunctions are not orthog- 
onal, the secular equation contains E in every 
element. It is more convenient for numerical work 
to orthogonalize the fourteen bond eigenfunc- 
tions. The set of orthogonalized eigenfunctions is 
shown in Table III. The first eigenfunction in this 
table is just the first of the set of bond eigen- 
functions of Table II. Each succeeding one is ob- 
tained by orthogonalizing the corresponding 
bond eigenfunction to the eigenfunctions occur- 
ring earlier in the table. 

It is not as simple to compute the matrix com- 
ponents between these eigenfunctions as it was 
for the bond eigenfunctions themselves. It was 
found best to express these eigenfunctions in 
terms of the Slater eigenfunctions and compute 
the matrix elements from the rules for forming 
matrix components between the Slater eigenfunc- 
tions. The results are shown in Table IV. In the 
calculation of these matrix elements each eigen- 
function was multiplied by the arbitrary factor 
1/(2)}. 

The greatest advantage gained by the orthog- 
onalization comes from the fact that in sym- 
metrical configurations and also some configura- 
tions in which some of the interchange integrals 
can be neglected a proper assignment of letters to 
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the atoms will give a secular equation which may 
immediately be factored into smaller deter- 
minants. A brief inspection of the various possi- 
bilities is usually sufficient to determine the cor- 
rect assignment of the letters in these cases. It is, 
of course, always convenient to make as many 
matrix components zero as possible and im- 
portant that the maximum number of these fall 
in a particular row or column. Having done this, 
three evaluations of the determinant for proper 
choices of E ordinarily suffice to estimate the 
lowest root of the equation. 

Not every choice of fourteen independent bond 
eigenfunctions would lead to orthogonal eigen- 
functions which have as simple matrix com- 
ponents as the set chosen here. To secure sim- 
plicity each bond eigenfunction should differ 
from the others as little as possible without losing 
independence. 

The seven electron problem is, of course, simply 
a special case of the eight-electron problem in 
which one of the electrons is at infinity. This 
amounts simply to setting equal to zero all the 
integrals containing one particular letter, say h. 

The problem still remains: to what systems do 
these simplified calculations apply? At present 
this question seems best answered by a detailed 
application to a variety of cases. One of these 
cases we have already referred to and others are 
being investigated. In any case it seems logical to 
consider first the approximations given here, 
since the results may readily be extended to in- 
clude the neglected multiple interchange inte- 
grals and orthogonality integrals should this be 
found desirable. 
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The Problem of Atomic Polarization 


CHARLES P, Smytu, Frick Chemical Laboratory, Princeton University 
(Received January 28, 1933) 


New values of atomic polarizations are obtained from 
the temperature variation of the dielectric constants of 
gases. Although the values are of a greater magnitude 
than those calculated for simple molecules from infrared 
intensity measurements, they show that the improbably 


high values previously found for several alkyl halides are 
incorrect. If the molecule is not very large, its atomic 
polarization is a relatively small quantity unless dipole 
rotation may occur within the molecule as indicated by 
the dielectric constant of the substance in the solid state. 





INTRODUCTION 


N an earlier paper! the writer collected the 

available data on atomic polarization and at- 
tempted to investigate such relations as might 
be found among them, a summary of this treat- 
ment being published subsequently. At that 
time Van Vleck* had calculated from infrared 
intensity measurements that the vibrational or 
atomic polarization was negligibly small for HCI, 
HBr, CO, CO., NH; and CH, and attributed to 
experimental error the small but not incon- 
siderable values found for them, indicating, how- 
ever, that this did not necessarily mean that 
other kinds of molecules might not possess sig- 
nificant atomic polarizations. In a recent treat- 
ment of these simple molecules,‘ he expresses 
scepticism as to ‘‘whether any information about 
the order of magnitude of the atomic polarization 
can be deduced from existing absorption data”’ 
and submits as “‘perhaps the best appraisal . . . 
that the effect of the atomic polarization on the 
dielectric constant is negligible in stable diatomic 
molecules, but not necessarily in molecules with 
more than two atoms.” Smallwood? has recently 
attempted to calculate atomic polarizations but 
has obtained values much smaller than the ob- 
served. Thus far we have failed to account in a 





1C. P. Smyth, J. Am. Chem. Soc. 51, 2051 (1929). 

*C. P. Smyth, Dielectric Constant and Molecular Struc- 
ture, p. 163, The Chemical Catalog Company, Inc., New 
York, 1931. 

*J. H. Van Vleck, Phys. Rev. 30, 31 (1927). 

‘J. H. Van Vleck, Electric and Magnetic Susceptibilities, 
P. 51, Oxford University Press, Oxford, 1932. 


*H. M. Smallwood, Zeits. f. physik. Chemie [B] 19, 
242 (1932). 
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general quantitative fashion for the atomic 
polarization. The present paper offers new and, 
for the most part, more accurate data on atomic 
polarization, which will be used to answer certain 
newly arisen questions and to eliminate certain 
incomprehensibly large values from among the 
earlier data. 


CALCULATION OF RESULTS 


It will be remembered that the atomic polariza- 
tion P4 is obtained by subtracting the electronic 
polarization Pr and the orientation or dipole 
polarization Py from the total polarization 
P=(e—1)M/(e+2)d, in which « is the dielectric 
constant, M, the molecular weight, and d, the 
density, that is, P4 =P—Pr—Py. If Py can be 
eliminated by fixing the dipoles in the solid state, 
Pa may be obtained by subtracting from P the 
value of Pz calculated by extrapolating the 
molar refraction to infinite wave-length. Some of 
the difficulties of this method have already been 
discussed and others which have been made 
evident by recent work on solids will be treated 
presently. Since, in the Debye equation, P=a 
+b/T, a=P2e+P., Pa can, perhaps best, be 
determined from a, which is obtained from the 
temperature variation of P. When the previous 
tabulation of values of P4 was made, few accurate 
measurements of the temperature variation of the 
dielectric constants of gases were available and 
the extreme inaccuracy of some was not fully 
recognized. Moreover, as the present work will 
show, the apparent temperature variation of the 
polarization in solution is sometimes misleading. 
As previously pointed out, the values of P. are, 
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at best, inaccurate, since they are determined as 
small differences between relatively large quan- 
tities, all the errors of experiment and method 
being thus accumulated in them. The values in 
Table I, assembled by Dr. Kenneth B. McAlpine, 
have been calculated from the temperature 
variation of the dielectric constants of the sub- 


TABLE I. 








Atomic polarizations 





Nitric oxide 0.4¢ Methyl fluoride 2.4¢ 
Nitrous oxide 0.4¢ Methyl chloride 2.2¢ 
Carbon disulfide 0.3/ Methyl bromide 1,2¢ 
Propylene 0.64 Methyl iodide 1.4¢ 
Propane 0.34 
Benzene 1.94 Ethyl chloride 1.4¢ 
Toluene 2.54 Ethyl bromide 2.9¢ 
n-Heptane i1.1¢ Ethyl iodide 1.5¢ 
Ethylene oxide 1.89 CHFCl, 0.9% 
Ethyl ether 3.99 CHF,Cl 3.9% 
3.34 CF.Cl, 4.1° 
CFCl; 3.0° 








*C, P. Smyth and K. B. McAlpine, J. Chem. Phys. 1, 
60 (1933). 

’C. P. Smyth and K. B. McAlpine, ibid. 1, 190 (1933). 

°C, P. Smyth and K. B. McAlpine, unpublished meas- 
urements. 

4K. B. McAlpine and C. P. Smyth, J. Am. Chem. Soc. 
55, 453 (1933). 

€C. H. Schwingel and J. W. Williams, Phys. Rev. 35, 
855 (1930). 

/C. T. Zahn, ibid. 35, 848 (1930). 

9H. A. Stuart, Zeits. f. Physik 51, 490 (1928). 

* OQ, Fuchs, ibid. 63, 824 (1930). 


stances in the vapor state, which seems to pro- 
vide the most accurate method of determination. 
Some recent data are omitted as being less accu- 
rate. An accurate value for ethyl either is taken 
from the old tabulation for comparison with a 
newly determined value. 


DISCUSSION OF RESULTS 


The atomic polarizations for the small mole- 
cules in Table I are small but not the negligibly 
small values that would probably be calculated 
from infrared intensity measurements. Although 
they are hardly larger than the possible experi- 
mental error, it is important to note that no nega- 
tive values have been obtained, nor have negative 
values been obtained from any dependable 
measurements. This shows that the magnitude of 
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the atomic polarizations cannot be due to experi- 
mental error. In accordance with a previously 
noted tendency, propylene and toluene, which 
have small dipole moments, have larger atomic 
polarizations than the closely related molecules 
of propane and benzene, which have no moments. 
The values for benzene and heptane vapors differ 
by less than the experimental error from those 
previously found for the liquids. 

In measurements upon solutions of poly- 
methylene bromides,® the moments of which in- 
crease with rising temperature because of the po- 
tential energy between the two principal dipoles 
in each molecule, it was found that the values of 
P. calculated from the apparent values obtained 
for the constant a of the Debye equation were 
very high, 60 in the case of pentamethylene 
bromide. As an increase of moment with rising 
temperature would reduce the slope of the 
P—1/T curve and thus give too large a value of a 
and of P,, the apparent high P,4 values were re- 
garded as practically meaningless. This suggested 
the possibility that the surprisingly high values 
of the atomic polarization found for a number 
of alkyl halides!’ from measurements on tem- 


perature variation of the polarization might be . 


due to an increase in moment with rising tem- 
perature. Although not highly probable, it ap- 
peared possible that there might be sufficient 
stretching of the carbon-halogen bond with rising 
temperature to cause an appreciable lengthening 
of the dipole with consequent increase in moment. 
In order to investigate this possibility, the values 
for the methyl and ethyl halides in Table I were 
acquired. The probable error in these values ap- 
pears to be somewhat less than 1. As there is no 
regular trend in the values, it appears probable 
that the true values lie between 1.5 and 2.0. 
The values for the four substituted methanes are 
about twice as great, with the exception of that 
for CHFCl:, in which the error would seem to be 
abnormally large. Although these molecules have 
moments smaller than those of the methyl and 


ethyl halides, they contain three or four principal . 


6C, P. Smyth and S. E. Kamerling, J. Am. Chem. Soc. 
53, 2988 (1931). 

7C, P. Smyth and H. E. Rogers, J. Am. Chem. Soc. 52, 
2227 (1930). 

S. O. Morgan and H. H. Lowry, J. Phys. Chem. 34, 
2385 (1930). 
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dipoles instead of one. Their larger atomic polari- 
zations are in accord with the previous observa- 
tion! that the atomic polarization is greater, the 
greater the number of dipoles in the molecule, 
probably because of the bending or displacement 
of these dipoles within the molecule. 

The value 2.2 for methyl chloride compares 
with the value 8.1 found by Morgan and Lowry 
in hexane solution, 1.3 in carbon tetrachloride 
solution and 8.1 in the solid state, the latter value 
falling to 2.0 as the temperature was lowered 
from — 100° to —190°. The value 1.2 for methyl 
bromide compares with 21.4 in hexane solution 
and 2.8 for the solid at — 100°, decreasing to 0.8 
at — 120°. For methyl iodide, Table I gives 1.4 as 
compared to values of 15.5 in hexane solution 
and 0 in the solid state. For ethyl bromide vapor, 
the value is 2.9 and for the iodide 1.5 as compared 
to 11 found by Smyth and Morgan® for the 
former in hexane solution and 12.3 by Smyth and 
Stoops’ for the latter in heptane solution. The 
low values of the atomic polarizations of the 
vapors, which are dependable and of the same 
magnitude as those found for the solids, show 
that there is no detectable change of moment in 
these molecules in the temperature region 
studied. The high values of the atomic polariza- 
tions found in the measurements upon solutions 
must be regarded as incorrect. 

In the measurements upon solutions each 
value of P used to determine the constants a and 
b for a substance was obtained by extrapolating 
to zero concentration a curve in which the 
polarization of the substance P, was plotted 
against its mole fraction ¢. in a nonpolar solvent. 
As the solvent reduces the forces between the 
polar molecules which lower the observed polari- 
zation, the value found for P2 increases with in- 
creasing dilution and the P,—c. curve rises 
rapidly as zero:concentration is approached. As 
rising temperature diminishes the intermolecular 
action, the curvature of the curves becomes less 
Pronounced so that the extrapolation to zero con- 
centration becomes more accurate, the higher the 
temperature. Since extremely dilute solutions 
were not used in these measurements on the alkyl 


*C. P. Smyth and S. O. Morgan, J. Am. Chem. Soc. 
50, 1547 (1928). 

*C. P. Smyth and W. N. Stoops, J. Am. Chem. Soc. 
51, 3312 (1929). 


halides and since the same general method of 
extrapolation was employed on every curve, it is 
probable that the extrapolated polarization 
values at the lower temperatures are too low and 
approach the correct values only at the higher 
temperatures. Such polarization values would 
give a P—1/T curve of too small a slope and, 
consequently, too large a value of P.4 just as in 
the case of a substance, the moment of which in- 
creased with rising temperature. 

A very pronounced case of the effect of solvent 
upon intermolecular action is shown by measure- 
ments upon butyl alcohol in heptane, cyclohexane 
and benzene solution.’ In benzene solution, the 
moments calculated at 10° intervals from 20° to 
70° deviate little from a mean value 1.74 10-8, 
which is but slightly higher than the value 
1.66 10-8 determined for the vapor.'® In the 
other two solvents, the polarizations in the dilute 
solutions increase with increasing temperature 
instead of falling off as required by the Debye 
equation. Smallwood! calculates from these re- 
sults in heptane solution an apparent moment of 
1.29 10-8 at 10° rising to 1.65 x 10-8 at 70° and 
seems to imply that the difference is due to the 
effect of the solvent upon the actual polarization 
of the alcohol molecule. Since, however, it is diffi- 
cult to picture the change in the moment of the 
alcohol molecule as occurring with sufficient ease 
to be affected by the difference in solvent, it 
seems more reasonable to adhere to the view 
originally advanced by Smyth and Stoops that 
the benzene is more effective than the heptane 
and cyclohexane in reducing the forces between 
the alcohol molecules. On this view the extra- 
polated polarization of the alcohol in heptane 
solution is much lower than the true value for the 
molecules isolated from one another at the low 
temperatures and approaches the true values with 
rising temperature, being approximately correct 
at 70°. The failure of the solvent to eliminate the 
effect of the forces between the polar molecules 
at the lower temperatures is so considerable as to 
be obvious in the case of the heptane solutions of 
butyl alcohol. A much smaller failure in the case 
of the solvent in the alkyl halide solutions could 
cause the large values observed for the atomic 
polarizations. 


10]. B. Miles, Phys. Rev. 34, 964 (1929). 
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DETERMINATIONS ON SOLIDS 


The highest atomic polarization seeming to 
possess any real significance thus far reported is 
the value 33.4 for dimethyl sulfate in the solid 
state.!! Later measurements of the dielectric con- 
stant of solid dimethyl] sulfate” give evidence of 
dipole orientation between the melting point, 
—31.4°, and a transition at —68.0—71.0°. If 
this is an orientation of dipoles within the mole- 
cule, as seems rather probable, the polarization 
arising from it must be classed as a true atomic 
polarization, which drops to a small value below 
the transition point. If, however, it arises from a 
limited orientation of the molecule as a whole, the 
atomic polarization must be estimated from the 
small dielectric constant below the transition 
point, 

The high atomic polarization, 24.8, calculated 
by Ebert for cane sugar,!* may arise from orienta- 
tion of the hydroxyl group dipoles within the 
molecule but may be due, at least in part, to the 
presence of water in the material. The same may 
be said of Ebert’s value 28.4 for citric acid con- 
taining a molecule of water of hydration, which 
probably contributes heavily to the value of Pa. 


In view of the rotation of the water molecule in 


J. Errera, Polarisation Diélectrique, p. 115, Les 
Presses Universitaires de France, Paris, 1928. 

12C, P. Smyth and C. S. Hitchcock, J. Am. Chem. Soc. 
54, 4631 (1932). 


13, Ebert, Zeits. f. physik. Chemie 114, 430 (1925). 
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ice and the great elevation of the dielectric con- 
stant by the presence of minute quantities of 
ionic impurities,” large values of P, for solids 
liable to contain traces of water must be viewed 
with scepticism. As small quantities of impurities 
may cause the presence of liquid in the solid for 
some distance below the freezing point, it is evi- 
dent that, unless the substance is very pure, the 
atomic polarization must be determined at tem- 
peratures far below the freezing point, although 
the use of high frequency in the dielectric con- 
stant measurements tends to reduce the effect of 
impurities.’ In the case of small polar molecules, 
like those of water and the hydrogen halides, 
ordinarily high frequencies and low temperatures 
would be useless in determining P4 because of the 
rotation of the molecules in the space lattice. 
However, P, may be determined for these small 
molecules from the temperature variations of the 
polarizations of the gases. 


CONCLUSIONS 


It may be concluded that improved methods 
of determining atomic polarization have demon- 
strated the incorrectness of the large values pre- 
viously reported except, perhaps, in the case of 
complex molecules where dipole orientation may 
occur within the molecule. In the calculation of a 
large dipole moment, the error caused by neglect 
of atomic polarization is normally almost 
negligible. 
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The Rate of Dissociation of Nitrogen Tetroxide' 


CHARLES EpwWIN TEETER, JR., Research Laboratory of Inorganic Chemistry, Massachusetts Institute of Technology,* 
Department of Chemistry, University of Californiat 


(Received January 17, 1933) 


It has long been supposed that the rate of dissociation 
of nitrogen tetroxide could be calculated from observations 
of the dispersion of the sound in the gas. Einstein's 
equations gave the velocity constant in terms of the 
dispersion and the properties of the substances, regarded 
as ideal, acoustically transparent gases. Dispersion has 
been reported and rate constants calculated by several 
observers. On the other hand an equal number of negative 
results have been reported. Recently the theory has been 
examined critically by D. G. C. Luck and modified to 
take account of the absorption of sound and the deviations 
from the gas laws. Moreover dispersion has been observed 
for nondissociating gases, particularly carbon dioxide. 
In this gas, time lag in the equilibrium between vibrational 
and translational-rotational degrees of freedom causes an 
increase of velocity exactly parallel with that caused by 
dissociation. Indeed any type of equilibrium in the system 
may be expected to produce a dispersion of sound in some 
frequency region related to its own frequency, and it is 
not possible to distinguish between the types of equilibrium 
by sound measurements. The use of very small measuring 


tubes also produces a dispersion of somewhat smaller 
magnitude. Under these circumstances it is difficult to see 
how any interpretation can be made or conclusions drawn. 
In a series of measurements made during the last five 
years, the author has found no measurable dispersion of 
sound in N,O, up to 53.8 kc. The velocity of sound at 
this frequency agreed closely with that determined by 
Griineisen and Goens at audible frequencies, although 
high absorption made observation difficult. At higher 
frequencies no evidence of transmission could be obtained 
although sound up to 860 ke was employed. This increasing 
absorption shows that we are entering a dispersive region 
but it is not possible to say whether dissociation, time 
lag, or a combination of both is responsible for the dis- 
persion. This applies equally to the velocity increases 
found by Richards and Reid in the same frequency region. 
Hence it seems that the sonic method for measuring rates 
of dissociation is not likely to fulfill the hopes of its 
proponents, and further experimental work will not be 
fruitful. 





I. THEORY OF THE CALCULATION OF THE 
RATE CONSTANT 


The dispersion equations of Einstein 


The speed at which the reaction N2O;—2NO, 
attains equilibrium has been a subject of in- 
vestigation for the last forty-seven years. Ordin- 
ary methods are useless for so rapid a reaction, 
and attempts have been made to determine the 
rate by measurements of the velocity of sound in 
the gas. 

The theory upon which this work was based 


* Contribution No. 26. 

{ This investigation was suggested by A. R. Olson, and 
the early portion of the work was carried out at the 
University of California with his counsel and cooperation. 
The author is indebted to the Parker Fund of Harvard 
University for a Fellowship in 1927-8, and to the National 
Research Council for a National Research Fellowship in 
1928-9, 

‘A preliminary notice of this communication has 
already appeared. (J. Am. Chem. Soc. 54, 4111 (1932).) 


was proposed by Einstein? in 1920. Equations 
were derived for the propagation of sound of 
different frequencies in a system of ideal, 
acoustically transparent gases one of which was 
the dissociation product of the other. At either 
very high or very low frequencies the velocity of 
sound is given by an exact equation not involving 
either the frequency or the rate constant of dis- 
sociation. At frequencies of the order of magni- 
tude of the dissociation frequency dispersion sets 
in, and the velocity is given by an approximate 
relation involving both the sound frequency and 
the rate constant. These equations are: 


Vo=[(p/p)(1+A/B)]}! (1) 
at high frequencies, V,,=[(p/p)(1+R/C,)]' (2) 
and at intermediate frequencies, v=w/27, 


p k2AB+C,Re\ 7} 
na ee 
p k2 B+ Cw? 


at low frequencies, 





(3) 


2 A. Einstein, Sitzungsber. d. Berl. Akad. 1920, 380. 
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By solving (3) for k; and substituting (1) and (2) 

for their values we can obtain 
ki=(22vC,/B)((V.2— V.)/(V2— Ve?) J}. (4) 


Applied to the N2O, dissociation, p is the pres- 
sure, p the density, &; the rate constant of disso- 


ciation, 
2D ao 1—a 2—a, 
A-(—-@.) +R ’ 
Z i+a a 


D? 1-a 
Ba +C, 
RT? i+a a 





_ 2-a 





D is the heat of dissociation at constant volume, 
a the degree of dissociation, C, is the weighted 
average of the heat capacities of NxO, and NO», 
C,=[C,"(1-a) +2aC,’]/(1+a), and Ris the gas 
constant. 

Several of the early experimenters*: *: *: ° 
thought they had reached the dispersive region 
at audible frequencies. From some measure- 
ments completed in 1929 the author’ believed he 
had found dispersion at 51.7 kc. Richards and 
Reid® report that they have found measurable 
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dispersion from 92 to 450 ke and calculate rate 
constants. 

On the other hand Griineisen and Goens’ 
found no dispersion up to 15.6 kc. Kistiakowsky 
and Richards” extended this to 80 ke. In the 
present paper new experiments are described 
which show no measurable dispersion up to 53.8 
kc, but do show the presence of strong absorption, 
increasing greatly up to 860 kc. This indicates 
the advent of dispersion, without affording 
quantitative evidence. 


The dispersion equations of Luck 

It seems a particularly bad approximation to 
apply the ideal gas laws to a system such as 
N.O.—NO2, where both gases deviate from 
ideality, but one (N2O,) deviates much more than 
the other. It becomes important therefore to 
examine D. G. C. Luck’s! modification of 
Einstein’s theory to apply to real gases and at the 
same time to take account of sound absorption. 
Since no equation of state was known for the gas 
mixture, Luck used a functional form, p=f(p, 
T, a). Following the general method of Einstein, 
he obtained 





V, pitex?? 7 1+6 
—-| ——) for velocity, (5) 
Vo 1+xX? (1-6)? 
d (1+e¢x?)?7i 1+6€X? 
and §=—= 1+ - for absorption. (6) 
2nl (e? —1)?X? (?@—1)X 


In these equations X (termed the relative frequency) =27v/bki, «= V../Vo, Vo=0p/dp(1+ga/b), 


V..=0p/dp(1+g), 

















2-—a Pees Cee D (= *) 
a=—+(1- ae mone | 
‘ “" ap/ao aT pM  ap/aT Co\ app 
_— DainK ap/aT Mp 

beantingd————, 9 


3E. and L. Natanson, Wied. Ann. 24, 454 (1886); 
27, 606 (1886). 

4 Argo, J. Phys. Chem. 18, 438 (1914). 

5 F, Keutel, Inaug. Diss., Berlin (1910). 

6 H. Selle, Zeits. f. physik. Chemie 104, 1 (1923). 

7A. R. Olson and C. E. Teeter, Jr., Nature 124, 444 
(1929). 

8W. T. Richards and J. A. Reid, J. Chem. Phys. 1, 
114 (1933). 





ap/dp pCr 


C,=(1—a)C,""+2aC,’ = heat capacity/mole of 
mixture, = molecular weight, / is the distance 
in which the amplitude of the sound is reduced by 


*E. Griineisen and E. Goens, Ann. d. Physik 72, 193 
(1923). 

0G. B. Kistiakowsky and W. T. Richards, J. A™ 
Chem. Soc. 52, 4661 (1930). 
u —D. G. C. Luck, Phys. Rev. 40, 440 (1932). 
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the factor e~!, and K is the dissociation constant. 
The equations reduce to the form given by 
Einstein if the equation of state becomes 
p/M-RT(i+a) and @ becomes negligible. By 
assuming values of e, Luck is able to show the 
form of the velocity and absorption curves but 
application to the calculation of rate constants 
cannot as yet be made, since the differential co- 
efficients, [0/dpla, 7, [0P/AT |a, », and [dp/da],,r 
are at present unknown and not measurable, and 
there is no equation of state from which they may 
be computed. The curve for V,/Vo differs in 
position but not in shape from that plotted from 
Einstein’s equation. That for absorption rises 
rapidly in the dispersive region and attains a 
maximum at X=1/e, thereafter dropping off 
again. 

If Luck’s equations are solved for k, the 
relation 


Qav V.2(1+62)?— V,2°(1 —62) 
b LV2(1—@)—V2(1+e)2 





is obtained. 


This is seen to reduce to the form derived from 
Einstein’s equations, 





1 


2rvC, V..7— V;? 


3 
if the gases are ideal, 


B V/?— Ve 


and to the form k,” = (27v/b)-(1+6)/(1—6) atthe 
absorption maximum. If a value of ¢€ is assumed, 
§and V,/V» at the absorption maximum can be 
calculated, while 6 and C,/B follow from the dis- 
sociation data. It is thus seen that a rough notion 
of the separate effects of the deviation from the 
gas laws and the absorption of sound can be ob- 
tained, if for K the variable K,=F(p, T, «) is 
used. The chief uncertainties are in the values of 
¢,a, and 0K/dT, in the case of ¢ since Vp and V,, 
calculated from Luck’s equations may not be the 
same as the values from Einstein’s equations. 
At the absorption maximum X=1/e and 0 
=(e-1)/(e+1). Also V,/Vo=(1+6)/2 (approx- 
imately, if eis near 1). Using values of Vp and V,, 
calculated from Einstein’s equations, we find 
that € varies from about 1.04 to 1.06 at 25°. 
With e= 1.05, V,/Vo=1.025 and 6=0.024. This 
means that the amplitude will be diminished 1/e 
in 6.6 wave-lengths. By substituting in the two 
equations for ky, k,“=2nv/bX1.050 and k,® 
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= 2nvC,/BX1.012, a difference of about 3.8 per- 
cent. If e= 1.06, 0= 0.029 and k,“= 2rv/b X 1.060, 
ki®=2nvC,/BX1.015, a slightly greater differ- 
ence. The other factors 6 and B/C, may be calcu- 
lated to show the effect of the deviations from 
ideality. At 25° and 565 mm pressure } = 86.1 and 
b! = B/C, =77.5, a difference of about 10 percent. 
At this same temperature and 260 mm b=75.7 
and b’= 65.3, and at 106 mm b= 58.1 and b’= 49.9, 
differences of roughly 14 and 13 percent, re- 
spectively. At lower temperatures the calculation 
becomes too uncertain as the values for a and K 
were experimentally determined” only at 25°, 
35° and 45°. 

Unfortunately no more information can be 
gleaned from Luck’s equations in the present 
state of our knowledge. It might be possible to 
guess at the values of the constants for an equa- 
tion of state, but the accuracy would not be suffi- 
cient to justify a calculation from Luck’s 
equations. We may state, however, that at the 
absorption maximum the values of k; from 
Einstein’s equations may be about 4 percent too 
low on account of absorption and perhaps 15 
percent too high because of deviations from the 
gas laws. 


The dispersion equations of Kneser 


At the end of Luck’s paper it was pointed out 
that dispersion due to other causes than dissocia- 
tion might render invalid any conclusions drawn 
from acoustic data. Such dispersion was ob- 
served for air and carbon dioxide by Pierce’ and 
Reid,'* and it was attributed by Herzfeld and 
Rice® to a time lag in energy transfer from trans- 
lational to internal molecular degrees of freedom. 
More lately Kneser’® studied this dispersion both 
theoretically and experimentally, and found an 
effect closely analogous to and of nearly the same 
magnitude as that predicted by the Einstein 
theory for the dissociation dispersion. In CO: the 
velocity rose about 3.5 percent in the frequency 


122F, H. Verhoek and F. Daniels, J. Am. Chem. Soc. 53, 
1250 (1931). 

13 G, W. Pierce, Proc. Am. Acad. Sci. 60, 27 (1925). 

14 C, D. Reid, Phys. Rev. 35, 814 (1930). 

%K. F. Herzfeld and F. O. Rice, Phys. Rev. 31, 691 
(1928). 

1H. O. Kneser, Ann. d. Physik (5) 11, 761 (1931); 
11, 777 (1931). 
























region 100-600 kc and again became constant at 
higher frequencies. It is instructive to compare 
the equations Kneser derived for the velocity at 
different frequencies (which gave curves checking 
his data closely) with the corresponding equa- 
tions of Einstein. For very low frequencies he 
found Vo=[(p/p)(1+R/C) ]', and for very high 
frequencies, V..=[(p/p)(1+R/C,) ]*. For inter- 
mediate frequencies in the dispersive region he 
obtained 


p C+ Ci w* \ 7} 
v.-|-(14R———)]. 
p C?+ C2 Bw* 


In these equations C is the total heat capacity per 
mole and C, is the part of the heat capacity due 
to translational and rotational degrees of free- 
dom. The symbol 8 represents the time charac- 
teristic of excitation of the vibrational degrees of 
freedom. If Kneser’s equation for V, is rewritten 


pb RC+RC,f'o\ 7} 
V,=|-—( 1+ ‘ (7’) 
C+ orc) | 


and Einstein’s equation is written 


p AB+RCh,*w"\ 7} ) 
v.=|- (1+ C )| 8) 
p B+ C*k 1?” 








the close similarity between the two dispersion 
formulas may be seen. By a fortunate choice of 
constants it would be possible to explain a given 
dispersion either as dissociation or as time lag. 

It should be pointed out that Kneser, like 
Einstein, used the ideal gas laws and neglected 
sound absorption in his derivation. Hence his 
equations should be modified in the same way to 
take account of deviations from these ideals. 
This would not be expected to alter the similarity 
of the two corrected formulas. 
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It can be readily seen that any type of equi- 
librium in the system, whether in the inter- 
molecular forces, in a chemical reaction, or in the 
states of the molecules themselves would be ex- 
pected to give a dispersion in some frequency 
region related to its own frequency. If the dis- 
persive regions overlap or coincide it is very diff- 
cult to see how any clear interpretation can be 
made or conclusions drawn. 


The dispersion equations of Richards and Reid 


Nevertheless an attempt has been made by 
Richards and Reid to modify the Einstein theory 
to include the ‘“‘heat capacity” dispersion of 
Pierce, Reid and Kneser as well as the dissocia- 
tion dispersion. The resulting equation however 
is used only to try to justify their neglect of the 
heat capacity effect. In their derivation, an ab- 
sorption term is included, but this also is neg- 
lected in all calculations. The ideal gas laws are 
still assumed throughout. With this procedure 
the equation 


VAL —B°V2/a*)/(1+6°V2/a*)*=(p/0)y (8) 


is obtained. Here 8, the ‘‘damping coefficient” 
used by Richards and Reid, is nearly identical 
with the “‘/’’ used in Luck’s equations, and 
y(=C,/C,), which is also the ratio of isothermal 
to adiabatic compressibilities (Richards and Reid 
seem to have this ratio inverted), is given by a 
complex involving the different heat capacities, 
the frequency, and also the time lag @ (Kneser’s 
B). When absorption and time lag are negligible 
(8) should reduce to the form given by Einstein. 
It is evident that there is an algebraic error in the 
reduction as carried out by Richards and Reid. 
Their Eq. (9) should read 





_ | CoRT*al_2+a(1—a) ] | V.2— a 
7 (2—a)[CoRT?(2+a(1—a@))+D*a(1—a) J | V2 — Ve? ; 


which readily reduces to our Eq. (4) above, instead of 





ee 
“1 Q-a)LORT?+D'a(t—a) J) Lv e— vel’ 


which is used in all their calculations. Fortunately, this error is independent of frequency and 


: ‘ -V 
affects only the first or smaller term of the denominator. The rate constants calculated in Tables \ I 
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and VII, however, should be decreased by about 
5 percent on this account. 


The tube dispersion 


A somewhat smaller change of velocity with 
frequency is caused, as has long been realized, by 
the use of very small tubes as the measuring 
chambers. The effect of small tubes is quite com- 
plicated. In the first place the close proximity to 
the walls causes side reflections which break up 
and distort the sound pattern and cause the ap- 
pearance of spurious nodes and antinodes. If an 
attempt to avoid these disturbances is made by 
examining only the region close to the source of 
sound, the anomalies noted by Pierce and Reid 
are likely to vitiate the results. In the second 
place the velocity, according to the Helmholtz- 
Kirchhoff formula, will decrease as the sound 
frequency or the radius of the tube diminishes. 
In a region of dissociation dispersion this change 
of velocity would be expected to increase greatly 
as Richards and Reid pointed out. At the most 
they suggested that the correction might double. 
This will not be very significant if the dispersion 
observed is large and the tube dispersion only a 
small proportion of the total. However if the 
two are of nearly the same order of magnitude it 
becomes extremely important to know the tube 
dispersion very accurately in order that the cor- 
rected value for the dissociation dispersion may 
not be meaningless. 


II. EXPERIMENTS 


Preliminary experiments 

The apparatus was particularly designed to 
avoid the complications caused by the use of 
small tubes. Sound waves were produced in a 
chamber, made as large as possible to cut down 
wall effects, and standing waves were set up be- 
tween the crystal which was the source of sound 
and a plane reflector. A sketch of the apparatus 
is given in Fig. 1 and a diagram of the electrical 
circuit in Fig. 2. The sound chamber was a bell 
of Pyrex tubing five inches in diameter, the open 
end of which was closed by a brass plate heavily 
coated with paraffin. Thecrystal,c, was contained 
in the male part of a 13 inch ground joint; the 
female part of which was sealed in a hole in the 
center of the brass plate. The lower electrode was 
a gold plate resting on a block of paraffin, the 
upper a square of platinum foil pressing lightly 
on the crystal. The seals s, s were made vacuum 
tight by applying beeswax and rosin and shellac 
on the outside. The reflector, 7, was a lantern 
slide cover supported by a glass tube and moved 
by a calibrated screw S. A vacuum-tight movable 
joint, 7, was made by stretching a piece of Gooch 
rubber tubing between rubber stoppers. To pro- 
tect the rubber and to prevent it from collapsing 
under vacuum the tube was partly filled with a 
mixture of paraffin and vaseline. The chamber, it- 
self thermostated, was in communication with a 





















































Fic. 1, Apparatus for measuring velocity of high-frequency sound in nitrogen tetroxide. Early form. 
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Fic. 2, Oscillating and resonance circuits for forcing vibration of quartz crystal. 


reservoir of liquid nitrogen tetroxide at a definite 
temperature. 

The quartz crystals were cut as suggested by 
Pierce. No success was had in inducing oscillation 
by means of his circuit, hence the power oscillator 
of Wood, Loomis and Hubbard" was employed. 
To the left of Fig. 2 is an ordinary Hartley oscil- 
lator and to the right a resonant circuit contain- 
ing the crystal (cr). The inductances L, capacities 
c, choke C and resistances R have values deter- 
mined by the tube and the crystal employed. In 
operation the Hartley oscillator was set to the 
frequency of the crystal with the aid of the stand- 
ard wave meter, WM, the circuit L’c2’ tuned to 
resonance with Lc, as indicated by the thermo- 
ammeter A;, and the two circuits varied simul- 
taneously until the crystal broke into oscillation. 


17 (a) Wood and Loomis, Phil. Mag. (VII) 4, 417 (1927). 
(b) Hubbard and Loomis, Nature 120, 189 (1927). 
(c) Hubbard and Loomis, Phil. Mag. (VII) 5, 1177 
(1928). 
(d) Hubbard and Loomis, J. Opt. Soc. Amer. 17, 295 
(1928). 





Moving the reflector had no effect on the plate 
meter of the Hartley oscillator (ZA), but caused 
a regular rise and fall in the current indicated by 
the thermomilliammeter A: coupled to the reson- 
ance circuit through a small coil Z,. Other 
methods of detection were tried but none was 
completely satisfactory. The frequency was 
measured to about 0.2 percent by the standard 
wave meter. To check the apparatus the velocity 
of sound in air was measured. Erratic results 
were obtained very near the crystal, as was noted 
above, but at distances greater than three or four 
half wave-lengths the velocity was 332 m/sec. 
at 0° in close agreement with the accepted value. 
This figure was the mean of four series, compris- 
ing seventeen measurements in all. 

The nitrogen tetroxide was prepared by heat- 
ing a mixture of sand and lead nitrate previously 
dried at 120°. A strong current of oxygen—dry 
and free from CO.—was passed through the 
combustion tube to minimize decomposition. The 
product was first collected in a trap cooled with 
solid CO;-alcohol, then distilled through a column 
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of P2O;, and finally collected in a bulb half filled 
with P.O; and cooled in liquid air or COs-alcohol. 
The product obtained from distillation was pure 
white as a solid, and reddish brown as a liquid at 
the boiling point. Permanent gases were removed 
after the bulb had been sealed onto the main ap- 
paratus, first by pumping on the solid N2O, with 
amercury diffusion pump, and then by boiling off 
two liters of the gas into an evacuated flask with 
the temperature of the liquid about 0°. The resi- 
due was analyzed for air by sealing off a 50-100 
cc bulb and noting the amount of undissolved gas 
after breaking the point of the bulb under KOH 
solution. 

The reaction on the reflector fell off rapidly 
with distance from the crystal and only a few 
half wave-lengths could be measured. In the 
single successful series of measurements, the 
value 192.6 was obtained for the velocity of sound 
of frequency 51.7 kc in N2O, at 25° and 565 mm. 
From this a rough value of the rate constant, 
ki=5600 was calculated. Analyses showed 0.1 
percent of air. 

One of the chief criticisms of the early appa- 
ratus was the use of large surfaces of paraffin and 
vaseline. It was later shown that N2O, exerts a 
slight but definite oxidizing action on these sub- 
stances, NO and N,O; being formed. This would 
tend to increase the velocity of sound. Another 
desirable modification was an increase of the 
sensitivity to permit measurements of a greater 
number of half wave-lengths and also greater 
distances from the crystal. It was also necessary 
to decrease the size of the reflector since at low 
frequencies large reflectors introduce periodic 
errors. 


Final experiments 


With all this in mind a new apparatus was 
built (see Fig. 3). The reflector was cut to 1} 
inches square. A new crystal holder was con- 
structed of nichrome IV. Brass parts were triple 
gold plated and burnished. A groove was cut in 
the end plate, and the bottom of the Pyrex tube 
ground and polished to fit in this. The joint was 
made vacuum tight by filling the outside of the 
groove with a beeswax-rosin-ferric oxide mixture 
(at s, s). For the movable joint a metal stuffing 
box was used, sealed to the upper end of the 
chamber (s, s). A rod of nichrome IV supporting 


the plane glass reflector at one end passed 
through a well fitting nichrome bearing into the 
stuffing box. This was of steel, silver brazed to 
the bearing (at a, a), and contained a disk of talc 
turned and polished to a close fit, surmounted by 
a second disk of ceresin wax, and a third disk of 
Bakelite. Pressure was applied to the topmost 
disk by the threaded top which was also tubu- 
lated to act as a guide for the rod. Even under 
vacuum it was possible to advance or withdraw 
the rod without any leak. The rod was moved by 
a standardized screw (not shown). 

The electrical circuit is shown in Fig. 4. The 
ordinary type of crystal controlled oscillator was 
employed. Values of Z, C, and C, varied with 
the crystal used. To observe the change in the 
plate current as the reflector was moved, the bulk 
of the current was balanced out by the dry cell 
and the resistances R; and Ry. Only a small frac- 
tion of the current flowed through the Ayrton 
shunt AS to the galvanometer G. The latter had 
fifty divisions, each corresponding to 1.1 micro- 
ampere. Although the amplitude of the sound 
must have been much less with this oscillator the 
reaction on the reflector in air was sufficient to 
produce full scale deflection of the galvanometer 
between nodes and antinodes. 

No change was made in the method of making 
nitrogen tetroxide, as that was deemed satisfac- 
tory. The last experiment however was made 
with gas prepared by distilling C.P. red fuming 
nitric acid in a current of oxygen through a 
column of P.O; and glass pearls and then through 
the apparatus previously used. It was found that 
practically all of the nitrogen tetroxide could be 
distilled off with the volatilization of only a small 
amount of nitric acid, and most of this was left in 
the first receiver. Further distillations under 
vacuum removed the remainder. 

After distilling the nitrogen tetroxide into the 
main apparatus it was purified even more care- 
fully than before. It was repeatedly distilled, 
frozen, left over P2O;, and pumped free of more 
volatile gases. After each experiment it was re- 
purified, and before use analyzed for permanent 
gases and examined for NO; by condensation. 
A trace of the latter gas was sufficient to show a 
slight blue color against the pure white nitrogen 
tetroxide, especially as the trioxide condensed 
first. 
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Fic. 3. Apparatus for measuring velocity of high-frequency sound in nitrogen tetroxide. Final form. 
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Fic. 4. Circuit of piezoelectric oscillator used in final experiments. 


Three crystals were employed, approximate 
frequencies (to 0.2 percent) 53.8, 100.6, and 860 
kc. These values were obtained by heterodyning 
the crystal oscillator with a second oscillator and 
measuring the frequency of the latter with a 
standard wave meter. More accurate measure- 
ments could easily have been made but were not 
deemed worth while. It should be noted that the 
frequency of a quartz crystal varies slightly with 
the temperature, and also changes with the 
mounting unless this is invariable. 

The apparatus was checked by measurements 
in air at the two lower frequencies. At 53.8 kc 
and atmospheric pressure the velocity of sound, 
corrected to 0°, was found to be 333.0 m/sec. as 
the mean of 296 determinations of the half wave- 
length, 4/2. At 100.6 kc and atmospheric pres- 
sure the mean of 411 determinations gave 332.0 
m/sec, for the velocity. Pierce obtained 332.4 
m/sec. at 50 kc and 331.8 m/sec. at 100 kc. Reid 
found values for the velocity of sound which 


decreased with increasing distance from the 
source. At a mean distance of 15 cm the velocity 
was 332.7 m/sec. at 50 kc. At 30 cm it was 332.1 
m/sec. and at 45 cm it was 331.8 m/sec. The 
mean distance in our experiment was about 10 
cm so the high result is not surprising. 

The results for nitrogen tetroxide are given in 
Tables I, II and III. The roman numerals repre- 
sent different fillings of the chamber (new samples 
of nitrogen tetroxide); the Arabic numerals 
represent series of measurements on the same 
sample. In all the experiments except VI, Table I, 
from 3 to 5 stopcocks were in contact with the 
apparatus. All were lubricated with an “‘acid- 
proof” grease. In addition to these a mercury 
manometer was used intermittently to check the 
constancy of the pressure. This was disconnected 
whenever a reading was not being made. The 
actual pressure in the chamber, however, was 
determined from the vapor pressure curve of 
nitrogen tetroxide, which seems to be fairly well 
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known." The amount of volatile decomposition 
product which reached the chamber during a run 
from these reactive parts of the apparatus must 
have been very small, as the amount of N.O; 
produced on standing overnight was barely de- 
tectable. In the experiments reported the percen- 
tage of permanent gases was always below 0.3 
percent, and no blue color of N:O; could be ob- 
served on condensation. The stopcocks and 
manometer were used simply as a convenience 
and to save time. If positive results had been 
obtained no possible precautions would have been 
spared. As conclusive evidence that the presence 
of the stopcocks and manometer did not affect 
the velocity of sound, Experiment VI, (Table I) 
was made with a redesigned apparatus eliminat- 
ing these reactive parts. The chamber and reser- 
voir were separated by a thin glass partition 
which could be broken when desired. Provision 
was made for evacuating, filling and sealing off 
the reservoir, and for evacuating and sealing off 
the chamber separately. The wave-length found 
checked very well with the results of the other 
experiments, and the average velocity of sound 
obtained from all the measurements with the 
53.8 kc crystal agreed with the low-frequency 
measurements of Griineisen and Goens.?° 

In Experiment IV with the 53.8 kc crystal a 
rough determination of the absorption was made. 
It will be recalled that in air the reaction on the 
reflector caused a change of fifty scale divisions 
on the galvanometer between node and antinode, 
even at some distance from the crystal. In nitro- 
gen tetroxide at a distance of less than } inch 
from the crystal or only one \/2 a change of less 
than fifteen divisions was observed. This fell off 
rapidly as the reflector was withdrawn from the 
crystal as shown in Table IV. In the space of 
some 6 inches or 45) the reaction fell practically 
to zero. 

At the higher frequencies no evidence of trans- 
mission could be obtained even very close to the 
crystal and although a range of pressures was ex- 

18 (a) Scheffer and Treub, Zeits. f. physik. Chemie 81, 
308 (1912). 

(b) Baume and Robert, Comptes Rendus 168, 1199, 
1201 (1919). 

(c) Macisevskii, J. Russ. Phys. Chem. Soc. (Chemical 
Part) 48, 1917 (1916). 

(d) Mittasch, Kuss, and Schlueter, Zeits. f. anorg. 
Chemie 159, 1 (1926). 
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plored. A reaction on the reflector of 0.5 scale 
divisions, or 1/100 of the effect in air, could easily 
have been detected. A confirmation that the ab- 
sorption in N2O, at 53.8 kc was much greater 
than that in air, the absorption at 100.6 greater 
than at 53.8 kc, and that at 860 kc greater than 
at 100.6 kc was the fact that at decreasing pres- 
sure limits for increasing frequencies enough 


TABLE I. Wave-length and velocity of sound in nitrogen 
tetroxide at.53.8 kc, 550 mm pressure, and 25°. 











No. of check No. \/2 Av. \/2 

Expt. series measured (scale divs.) 
I 6 93 1.348 
II 11 393 1.361 
III 7 303 1.359 
IV 5 317 1.364 
V 3 114 1.354 
VI 8 353 1.357 
Total 1220 Av. 1.360 


Velocity of sound 186.9 m/sec. 








TABLE II. Attempts to measure wave-length and velocity of 
sound in nitrogen tetroxide at 100.6 kc and 
various temperatures and pressures. 











Temp. Press. Temp. 
Expt. liquid (mm) gas Observations 
Series I 
11 —10.09 150 = 21.4 ) 
2 — 3.65 216 21.5 | No transmission; crystal 
3 + 0.07 265 21.2 oscillated strongly. 
4 + 5.01 343 ZA2 
Series II 


1 + 4.71 338 19.8 | No transmission; crystal 
2 10.52 453 20.7 oscillated strongly. 

3 15.03 565 20.3 | No transmission; crystal 
+ 19.93 714 19.9 oscillated feebly. 


II 








At pressures greater than about 450 mm the crystal 
could not be made to oscillate strongly. 


TABLE III. Attempts to measure wave-length and velocity of 
sound in nitrogen tetroxide at 860 kc and various 
temperatures and pressures. 











Press. Temp. 
Expt. (mm) gas Observations 
Series I 
I 97 25.0 \ No transmission; crystal 
II 131 25.5 oscillated strongly. 
Series II 
III 123 25.0 } No transmission; crystal 
IV 78 25.0 oscillated strongly. 








At pressures greater than about 140 mm the crystal 
could not be made to oscillate. 
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TABLE IV. Readings of screw, half wave-lengths, and change 


in galvanometer deflection in nitrogen tetroxide 
at 53.8 kc, 550 mm and 25.0°. 











Screw Av. galv. 
reading Difference No. A/2 d/2 defi. div. 
Series I 
28.95 
13.20 10 1.32 14.8 
42.15 
16.51 12 1.37 13 
58.66 
12.59 9 1.40 10 
71,25 
10.53 8 1.32 7 
81.78 
9.47 7 1.35 5 
91.25 
11.50 8 1.44 3 
102.75 
10.56 8 1.32 2 
113.31 
12.04 9 1.34 1.5 
125.35 
19.74 15 1.32 0.5 
145.09 
Series II 
31.26 
15.94 12 1.33 14 
47.20 
13.44 10 1.34 11 
60.64 
12.24 9 1.36 8 
72.88 
12.57 9 1.39 6.5 
85.45 
11.74 9 1.31 4 
97.19 
12.37 9 1.37 2.6 
109.56 
13.77 10 1.38 1.6 
123.33 
9.75 7 1.39 0.5 
133.08 
11.41 8 1.42 0.6 
144.49 








energy was absorbed to stop or greatly diminish 
the oscillation of the crystal. These limits were 
575 mm for the 53.8 kc crystal, 450 mm for the 
100.6 kc crystal and about 140 mm for the 860 
ke crystal. Of course it might have been possible 
to increase the sensitivity of detection or the 
amplitude of the sound waves, and hence force 
some of the sound through the gas. But with the 
present apparatus it seems to be clearly demon- 
strated that the absorption is too great to be 
neglected. Since greatly increasing absorption 
accompanies dispersion, as shown by the theory 
of Luck and the experiments of Kneser, it may 
be assumed that we are entering a dispersive 
region, but no clue is afforded as to the cause of 
the dispersion. 


III. Discussion 


Kistiakowsky and Richards calculated theo- 
retical values of Vo and V,, to check against their 
experimental points. Unfortunately in making 
use of Einstein’s equations for this calculation 
they failed to check his algebra and hence missed 
the slip in sign in the equation for Vo. Hence a 
new calculation was desirable. This is made in 
Table V for 25° and various pressures, and the 


TABLE V. Observed velocities of sound in nitrogen tetroxide 
at 25° and various pressures. 











(mm) uz obs.(m/sec.) b(mm) 4 obs.(m/sec.) 
Kistiakowsky and Richards Natanson 
9.9 kc 617 187.2* 
776 183.0 492 189.9* 
761 183.8 231 197.3* 
364 192.4 102 208.1* 
270 197.0t 43.6 220.2* 
41.6 kc Nernst and Keutel 
760 * 183.6T 06 188.7 
491 188.8 308 196.6* 
275 195.5f 221 199.8* 
146 203.0f 106 210.4* 
80.3 kc Argo 
760 183.5 685 189.0* 
486 189.4f Griineisen and Goens 
333 193.5 760 184. 
304 195.3 565 187.0 
Teeter 440 189.6 
53.8 ke 265 194.5* 
550 186.9 
Richards and Reid 
9 ke 92 kc 
668 185.2 185.3 
260 195.1f 195.4 
162 202.0 202.5 








* Extrapolated over 1°-5°, using temperature coefficient 
of 0.95 m/sec. /degree. 
t Average values. 


curves are plotted in Fig. 5. On this diagram are 
given all recorded points obtained at or near this 
temperature at all frequencies. In making the 
calculation, the new data of Verhoek and 
Daniels” were used. These authors give D 
= 13,960 cal./mole at constant volume" (not 
14,600 as used by Richards and Reid) and they 
use the values C,’’(N2O,)=12.75 and C,’(NOz) 
=6.8. These are not very accurate but do not 
greatly affect the results. Verhoek and Daniels’ 
data do not permit the ready calculation of 
either K or a at a given pressure. They give in- 
stead equations involving Cnooy°=g/mv, the 


19F, H. Verhoek and F. Daniels, J. Am. Chem. Soc. 
53, 1186 (1931). 
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Fic. 5. Experimental and theoretical velocities of sound in nitrogen tetroxide at 20° and various 
pressures. Upper curve, values of %,., lower curve values of u%, both calculated from Einstein’s 
theory, by using Verhoek and Daniels’ dissociation data. N, data of Natanson; K, Nernst and 
Keutel; A, Argo; G, Griineisen and Goens; R10, R40, R80, data of Kistiakowsky and Richards at 
10 kc, 40 kc, and 80 kc, respectively; 7, Teeter, R9, R92 data of Richards and Reid at 9 kc, 


and 92 kc, respectively. 
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“concentration”’ of the gas, assumed all NeO,, 
in moles per liter. This may be expressed in terms 
of P, the total pressure in atmospheres, and a 
by the relation a=(P— pnoo,°)/pneo4°, where 
pxe2og° is the pressure which would be exerted by 
undissociated N2O,, which Verhoek and Daniels 
set equal to gRT/Mv. Here again we have the 
use of the ideal gas laws for NsO,, so that the 
resulting values for a will still be somewhat in 
error. However, assuming ideal conditions to ap- 
ply, Cn2ou°= prooy?/RT= P/RT(1+a). This gives 
two equations for K, K=4a?P/(1~—a?’), and 
K=a—bP/RT(1+a), where a and }b are em- 
pirical constants found from the experimental 
curves. Equating the two values of K and solving 
for a, 


_bP/RT+ [(bP/RT)*+4(a+4P) (a—bP/RT)}} 
7 2(a+4P) 





a 


The values of a thus calculated are always 
slightly lower than the corresponding values from 
previous calculations. 

It should be noted that the points obtained by 
the authors, Griineisen and Goens, and Richards 
and Reid, define a fairly smooth curve which 
crosses the experimental curve for u», deviations 
being negative at atmospheric pressure and be- 
coming positive below 260 mm. Some of Kis- 
tiakowsky and Richards’ points come fairly close 
to this experimental curve, but others as well 
as most of the early values lie above it. Exactly 
the same type of divergence between theory and 
experiment is seen in Richards and Reid’s paper, 
Table V, where at constant pressure negative 
deviations of experiment from theory are found 
at 1°, and positive deviations at 30°. This type of 
behavior might be expected of an equilibrium 
mixture of two gases, one of which (N2O,) devi- 
ates much more from the ideal gas laws than 
the other. 

This brings us to the paper of Richards and 
Reid. Are these authors really justified in neglect- 
ing absorption, deviation from the gas laws, and 
time lag? Are their values for the velocity con- 
stant really significant? Their arguments are 
complicated and not at all clear. Furthermore the 
confusion is increased by a number of errors 
throughout the paper. The algebraic mistake in 
their final equation has already been noted. In 
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addition Tables III and IV giving the experi- 
mental values of the velocity of sound seem to be 
somewhat confused. All of the velocities are said 
to be ‘fully corrected,” yet in the last three lines 
of Table III and all of Table IV it is not possible 
to obtain the sixth column (V,— Vo) by sub- 
tracting the third column (V)) from the fifth 
column (V,). If the sixth column be assumed 
correct (and if it is incorrect their whole argument 
falls) then a number of the values of V») have not 
been corrected. It also seems that values of 
Visi— Vox) have been labeled (V.— Vo) and so 
used in the calculation whereas it is more logical 
to use values of (Vss1— V9) for the dispersion. 
The heat of dissociation at constant pressure was 
also used instead of that at constant volume. 
Let us examine their arguments in order. To 
show that @ is negligible they first compare rate 
constants calculated from Eqs. (8) and (9). It is 
not clear just how Eq. (8) is used in this calcula- 
tion, since this equation includes the effect of 
absorption as well as time lag. The magnitude of 
the former is entirely unknown and no mention is 
made of any assumption concerning it. From 
their calculations they conclude that if the con- 
stant k; calculated from experimental dispersions 
by (9) does not vary with frequency at frequen- 
cies greater than or equal to its own order of 
magnitude, the time lag can be completely neg- 
lected. They state, rather arbitrarily, that any 
dispersion as large as 1 m/sec. may be attributed 
at least in part to the dissociation. In order to 
establish further that 6 may be neglected, they 
state the necessity of obtaining concordant 
values of k; at low frequency (and dispersion) and 
at higher frequency (and dispersion). The dis- 
persion at the higher frequency must be ‘“‘sig- 
nificantly greater than the maximum which can 
be produced by heat capacity change.’’ What is 
this maximum? The maximum for CO, is about 
9 m/sec. while the greatest dispersion observed 
by Richards and Reid for N2O, was 5 m/sec. 
Furthermore a dispersion as great as 1 m/sec. is 
observed at 92 kc only at low temperatures (1°). 
The conclusion would seem to be that where ab- 
sorption can be neglected, @ cannot be and vice 
versa. Richards and Reid rext attempt to show 
that even small dispersions can be interpreted as 
wholly due to dissociation. Their argument is 
here particularly complicated and unconvincing, 
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but finally reduces to the agreement of values of 
k, at high and low frequency. 

The apparatus used does not differ very 
greatly from that used by Kistiakowsky and 
Richards. An important exception is the use of 9 
and 92 kc magnetostriction oscillators and dual 
reflectors in the same tube. This eliminates the 
difficulty of comparing different samples which 
accompanies the older method. At higher fre- 
quencies (94 and 451 kc) however separate tubes 
were used, and these were made still smaller. This 
brings us to the highly uncertain “tube correc- 
tion.’”’ This was determined by measuring the 
difference for various gases between the 9 and 92 
kc magnetostriction oscillators in 1 inch tubes and 
the 451 and 94 kc oscillators in 1 cm tubes. There 
is some question as to how the corrections of 0.4 
and 0.8 m/sec. for the different tubes were calcu- 
lated, and as to the justification of the assump- 
tion that no correction need be applied to the 
92 kc oscillator in the 1 inch tube, and that — 0.2 
m/sec. is applied to the 451 kc oscillator in the 
1 cm tube. The former is somewhat explained by 
the measured value of the velocity of sound in 
argon at 30° and 760 mm, 324.0 m/sec. “‘accord- 
ing well’’ with the value ‘324.2 calculated from 
data in International Critical Tables.” The differ- 
ence, however, is of the order of magnitude of the 
tube corrections (and of some of the dispersions 
observed). It should also be noted that in NO. 
the value of Vo, in the 1 cm tube is about 0.6 
m/sec. less than that of Vo in the 1 inch tube 
even after the corrections have been applied. 
Besides this there is the question, raised pre- 
viously, of the increased tube correction in the 
dispersion region. This correction seems especially 
important for the calculation of the rate con- 
stants at 30° at 92 and 451 kc, for it is upon their 
equality that the whole argument depends. 
Different corrections at the two frequencies 
would throw the constants apart. The same ap- 
plies to the correction for absorption which 
would increase as the absorption maximum is ap- 
proached. A very rough calculation shows that 
the absorption maximum comes somewhere near 
700 ke. Furthermore if the measured dispersion 
at 92 kc (0.2 m/sec.) is only 0.1 m/sec. too high 
(and all the velocities are accurate only to 
+0.1 m/sec. at best) the rate constant will be 
increased by nearly 50 percent. 








The final argument for neglecting all the devia- 
tions from the simple theory is a comparison of 
observed velocities with those “calculated from 
(8) when absorption and time lag are negligible’’ 
and w=0. This means the simple Einstein equa- 
tions must have been used, although this is not 
pointed out. How good is the agreement? At 
higher pressures than 260 mm it is “‘unsatisfac- 
tory” and no data are given. At 260 mm and 30°, 
the only temperature for which data at 451 kc 
are available, the observed velocity of sound is 
about 2 m/sec. greater than the theoretical. At 
1°C where dissociation is much less the observed 
velocity is about 1.1 m/sec. less than the theoret- 
ical. This latter value should be particularly 
noted as later Richards and Reid state that at 
1°C, Vo calculated from (9) agrees with the ob- 
served velocity at low frequencies within 0.1 
m/sec. We have found (9) to be in error, but since 
this equation does not give directly a value of Vi, 
they may have used some intermediate equation, 
making a direct check on their statement diff- 
cult. However the values calculated from (8) in 
Table V seem to be free from error (cf. the value 
at 25°, 194.4, with that of mu in Fig. 5 of the 
present paper). In making the calculations 
Richards and Reid used interpolated experi- 
mental data of Verhoek and Daniels for the 
highest temperatures, but the figure at 1°C in- 
volves a long extrapolation of the dissociation 
constant and the use of the ideal gas laws. 
Richards and Reid consider the agreement be- 
tween theory and experiment— +2 to —1.1 
m/sec.—“‘sufficiently good’’ to base conclusions 
on dispersions of 0.2 to 1.8 m/sec. In view of the 
large discrepancies it seems rather questionable 
to use the measured rather than the calculated 
values of Vo in Einstein’s equations, or as a 
matter of fact to use these equations at all. 

If Kneser’s equations are solved for B' we 
obtain 


B= (2nvC./C)[(V.2— V,?)/(V2— Vo") }?. 


Since C, is wholly unknown it is not possible to 
carry out this calculation accurately. Making a 
rough guess at C, and using Richards and Reid's 
dispersion data, B comes out about 5X10~° at 
92 and at 451 kc (Kneser found 1107 for 
CO). It seems therefore that an intelligible re 












RATE OF DISSOCIATION OF 


sult might be obtained by treating the dispersion 
as wholly due to time lag. 

Another point which may be raised is the dis- 
agreement between the original results of Kis- 
tiakowsky and Richards over practically the same 
frequency range and the results of Richards and 
Reid. In the former no definite frequency trend 
in the region 10-80 kc was observed. At about 
260 mm the velocity at 10 kc was nearly 2 m/sec. 
greater than that at 80 kc and that at 40 kc lay in 
between. Similar divergences were noted at other 
pressures. There is little difference apparent in 
the technique of the two series of measurements, 
but in the new series the tube correction of 0.4 
m/sec. is applied. Application of this correction 
to the old results would not establish a definite 
frequency trend. It is difficult to see why the new 
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measurements are so much better than the old, 
and Richards and Reid offer no explanation. 

To sum up it seems that the interpretation of 
the measurements of the velocity of sound is not 
nearly as easy as at one time supposed. It seems 
hazardous to neglect any of the three disturbing 
factors—absorption, deviation from the ideal 
gas laws, and time lag in equilibrium between the 
degrees of freedom. Not only does it seem that the 
time lag should not be neglected but it also seems 
possible that it can account for all the results ob- 
tained so far by Kistiakowsky and Richards, 
Richards and Reid and the present author: all 
the velocity increases being anomalous dispersion 
due to the time lag and the high absorption that 
accompany the dispersion. 
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Allotropy of Liquid Nitrobenzene 


EpMuND L. LIND AND T. FRASER YOUNG, University of Chicago 
(Received January 30, 1932) 


New measurements have been made of the surface 
tension and density of nitrobenzene between 6°C and 
20°C. No evidence of allotropy was found. Neither in- 
tensive drying with phosphorus pentoxide nor the addition 


of moisture produced any anomalous density changes, 
although some samples, cooled from higher temperatures, 
were kept for twenty-four hours at 6°C. 





ECENT work suggests that certain organic 

liquids may exist in two or more allotropic 
forms. Mazur,! and Wolfke and Mazur? report 
that the density, rate of cooling and dielectric 
constant of carbon disulfide, ethyl ether and 
nitrobenzene exhibit anomalous changes at 
various characteristic temperatures. 

Menzies and Lacoss,* studying temperature 
plots of the specific volume, refractive index, 
pressure for melting, dielectric constant and 
specific heat of benzene, show that the smooth 
curves drawn by the various experimenters to 
represent their data might better be replaced by 
broken lines. These discontinuities are at about 
40°C except in the pressure for melting. The 
vapor pressure and viscosity of benzene were 
found to exhibit similar irregularities. 

Probably because of the magnitude and the 
convenient temperature (about 9.6°C) of the re- 
ported change for nitrobenzene, a number of 
workers have investigated this substance since 
Mazur’s first notes appeared, and much contra- 
dictory evidence has resulted. Wolfke and 
Ziemecki* found that their plotted values of the 
index of refraction of nitrobenzene between the 
temperatures of 6.6°C and 14.6°C fell on a 
straight line. The discontinuity in their curve for 
molecular refraction arises from their use of 


1J. Mazur, Nature 126, 993 (1930); 127, 741 (1931); 
127, 893 (1931); Comptes Rendus des Seances de la Soc. 
Pol. de Phys. V, 349 (1931); Acta Phys. Pol. 1, 53 (1932). 

2M. Wolfke and J. Mazur, Nature 128, 584 (1931); 
Acta Phys. Pol. 1, 71 (1932); Zeits. f. Physik 74, 110 
(1932). 

3 Menzies and Lacoss, Proc. Nat. Acad. Sci. 18, 145 
(1932). 

4M. Wolfke and S. Ziemecki, Acta Phys. Pol. 1, 271 
(1932). 


Mazur’s density data in their calculations. 
Friend® likewise failed to find evidence of a dis- 
continuity in the index of refraction. Repeating 
Mazur’s density work, he obtained a smooth 
curve, whereas there is a difference of more than 
3 percent between Mazur’s density measurement 
at 5.5°C and the value extrapolated from the high 
temperature portion of his curve. Friend’s data 
of course led to a continuous function for molec- 
ular refraction. The density data of Tyrer® and 
of Massy, Warren and Wolfenden’ are in agree- 
ment with those given by Friend. 

Piekara,® in an investigation of the dielectric 
constant of carefully dried nitrobenzene, found 
no evidence of any transition above the freezing 
point, and further found that values for the su- 
per-cooled liquid fell on a curve continuous with 
that representing the data for temperatures 
above the freezing point. 

Van Itterbeek,® using nitrobenzene furnished 
by the Bureau Internationale d’Etalons physico- 
chimique of Brussels, investigated the viscosity. 
His results are very interesting, since he reports 
no break in the temperature curve of intensively 
dried nitrobenzene, whereas he does find one at 
about 9.6°C for material not so dry. Massy, 
Warren and Wolfenden’ found no discontinuity 
either in viscosity or in rate of cooling. 

In investigations of the Raman spectrum for 
evidence of allotropy, neither Wolfke and 


5 J. N. Friend, Nature 129, 471 (1932). 

6 D. Tyrer, J. Chem. Soc. 105, 2534 (1914). 

7N. B. Massy, F. L. Warren and J. H. Wolfenden, 
. Chem. Soc. 1932, p. 91. 

8 A. Piekara, Nature 130, 93 (1932). 

® A. Van Itterbeek, Nature 130, 399 (1932). 
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Ziemecki‘ nor Thorne and Bayley’? found any 
indications of a discontinuous variation with 
temperature. Stewart! reports that his x-ray dif- 
fraction measurements indicate ‘‘noticeable 
changes in the cybotactic groups” for nitroben- 
zene at 9.5°C, and rapid changes for ethyl ether 
at —105.4°C. The experiments which he reports, 
however, cannot confirm the existence of dis- 
continuities characteristic of changes of state. 

Allotropy of liquid nitrobenzene, whether due 
to changes in molecular aggregation or to other 
causes, would probably affect surface tension. 

Measurements of this property in the neigh- 
borhood of 9.5°C seemed desirable to supplement 
the existing evidence for and against a transition 
at this temperature. 


EXPERIMENTAL 


A sample of Eastman’s best grade of nitro- 
benzene was fractionated. A large portion was 
obtained which boiled within a range of 0.05°C 
at 211.,°C, and froze within a range of 0.06°C 
at 5.5°C. 

For the measurement of surface tension, a drop 
weight apparatus of the type described by Hark- 
ins and Brown,” Gans and Harkins,® and others, 
was adopted. The tip used was of Pyrex glass. 
Its diameter, as determined with a travelling 
microscope, was 0.4980 cm. The supply bottles 
were also made of Pyrex glass. 

Fresh calcium chloride was kept in the drop 
weight box at all times, and any air entering due 
to temperature changes when the box was im- 
mersed in the thermostat, passed through a cal- 
cium chloride tube. Since much of the work was 
done at temperatures lower than that of the 
room, the box was not opened after removal from 
the thermostat until it had come to temperature 
equilibrium with the room. This was necessary in 
order to prevent the condensation of moisture 
upon the contents of the box. Thermostat tem- 
peratures, maintained constant to +0.01°C, 
were determined to +0.05°C. 


10 A, M. Thorne and P. L. Bayley, Phys. Rev. 41, 376 
(1932). 

u G, W. Stewart, Phys. Rev. 39, 176 (1932). 

12 W, D. Harkins and F. E. Brown, J. Am. Chem. Soc. 
41, 499 (1919). 

13D, M. Gans and W. D. Harkins, J. Am. Chem. Soc. 
52, 2289 (1930). 
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For the determination of the factor required 
for calculating surface tension from drop weights, 
the density is required. In order that our density 
data might be compared with those of Mazur and 
of Friend, they were determined carefully with a 
Pyrex glass pycnometer. Surface tension and 
density values are given in Table I, and are 














TABLE I, 

Temperature Density Drop weight Surface tension 
6.00 1.2168 0.04369 45.33 
7.00 1.2158 0.04353 45.17 
8.00 (1.2148) 0.04347 45.11 
9.05 (1.2138) 0.04327 44.91 
9.4, (1.2134) 0.04326 44.90 
9.7 1.2131 + 

10.05 (1.2128) 0.04317 44.81 
15.00 1.2078 0.04252 44.15 
20.00 1.2025 0.04190 43.52 








represented graphically in Fig. 1, together with 
the densities reported by Mazur and by Friend. 
Density values enclosed in parentheses were in- 
terpolated. 

Density studies were continued as follows. 
Dilatometers were filled, one with nitrobenzene 
as used in our previous investigations, a second 
with slightly moist nitrobenzene, and a third with 
nitrobenzene intensively dried over phosphorus 
pentoxide; these were sealed and immersed in 
the thermostat. The last sample darkened 
somewhat on standing. Each sample contracted 
uniformly between 15°C and 6°C. Because of 
the possibility of a delayed transition, the 
dilatometers were kept at 6°C for twenty-four 
hours, but no further contraction in volume 
occurred. 


DISCUSSION OF RESULTS 


It will be observed that our density curve 
differs but little from Friend’s, and, like his, 
shows no such break as that of Mazur. Neither 
does the surface tension curve exhibit any irregu- 
larities (all points lie within less than 0.1 percent 
of the curve).!5 It must be recognized, of course, 


4 International Critical. Tables, Vol. IV, p. 435, Mc- 
Graw-Hill, New York (1928). 

15 Our value, 43.52 dyne/cm at 20°C, is in good agree 
ment with the 43.38 of Harkins, Clark and Roberts, 





ALLOTROPY OF LIQUID NITROBENZENE 


that it is possible that our material lacked a 
catalyst necessary to produce transition from a 
metastable to a stable state, or that it contained 
an impurity sufficient to mask the effect, as sug- 


J. Am. Chem. Soc. 42, 702 (1920), and the 43.35 reported 
by Hennaut-Roland and Lek, Bull. soc. chim. Belg. 40, 
177 (1931). The estimate in the International Critical 
Tables, Vol. IV, p. 454, for 20°C was based largely on 
less precise data, and now appears to be about 0.4 dyne/cm 
too high. 
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gested by Dobinski.'® On the other hand, Mazur’s 
liquid may have contained an impurity which 
produced a spurious effect? (possibly due to 
limited solubility). Further work is necessary to 
establish the existence or non-existence of dis- 
continuities in any of the properties of pure 
liquid nitrobenzene. At the present time, it 
seems to us that there is a preponderance of 


evidence against the contentions of Mazur. 


16S, Dobinski, Nature 130, 662 (1932). 
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Dispersion and Polarizability and the van der Waals Potential in the Alkali Halides 


JoserH E. MAyEr, Chemical Laboratory, Johns Hopkins University 
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It is shown that the ultraviolet absorption of NaCl, 
KCI and KI, which one may estimate from recent experi- 
mental work, is in agreement with the dispersion of these 
salts. With these assumed absorption curves, the dipole- 
dipole potential constant for the van der Waals attraction 
between negative ions can be calculated with considerable 
accuracy. The constants can be estimated for the other 
alkali halides. The quadrupole-dipole constant is also 


approximated for all alkali halides. The resulting van der 
Waals potential is much greater than that previously 
calculated and accounts for the stability of CsCl type 
lattice for CsCl, CsBr and CsI. The polarizability of a 
given ion depends on the crystal. It is here assumed that 
the polarizability varies inversely as the ‘‘main frequency” 
in the crystal. This is shown to be in approximate agree- 
ment with experiment. 





INTRODUCTION 


ONDON! has shown that a dipole-dipole 
resonance potential varying as the inverse 
sixth power of the distance exists between all 
atoms and ions. An approximation formula was 
given by London for calculating the potential. 
Born and Mayer? have used this approximation 
to compute the van der Waals potential in alkali 
halide crystals. It was there remarked that the 
approximation used appeared to give consider- 
ably too small a value. Indeed, the existence of 
caesium chloride lattices leads to the conclusion 
that the actual van der Waals potential must be 
nearly twice that calculated. Recent unpublished 
investigations indicate that this potential is 
capable of quantitatively explaining the anom- 
alous solubility and lattice constants of the silver 
halides, and that it is very large in these crystals. 
It becomes of interest to investigate the pos- 
sible methods of calculating the complete van 
der Waals potential in ionic crystals to ascertain 
if more reliable values can be obtained than those 
of Born and Mayer. Any such investigation must 
examine critically the values of the polarizabili- 
ties and absorption frequencies of ions in crystals. 
Furthermore, the dipole-quadrupole potential, in- 
versely proportional to the eighth power of the 
distance, must be considered. Equations which 
may be used to evaluate this latter term have 

been obtained by Margenau.® 

1 London, Zeits. f. physik. Chemie B11, 222 (1930). 


2 M. Born and J. E. Mayer, Zeits. f. Physik 75, 1 (1932). 
3H, Margenau, Phys. Rev. 38, 747 (1931). 
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The approximation formula for the dipole- 
dipole potential most frequently used is, for two 
ions 7 and j at a distance apart 7;; 


AE aa=—Cij/Pij®} Cig D3 (€v€jesae;/ecte;). (1) 
In this equation e;=hyv;, ¢;=hv; are energies cor- 
responding to the ‘‘main frequencies”’ »;, v; of 
the ions 7, 7. a; and a; are respectively the polar- 
izabilities of ion 7 and 7. The equation is closely 
valid if absorption frequencies of considerable in- 
tensities are found only in the immediate neigh- 
borhood of the main frequencies. e; and e; are 
usually assumed to be, and are assumed by Born 
and Mayer to be, the ionization potentials of the 
ions. In the case of singly charged positive ions 
this is the second ionization potential of the ele- 
ment. In the case of negative ions it is the electron 
affinity. 

The most obvious and serious objection to this 
use of the formula is that generally in crystals 
there is no absorption at all in the neighborhood 
of the assumed main frequencies of the negative 
ion. This paper shall be particularly concerned 
with the choice of values for ¢ and a in Eq. (1). 
A calculation of the potential constant ¢ by 
means of a more accurate method is used to help 
in this choice and in the estimation of the error 
involved in the approximation. Furthermore, 
calculations are made of the dipole-quadrupole 
constant. In addition it is pointed out that esti- 
mates from experiment and theory of the relative 
absorption coefficients are in complete agreement 
with the dispersion data for some crystals. Some 
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proposals are made regarding the calculation of 
the polarizabilities of ions in crystals. 


THE UsE oF ABSORPTION FREQUENCIES 


Experimental data 


Pfund* has measured some reflection maxima 
in the region near 1600A for NaCl. Hilsch and 
Pohl> have carried out more extensive measure- 
ments of optical absorption down to 1600A for 
many alkali halides. The results are summarized 
in the latter of the two quoted articles. Hilsch® 
has used a most interesting method of electron 
collisions, by which he obtains curves which run 
parallel to the optical absorption over the range 
where this is known. Smith’ has reported the 
absence of heavy absorption between the region 
investigated by Hilsch and by Hilsch and Pohl 
and a band at approximately 400A. Although 
this investigation was carried out only for NaCl 
the result may probably be safely generalized. 
It appears obvious, then, to associate the low- 
frequency band at approximately 1600A with the 
negative ion, and the one at around 400A with 
the positive ion. The excellent agreement between 
Hilsch’s electron collision curves and the optical 
absorption for the range above 1600A may 
further incline one to assume that the optical ab- 
sorption on the short wave-length side of this 
limit is approximately given by Hilsch’s curves. 

For several salts one may then assume an ap- 
proximate experimental knowledge of the absorp- 
tion between about 2500A and 1000A. To the 
high-frequency side of this region one may as- 
sume no or little absorption until approximately 
400A, and there some sort of absorption due to 
the positive ion. 

The assumptions are capable of fairly vigor- 
ous test in a few cases. The contribution of the 
positive ion to the dispersion curve can be esti- 
mated from the polarizability and frequency. In 
general the contribution is not large. The 
“frequency”’ of the positive ion may be assumed 
to be about that corresponding to its ionization 
potential in the gaseous state. It may more 





‘A. H. Pfund, Phys. Rev. 32, 39 (1929). 

*R. Hilsch and R. W. Pohl, Zeits. f. Physik 57, 145 
(1929) and 59, 812 (1930). 

*R. Hilsch, Zeits. f. Physik 77, 427 (1932). 

7A. Smith, Phys. Rev. 39, 1013 (1932). 


reasonably be expected to be somewhat lower, 
but even a considerable percentual difference in 
polarizability and frequency of the positive ion 
would not greatly affect the dispersion curve. It 
then remains to be seen if the experimental ab- 
sorption in the 1600A region plus the small effect 
of the positive ion will account quantitatively for 
the known dispersion. Actually in only three 
cases are the complete data at hand. The disper- 
sion has been measured for NaCl and KCI with 
great accuracy by Martens,*® and good measure- 
ments by Gyulai® exist for KBr, NaBr, KJ, 
RbCl and LiF. Of these salts a considerable por- 
tion of the absorption is known optically only for 
KI. Hilsch’s electron measurements show a 
sharply defined “‘absorption’”’ which extends con- 
siderably farther into the ultraviolet. For NaCl 
and KCI satisfactory sharply defined curves are 
given by Hilsch, although only the long wave- 
length edge of the band has been optically 
measured. 

For these three salts it is actually the case that 
analytical expressions which approximately agree 
with what is known of the absorption can be 
found, which also are in satisfactory agreement 
with the dispersion data. There appears to be no 
reason to doubt that exact agreement with the 
dispersion could be obtained without assuming 
an absorption in contradiction to the experi- 
mental values. Indeed such agreement can prob- 
ably be found with so many different absorption 
functions that it appears to be of no value to 
continue further than has been done in this paper 
until more accurate experimental absorption 
curves far into the ultraviolet are at hand. 

With a knowledge of the absorption curve, and 
of the dispersion (which latter is needed to fix 
the absolute value of the transition probabilities, 
since the absorption is measured in arbitrary 
units), one is in a position to make more exact 
calculation of the van der Waals potential be- 
tween the negative ions in the crystals whose ab- 
sorption is known. The relatively certain values 
so obtained can then be compared with the ap- 
proximation formula (1) with a view to ascertain- 
ing the energies which one may best choose in 
such an equation. 


8F, F, Martens, Ann. d. Physik [4], 6, 603 (1901); 
8, 459 (1902). 
*Z. Gyulai, Zeits. f. Physik 46, 80 (1928). 
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The exact equations 

In crystals one has to deal with more or less continuous absorption. If we define as q?(v) the sum 
of all the amplitudes square up to the wave number », for one kind of ion in the crystal, say the 
negative, 

VoR=¥ 
P(v)= LD Goxquos (2) 
k, v4=0 

then the derivative with respect to v, d[q(v) ]/dv will be proportional to gogo at vorx=v. The ab- 
sorption coefficient at v due to the negative ion, will be proportional to voxgoxgxo, that is, to 














vd[_q?(v) ]/dv. 
The following equations hold: 
dl g’(v) 
Out= [ dp a= [7% Wo(>09)WVodr, (3) 
' coordinate all 
space electrons 
822mc pr d[q(vr) ] 
= wl vy ———— dv=number of electrons, (4) 
3h 0 dv 
Bi d{q*(v) } 
a=— —_ dv= polarizability (5) 
3heoJ4y v dv 
= f {dL.g?(vi) ]/dvs} {dL g?(v;) ]/dv;} — 
Ci= dv,dv;=dipole-dipole potential constant for like ions (6) 
= | Vitr; 


in which m is the mass of one electron and e¢ its 
charge, h Planck’s constant, and c the velocity of 
light. v is wave number and has the dimensions 
cm~!. Qoo” is the 00 matrix element of the sum of 
the coordinates squared for all electrons in the 
ion. In Eg. (6), were the constant c for unlike 
ions, it would be changed only by using two dif- 
ferent functions d[q*(v) ]/dv. 

If the integrations indicated are carried out 
from zero to infinity the equations are exact, 
providing the function d[q*(v)]/dv is exactly 
known. If the integrations are carried out only 
over the range of the ultraviolet where apparently 
almost all of the absorption due to the negative 
ion, except at x-ray frequencies, exists, a certain 
approximation is introduced. First obviously Qo” 
will be the 00 matrix element of 2g’ for at least 
only the outer electrons. Further a small propor- 
tion of the transitions of the outer electrons will 
be expected to be ‘‘double jumps” and will be 
far beyond in the ultraviolet. That is, Qo? so 
calculated will be somewhat smaller than the 
(Zq*)oo for all the outer electrons. The electron 
number p so calculated may be expected to be 
considerably smaller than the actual number of 








outer electrons since the small transition proba- 
bilities at high frequencies omitted become of 
considerable importance when multiplied by the 
wave number. 

However the crystals are actually transparent 
for all frequencies below the range of considera- 
tion. Those functions, a and ¢;;, as well as the dis- 
persion, which have » in the denominator will be 
given accurately in spite of the omission of some 
high-frequency transitions. 

The dispersion due to the negative ion will be 
given by 


D*(v)n 





© yild vi) \/dv; 
-c f = =s : ne (D 


If however in this equation D*(v),, is interpreted 
as being the additive contribution of the negative 
ion to the dispersion function D(v) =(n?—1) 

(n?+2), although an equation of the above form 
would be indicated, the frequencies would not be 
those experimentally observed in absorption. 
The function D*(v) =n?—1 will also be given by 
an equation of the form of (7) and here the 
frequencies will actually be the experimentally 
determined absorption frequencies, although the 
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constant C will not have a simple significance; 
that is m?—1 is not additively composed of the 
contributions of each ion.’ D*(v), is then to be 
compared to a function corresponding to n?—1 
due to the negative ion, which is obtained in the 
following manner. 

D(v)», the contribution to D(v)=(n?—1)/ 
(n?+-2) of the positive ion, is assumed to have the 
“single frequency” form D(v),=a/(v,?—v*). vp is 
taken as the frequency corresponding to the 
second ionization potential of the element. a is 
chosen to be in agreement with the polarizability 
of the positive ion, Va/vy,=42Na/3, with V the 
mol volume of the crystal and N Avogadro's 
number. The experimental D(v),=D(v)—D(v)>p 
is then found, the contribution to (m?—1)/(n?+2) 
of the negative ion, and from this D*(v), 
=3D(v),/(1—D(v)n) is calculated, the value of 
n’—1 if only the negative ions were present in the 
crystal. Since the absorption coefficient is only 
measured in arbitrary units, no attempt is made 
to define or determine further the constant C. 


Calculations 


It is now desired to select an analytical expres- 
sion for d[_q?(v) |/dv for which all the integrations 
can be conveniently performed, and such that 
vd{q’(v) |/dv is approximately equal to the ex- 
perimentally determined absorption coefficient. 
The expression selected can be tested by com- 
parison of D*(v), calculated from (7) with the 
“experimental” values obtained by the method of 
the preceding paragraph. 

The expression d[q?(v) ]/dv is assumed to be 
zero up to a certain wave number A, from A to 
increase linearly with vy to a maximum value 3d at 
v= B (except that at v= lying near A it takes an 
additional value bf) due to a sharp line at v), 
and then to decrease linearly with v from »=B 





to zero at v=C. From C to infinity it is assumed 
to be zero. 

The absorption assumed is shown in Fig. 1 
for NaCl, KCl and KI, along with what is known 
experimentally. The relative value of the in- 
tensity of the line at ») is obtained by forcing the 
function to give the correct ratio to two extreme 
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Fic. 1, Absorption curves for NaCl, KCl and KI. 
Unbroken curve, analytical expression; broken curve, op- 
tical absorption; curve with bars, electron collision method. 


















































points on the dispersion curve. It is seen that this 
demand does not lead to an unreasonable value of 
the relative absorption at vp; indeed the complete 
omission of this line would change the results but 
little. 
Using, then, this analytical expression 
d[q?(v) |/dv=0 from v=0 to v=A 
=[b/(B—A) ](v—A) 
from v=A to v=B 
=[b/(C—B) |(C—») 
from v=B to v=C 


=0 from v=C to v= @ 


plus bfo at v= vo, one obtains 


di q*(v) ] 

f —— dy=0[(C—A)/2+f] (3’) 
om dl gv b 

f , dO] =! (C-A)(A+B+O), (4’) 
0 dv 6 

“1 dl av 8 C A B 0 

f 2. db | —— tog "tog = +“ (S’) 
o » dy C-8 “8 B-4°A4 & 


Gtteeteescneeee 


See K. F, Herzfeld and K. L. Wolf, Ann. d. Physik 78, 35 (1925). 
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ei {dl q?(vs) ]/dvs} {dL g?(»;) J/dy;} So 2fo 
0 Jo vit; 2% (C—B)(B-A) 
+(A—C)(B+ 9) log (B+) +(B—A)(C+ 0) log (C+ v0) | 
4[1/6(C—B)*(B—A)*][(C—B)*(2A)? log 2A +(A —C)*(2B)? log 2B 
4. (B—A)?(2C)* log 2C-+-2(B—A)(A —C)(B+C)? log (C+B) 
+2(C—B)(B—A)(A+C)' log (A+C)+2(C—B)(A —C)(A +B)? log (A +B)]}, (6’) 
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dndr,-P [(C—B)(A+ 9) log (A +79) 































© y:{d[ q?(v;:) |/dy; b 1 C+y C-v 
f {dLg’(vi) ]/ ane | [ (c+) log“ ) (C-») pale ] 
Jo ve—yv 2\(C-B (B+y) (B—v) 
1 (B+) (B—v) 2vofo 
so v) 1 A-yvy)1 : ; 
7 [a+ ) log aaa ) log a (7’) 













The actual magnitude of the absorption, the In Table I are all assumed and calculated 
value of b, is determined by the experimental values. The constant fy is determined from the 
value of the polarizability. With the definition of _ ratio of D*(v), at the highest and lowest frequen- 
D(v), as the contribution of the negative ion to cies, and the undefined constant C in Eq. (7) 
(n?—1)/(n?+2), the polarizability of the negative from the value of D*(v), at one of these frequen- 
ion, cies. The value of D*(v), calculated for an arbi- 
cn =(3/4eN)-V-D(v), (at v=0). (8) trary frequency in the middle of the experimental 
‘range is compared with the experimental value 
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Salt 
















































































1 | AX104 cm 5.70 5.56 4.17 
2 } Assumed characteristic wave numbers of | BX 104 cm7 8.00 7.71 5.74 
3) the amplitude—squared function CX 104 cm 13.00 12.17 10.45 
4\ vo X 104 cm- 6.45 6.15 4.555 
5 fo (calculated from the ratio of (10) to (14) cm cm7 1370 274 960 

6 6X10 (calculated from X,) cm 0.895 1.06 1.69 
7 apos. X 10% (Pauling)*® cm? 0.182 0.844 0.844 
8 vp, X10 (from ionization potential) cm7 38.2 25.7 25.7 

9 aneg.X10* (from exp. dispersion) cm? 3.09 3.30 6.57 
10 D*(v), experimental 2.411 1.830 2.602 
11 atAinA 1854 1854 2540 
12 D*(v), experimental 1.720 1.245 1.914 
13 atAinA 2312 2312 3130 
14 D*(v), experimental 1.268 0.901 1.460 
15 atAinA 6563 6563 6150 

16 D*(v), calculated 1.708 1.258 1.910 
17 atAin A 2312 2312 3130 

18 Qoo?X 101° cm? 3.40 3.54 5.46 
19 » (electron number of neg. ion) 3.25 3.24 4.00 
20 ¢ X10 (dipole-dipole constant for like neg. ions) 117 131 400 
21 X10"? ergs/molecule 16.4 16.1 12.3 
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(see rows 12 and 16). The agreement in all cases 
is within 0.3 percent in , not, it is true, within 
the experimental error in 1, but well within the 
expected error of the assumptions. It appears 
useless to attempt more complicated assumptions 
about the form of the absorption to get better 
agreement until more extensive experimental 
data are at hand. 

The electron number of the chloride ion in 
NaCl and KCI have the same value, not however 
the value found by Herzfeld and Wolf" for the 
corresponding rare gas atom Ar, p=4.61, calcu- 
lated by a similar method from the dispersion 
curve for the gas. The value found here for iodine, 
4.00, is likewise not in agreement with their 
value 5.61 for Xe. Furthermore if the electron 
number of the positive ions is calculated from the 
polarizability, assuming a single main frequency 
at or near the ionization value, a poor approxima- 
tion but the only available one, values much 
higher than those of the corresponding rare gases 
are obtained. It appears then inadvisable to as- 
sume a constant electron number value for a 
given configuration. 

In row 21 and 22 are given values for e, and »,, 
the energy and corresponding wave number, 
which, with the experimental polarizability in 
the respective salts will give the values of c, 
the dipole-dipole potential constant, as calculated 
here, if inserted in Eq. (1). That is, Eq. (1) will 
give the correct constant c if the e used is the e, 
of Table I. One of the main purposes of this 
paper is the determination of «,. 


The determination of the frequencies 


For the three salts, NaCl, KCI and KI, by 
means of exact equations, using assumptions 
which are unlikely to be in great error, the dipole- 
dipole potential constant has been calculated for 
the negative ions. By a correct choice of “‘main 
frequencies” the approximation Eq. (1) will give 
the same result. What are the corresponding 
frequencies in other salts? 

To this problem it seems difficult to give a 
generally satisfactory answer. For the alkali 
halides the absorption frequencies observed ap- 
pear to show sufficient empirical regularity to 
enable reliable guesses to be made. The obtaining 


"'K. F, Herzfeld and K. L. Wolf, Ann. d. Physik 76, 
71 (1925), 
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of any sort of semi-theoretical basis which would 
enable extension to other salts is much more 
difficult. 

Hilsch and Pohl have observed that the first 
pronounced absorption maxima at the long wave- 
length end agree well with the energy ae?/r 
+J—E. a is the Madelungs number, ae?/r the 
electrostatic potential on one ion of the crystal, J 
the electron affinity of the negative ion, and E 
the (first) ionization potential of the metal. 
Born” has pointed out that this equation corre- 
sponds to the neutralization of two ions on the 
crystal surface, or on the surface of an interior 
split in the crystal. The corresponding (2a—1)e?/r 
+ J—E for the neutralization of two neighboring 
ions within the lattice'* appears not to correspond 
to any intense absorption. The use of ae?/r 
+J—V in which V represents the electron 
affinity of the crystal and may be calculated from 
magnetic susceptibilities has been suggested by 
Pauling.“* The agreement is by no means con- 
vincing. 

In Table II is given a list of the frequencies ex- 
perimentally observed. A, the lowest frequency 
at which absorption appears, and vw the first 
sharp line are taken from the optical data. The 
nomenclature is in agreement with the equation 
assumed previously for d[q?(v)]/dv. B when 
bracketted is guessed at from the optical data, 
although it is never certain that the curve may 
not later rise; when not bracketted it is taken 
from Hilsch.® C is always taken from the data of 
Hilsch. In column 6 is given the value proposed 
for e.. For comparison ae?/r, J and E are given in 
columns 7, 8 and 9, and ae®?/r+J—E and 
(2a—1)e?/r+J—E in columns 10 and 11, 
respectively. 

Use of the ratio between v. and » found for the 
three salts and recorded in row 23 of Table I has 
been made in estimating e,. The values of e¢. may 
probably be relied upon to give the dipole-dipole 
constant for the potential varying as the inverse 
sixth power of the distance in Eq. (1) to within 
about 15 percent. 

The positive ion frequencies are almost en- 


tirely unknown experimentally. The value of 


2 M. Born, Zeits. f. Physik 79, 62 (1932). 

3K. F. Herzfeld, Handbuch der Experimentellen Physik, 
Vol. VII, Part 2, “‘Gittertheorie.” 

4L. Pauling, Phys. Rev. 34, 954 (1925). 
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TABLE II. Frequencies of the alkali halides ergs /molecule X 10". 


























1 2 3 4 5 6 8 9 10 11 
Salt A vo B Cc Ec ae*/r J E (ae?/r) +J—E (2a—1)(e?/r)+J—E 
LiF (18.3) 23.2 24.3 19.81 6.58 8.55 17.84 26.31 
Cl 17.5 15.51 5.97 ” 12.93 19.56 
Br 8.7 10.62 15.3 14.42 5.62 ” 11.49 17.65 
I 7.4 8.89 (11.69) 13.3 13.10 5.12 * 9.67 15.27 
NaF 17.05 19.6 31.9 22.2 17.22 6.58 8.15 15.65 23.01 
Cl 11.2 12.68 15.7 29.9 16.4 14.18 5.97 ss 12.00 18.06 
Br 9.2 10.33 14.5 13.39 5.62 10.86 16.58 
I 7.4 8.56 (11.55) 12.8 12.35 5.12 - 9.32 14.60 
KF 19.6 14.93 6.58 6.87 14.64 21.02 
Cl 11.0 12.05 15.1 23.9 16.1 12.71 5.97 ¥“ 11.81 17.24 
Br 9.4 10.46 14.2 12.11 5.62 ay 10.86 16.04 
I 8.2 8.96 11.25 20.5 12.3 11.30 5.12 = 9.55 14.38 
RbF 18.5 14.14 6.58 6.62 14.10 20.14 
Cl 10.68 11.76 15.7 12.20 5.97 11.55 16.77 
Br 9.43 10.22 13.8 11.60 5.62 3 10.60 15.56 
I 7.81 8.73 (10.42) 12.0 10.86 ele 9.36 14.00 
CsF 19.0 13.26 6.58 6.16 13.68 19.35 
Cl 10.30 12.04 16.0 11.28 5.97 = 11.09 16.40 
Br 9.32 10.52 13.9 10.81 5.62 96 10.27 14.95 
I 7.63 8.96 (10.61) 12.1 10.19 5.12 9.15 13.56 








about 400A given by Smith’ for the NaCl crystal 
is some 65 percent of the series limit for the 
gaseous ion. One theoretical remark can be made, 
whereas removal of an electron from a negative 
ion neutralizes a lattice point and requires the 
work ae?/r, removal of an electron from a positive 
ion creates a doubly charged lattice point and 
releases the work ae?/r. Among the effects tend- 
ing to raise the frequency of the negative ion toa 
higher value in the crystal than in the gas is 
probably a term ae?/r, this same term will prob- 
ably tend to reduce the frequency of the positive 
ion in the crystal below the value in the gas. The 
relatively high frequencies of the positive ion 
will probably lead to a much smaller relative 
change in frequency with the change from gas to 
crystal, and a greater relative constancy of the 
frequency of one ion in different crystals. 


POLARIZABILITIES 


In the previous discussion it has been assumed 
that the polarizabilities of the positive ions are 
known, and are constants of the ion. The polar- 
izabilities of the negative ions have been obtained 
by subtracting from the extrapolated value of 
(m?—1)/(n?+2) at v=0 for the crystal, the con- 
tribution due to the positive ion. This method is 





justified only by the relatively smaller contribu- 
tion of the positive ion and its probably smaller 
variation with the crystal than that of the nega- 
tive ion. 

Actually considerable doubt exists as to the 
values of the polarizabilities of the positive ions. 
Many measurements of the sum of the polariza- 
bilities of positive and negative ions in solutions 
and crystals exist, and various estimates of the 
polarizabilities of the individual ions have been 
made. The values of Fajans and Joos" are doubt- 
less quite good, and they have observed and re- 
marked on the variation of these values with the 
state. 

In this paper the polarizabilities calculated by 
Pauling'® which are ostensibly those of the 
gaseous ions, are taken as fairly reliable average 
values for the positive ions in the crystals. 

If Eqs. (3) and (5) are examined it is seen that, 
if the usual approximation of a single effective 
frequency is assumed, one may write 


a = (2e?/3hcv) Qoo? = (2e?/3) (Qoo?/e). (9) 


In the three salts investigated, NaCl, KCI and 
Ki the correct ¢ to use in this approximation was 


found to be one to two percent higher than the ¢ 


16K, Fajans and G. Joos, Zeits. f. Physik 23, 1 (1924). 
161, Pauling, Proc. Roy. Soc. A114, 181 (1927). 
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to be used in Eq. (1). Herzfeld and Wolf" have 
used a similar approximation equation a =h’e?p/ 
4nr°me*, assuming the electron number p to be a 
constant for a given ion (indeed for a given con- 
figuration) and the polarizability a to vary in- 
versely as the square of the main frequency. Al- 
though the excellent constancy of the electron 
number # for Cl- in NaCl and KCl lends support 
to this view it appears to give too great a varia- 
tion of a with the frequency in most salts. It 
seems reasonable however to assume that in 
Eq. (11) the Qo?, which represent a sort of size of 
the ion, will be unaffected by inclusion in a crys- 
tal, and will be a constant for any given ion. The 
polarizability would then vary inversely as the 
frequency. The frequencies of the positive ions 
are probably appreciably decreased from that of 
the gaseous state by inclusion in a crystal. The 
frequency in the crystalline state is however likely 
to be relatively independent of the crystal. For 
this reason the polarizability of the positive ion in 
the crystal is very probably fairly constant from 
crystal to crystal and somewhat higher than in 
the gaseous state. The values given by Pauling 
are probably better approximations for the 
polarizabilities in the crystalline state than in 
the gaseous. 

Assuming that the polarizabilities of the posi- 
tive ions are constant, and those given by Paul- 
ing, we may test the assumption that Qg? is a 
constant for a given ion by calculating the extra- 
polated (m?—1)/(n?+2) function at v=0 for all 
alkali halides, and comparing with the experi- 
mental values. A value of Qo? for each negative 
ion is chosen so as to give the best average agree- 
ment. Table IV gives the results of such a calcula- 
tion. The values of (m?—1)/(n?+2) at »=0 are 
taken from Born and Heisenberg!® and the ex- 
perimental sum of the two polarizabilities calcu- 
lated by means of Eq. (8). The values used for 
Qoc? are found in Table III. It is to be seen from 
Table IV that the agreement is not very good. 
Nevertheless deviations from pure additivity of 
the polarizabilities are between three and four 
times as great as those shown in Table IV. In 
view of the doubtful positive ion polarizability 





“K. F. Herzfeld and K. L. Wolf, Ann. d. Physik 78, 
195 (1925). 

'*M. Born and W. Heisenberg, Zeits. f. Physik 23, 
407 (1924), 


















































TABLE III. 

Ion F- cr Br- Sag 
Qoo? X 1016 1.44 3.40 4.08 5.27 
TABLE IV. 

1 2 3 4 5 6 7 

Experimental Calculated 

mM—1 Ya X10% apog, X10" aneg. X10% La X10™ Diff. 3 
Salt n2+2 cm’ cm3 cm cm3 and 6 
LiF 0.234 0.906 0.030 0.90 0.93 +0.02 
Cc 0.358 2.91 2.95 2.98 + .07 
Br 0.404 4.00 4.06 4.09 + .09 
I 0.450 5.93 6.00 6.03 + .10 
NaF 0.198 1.17 0.182 0.98 1.16 —0.01 
cl 0.306 3.27 3.14 3.32 + .05 
Br 0.346 4.38 4.29 4.47 + .09 
I 0.394 6.36 6.24 6.42 + .06 
KF 0.218 1.97 0.844 1.11 1.95 —0.02 
Cl 0.280 4.14 3.20 4.04 — .10 
Br 0.310 5.30 4.38 5.22 — .08 
I 0.354 7.41 6.49 7.33 — .08 
RbF 0.237, 2.59 1.42 1.18 2.60 +0.01 
Cl 0.282 4.71 3.29 4.71 — .00 
Br 0.308 5.93 4.50 5.92 — 01 
I 0.346 8.15 6.65 8.07 — .08 
CsF 0.278 3.60 2.45 1.15 3.60 —0.00 
0.244 5.70 3.22 5.67 — .03 
Br 0.373 7.03 4.47 6.92 — .il 
0.404 9.14 6.60 9.05 — .09 








and the uncertain frequencies it is probable that 
better agreement could not be expected, even if 
the assumption were accurate. Better agreement 
could obviously have been obtained by assuming 
higher polarizabilities for the heavy positive ions, 
K+, Rb* and Cst*. 


THE VAN DER WAALS POTENTIAL 


It is now desirable to alter the approximation 
Eq. (1) used in calculating the dipole-dipole at- 
tractive potential to allow for the variable polar- 
izability. An approximation equation is also de- 
sired for the quadrupole-dipole potential. These 
potentials are then calculated for the ions in the 
alkali halides. 

Using Eqs. (1) and (9) one obtains: 


AE aa= — (Cij/1ij®) 3 Ci; = Fe4(Qioe?Qjo0?/ec+e;), (10) 


in which the Qj? is supposed to be a constant of 
the ion 7 and e¢; is dependent on the ion and on 
the crystal. 

The values of ¢ for the negative ions are the e, 
values tabulated in Table II, column 6. For thé 
positive ion a fair value is probably about 
€=0.75¢ if €& is the ionization potential of the 
gaseous ion. 
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The complete potential in the crystal is 


Cc C44+-C_- 
er iia ace iale wai (11) 


r 


in which the sums S,’ and S,’’ are sums, respec- 
tively, over the odd and even lattice points of 
the crystal.? In Table V are tabulated the values 


TABLE V. van der Waals constants for the alkali halides. 








C X 106, ergs-cmé d X10, ergs-cm® 





Salt 
D 


11. 
104. 
190. 
470. 


31. 
180. 
300. 
630. 


150. 
560. 
800. 
1420. 


290. 
960. 
1340. 
2240. 


600. 
2600. 
3600. 
5800. 


i) 
+ 


MORN BwWnod | 


6 as em 


14.5 
111. 


C44 
0.073 


185. 
378. 








of c,_, the potential constant for a single positive 
and a single negative ion, c,4 the potential con- 
stant for two positive ions, and c__ that for two 
negative ions. C for the crystal is also tabulated. 
Margenau*® has calculated the quadrupole- 
dipole and quadrupole-quadrupole potentials, 
proportional to the inverse eighth and tenth 
powers of the distance, respectively, for like 
atoms. The equations can be easily generalized 
for unlike atoms or ions, 7 and 7. The quadrupole- 
quadrupole term when evaluated is found to be 
altogether negligible and will not be here treated. 
The approximation corresponding to Eq. (10) 
for the quadrupole-dipole potential between 7 

and j at a distance apart 1;; is given by 

ME a= — 4ii[7;*; 
dij=e LQ inc’ Qjoot+Q:*Qjo0? J. (12) 
e+ €;j 

The energies ¢; and ¢; in Eq. (12) are not 
necessarily the same as in Eq. (10), the correct 
averages which should be used as probably some- 
what higher than those to be used in Eq. (10). 


E. MAYER 


Nevertheless for the purpose of the present ap- 
proximation it will be assumed that they are the 
same. Qoo' is, analogously to Qgo?, the 00 matrix 
element of the sum of g¢‘ for all electrons. The 
following equations are then nonrigorous, but 
will hold approximately: 

Qoo? = Lgoo? = Pgoo’, 

a 


ll 
electrons 


Qoot = > go0t = Pgoo', 
all 
electrons 


(13) 


(14) 


(15) 
(16) 


Qoo* =(doo")?, 
Qoot =(1/P)(Qoo*)?, 


in which # is the number of electrons. Eq. (12) 
may then be rewritten, 


dij = 3C:;(Qin0?/ Pit Qjo0?/D;). (17) 


For p the geometrical mean of the electron 
number calculated from the polarizability and 
the frequency, and theoretical number of outer 
electrons (8 or 2) is taken. The complete quadru- 
pole-dipole potential in a crystal becomes 


Uqa= —D/?r; 


D=S¢'ds~-+ Ss" (d44+d__/2). (18) 


The values of d,_, d,4 and d__ as well as D are 
tabulated in Table V. The values of the sums 
Se’, Se’, Ss’ and S3’’ for both sodium chloride and 
caesium chloride type lattices are given in 
Table VI. The arithmetical values of the sums 


TABLE VI. 








CsCl 


8.7088 
3.5445 
8.2077 
2.1476 


Lattice type NaCl 


So’ 6.5952 
Se” 1.8067 
Ss’ 6.1457 
Ss!” 0.8002 











given correspond to the use of the shortest lattice 
distance, the distance between unlike neighbor- 
ing ions, for the r of Eqs. (11) and (18). 

The values of the constants d and D are, of 
course, highly uncertain. They may be in error 
by as much as 50 percent of their value. How- 
ever, no more accurate method of estimation ap- 
pears to be available at present. The total 
quadrupole-dipole potential in the alkali halide 
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crystals, D/r’, is in all cases between 10 percent 
and 20 percent of the dipole-dipole potential 
C/r®. Since this latter term is itself in the nature 
of a small correction to the electrostatic at- 
tractive potential, the error introduced in calcu- 
lations of the total potential by uncertainties in 
D is small. 

A comparison of the values for CsI here calcu- 
lated, C=2970X10-® erg Xcm*® plus D/r? =370 


HALIDES 279 
x10-®, giving 334010-® with the value 
C=1740X10® used by Born and Mayer is in- 
teresting. The newly computed values are fully 
capable of explaining the stability of the caesium 
chloride type lattice for CsI, CsBr and CsCl. 

The author is greatly indebted to Professors 
Max Born and K. F. Herzfeld who have both 
helped greatly in developing some of the ideas 
expressed here. 
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Starting from the discovery by Rumer that the eigen- 
functions corresponding to different distributions of valence 
bonds in a molecule can be represented by plane diagrams 
which provide information regarding their mutual linear 


independence, a very simple graphical method is developed 
for calculating the coefficients of the integrals occurring 
in the matrix elements involved in Slater’s treatment of 
the electronic structure of molecules. 





N his valuable paper! on Molecular Energy 
Levels and Valence Bonds Slater developed a 
method for formulating approximate eigenfunc- 
tions for molecules and for constructing the cor- 
responding secular equation. This method has 
been successfully applied by a number of in- 
vestigators in the discussion of general valence 
theory and of the structure of individual simple 
molecules, its application to more complex mole- 
cules having been retarded in part by the labori- 
ous nature of the calculation of the coefficients of 
the various exchange integrals in the matrix ele- 
ments. It is shown in the following paragraphs 
that these coefficients can be easily calculated by 
a simple procedure involving the use of plane 
diagrams similar to those introduced by Rumer? 
in his recent discussion of the linear independence 
of electronic structures. 





SINGLET STATES WITH SPIN DEGENERACY ONLY 


For a system involving WN electrons, Slater 
constructed the function 


1 
¥= Oy! Die(—1)*"P"A(1)B(2)--- EW), 


which is completely antisymmetric in the elec- 
trons (satisfying Pauli’s principle). Here A, 
B.--E are one-electron spin-orbit functions, and 
the symbol P” represents the NV! permutations of 
the functions A---E among the electrons, the 
coefficient being —1 for odd, +1 for even permu- 
tations. For two such functions WY and W’ the 
matrix element corresponding to a dynamical 
quantity F is 


(1/N') f {Ze(—1)F"P"A(A)- »-E(N)} F{Dp(—1)?'P'A'(1)- + -E'(N) fda, 


which reduces to 


fad)---EQNFtDe(- 


with P’=P’’P, in case that F is completely sym- 
metric in the electrons. This may be written as 


Yp(—1)?(AB---E/F/PA’B’:--E’), 


in which each parenthesis (A---E/F/PA’---E’) 

represents one integral in the above sum. If F 

does not involve spin interactions, the integral 
1J. C. Slater, Phys. Rev. 38, 1109 (1931). 


*G. Rumer, Nachr. d. Ges. d. Wiss. zu Géttingen, 
M. P. Klasse, p. 337 (1932). 
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1)?PA’(1)-+-E’(N)}do, 





(A---E/F/A'---E’) vanishes unless the spin of 
A is the same as that of A’, B as B’, and so on, 
in which case it reduces to (ab---e/F/a’b’---e’), 
with a---e’ the orbital parts of A---£Z’. 

For four electrons, for example, with only spin 
degeneracy (the number of occupied orbits 
equalling the number of electrons), Slater gave 
the function }(W1— Vi1— Vii1+ Viv) as represent- 
ing the structure in which orbits a and 0 are 
bonded together, and also c and d. Here W1- - - iv 
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are functions of the type given above with the 
following distributions of spins among the orbits: 


a b Cc d 
I - + - + 
i+ - - + 
III _ + ao _ 
IV ao ~ aa _ 


The energy matrix element for this function is 


(abed/H/abcd) + (abcd /H/bacd) + (abcd /H/abdc) 
—43(abcd/H/cbhad) —}(abcd/H/dbca) 
—4}(abcd/H/achd) —}3(abcd/H/adcb) +higher 

exchange integrals. 


With the single exchange integrals negative, as 
they usually are for orbits on separate atoms, this 
leads to attraction between a and 6 and between 
c and d, and repulsion between other pairs. 

Similarly a function can be formulated repre- 
senting a—d and b—c bonds, and one representing 
a—c and b—d bonds. These three functions, 
corresponding to three separate electronic struc- 
tures of the Lewis type, are, however, not linearly 
independent; any one can be represented as a 
combination (sum or difference) of the other 
two. This linear dependence of structures con- 
tinues as the number of bonds increases; for 
single bonds between 2” atoms (2n)!/2"n! 
structures can be drawn, of which only (2)!/ 
n\(n+1)! are independent. In the paper men- 
tioned above, Rumer made the very interesting 
and important observation that if the orbits are 
formally arranged in a ring or other concave 
curve, and straight lines are drawn between 
orbits bonded together, the structures repre- 
sented by diagrams in which no lines intersect 
are independent. This observation forms the 
starting point for the following treatment. 

Let us arrange the orbits in a ring, and then 
draw arrows between pairs of orbits bonded 
together, the head of each arrow representing 





((2m) !)4 


































positive spin and the tail negative spin. The 
wave function for the corresponding structure is 
then given, except for a normalizing factor 2”, 
by the sum of the functions corresponding to this 
distribution of spins and to those obtained by 
reversing the arrows, each function having the 
coefficient (—1)*, in which R is the number of 


reversals. Thus the 7 structure would be rep- 


resented by the vector-bond diagram =, corre- 
sponding to the function 


a(s*)-(F=)-AS)+ Ca) 


It is verified at once that 
xX=It-2, 


so that an intersecting pair of arrows can be re- 
placed by the difference of two non-intersecting 
pairs, the pair connecting head with head and 
tail with tail having the positive sign. By succes- 
sive applications of this treatment any structure 
can be resolved into structures involving no in- 
tersecting bonds. 


Canonical sets of independent structures 


Of the yarious complete sets of independent 
structures which may be formulated, certain ones 
may be called ‘‘canonical’’ because of the relative 
simplicity of the calculation of matrix elements 
based on them. For a given order of the orbits in 
a ring, the canonical set of independent struc- 
tures comprises those with no intersecting bonds; 
moreover, numbering the orbits in order about 
the ring, the canonical vector-bond structures for 
singlet states with spin degeneracy only have 
arrows drawn from odd to even orbits (the func- 
tions for the corresponding distributions of posi- 
tive and negative spins thus occurring with the 
positive sign). The eigenfunction for a structure 
represented by a vector-bond diagram a—b 
c—d, etc., may be written 


1 1 
ma Le(— VR L(—1)*Pa(1)B(1)b(2)a(2)¢(3)B(3)d(4)a(4)- + |. 


in which @ represents positive spin and B nega- 
lve spin, P represents the (27)! permutations of 
the orbits and their associated spins among the 





electrons, and R the 2” operations of interchang- 
ing the spin functions @ and 6 for orbits (such as 
a and b) which are bonded together, that is, of 








reversing the arrows in the vector-bond diagram. 
The factor (—1)*” equals +1 for an even number 
of reversals, —1 for an odd number. 


The Coulomb coefficient 


To find the coefficient of the Coulomb integral 
for two structures, superimpose their vector- 
bond patterns to form the superposition pattern 
(Fig. 1). The Coulomb coefficient is 2~” times the 
sum (—1)* for the different patterns in which 
each orbit serves either as the head or as the tail 


~~ et ie. 


a ae 


= oo att 
IIOUD 


Fic. 1. The vector-bond diagrams for three structures of 
the canonical set of fourteen for »=4, and some of their 
superposition patterns. 


of two arrows which can be made by reversing 
the arrows, R being the number of reversals. The 
superposition pattern consists of ‘‘islands,” or 
closed polygons formed by an even number of ar- 
rows. For each island there are only two satisfac- 
tory orientations of the arrows, so that the 
Coulomb coefficient may be written as (—1)’ 
2-“—®, in which 7 is the number of islands and r 
the number of reversals required for an initial 
satisfactory orientation. 

For canonical structures the sign is positive, 
the Coulomb coefficient being 1/2"-‘. In drawing 
the superposition pattern for canonical structures 
the arrows may be replaced by lines. 


Coefficients of exchange integrals 

Each of the N! permutations may be repre- 
sented by a permutation diagram, made by 
drawing an arrow from each orbit in Pab---e to 
the orbit whose place in ab- - -e it occupies. Thus 
if Pab---e=bd---e, an arrow is drawn from b 
to a, from d to J, etc. A closed polygon of arrows 
is called a cycle of the permutation.’ A cycle may 


*The theory of permutation groups is discussed in 
standard treatises such as W. Burnside, The Theory of 
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be called odd if it involves an even number of 
arrows, otherwise even. The permutation is odd 
if it contains an odd number of odd cycles; other- 
wise it is even. The possible permutations are ob- 
tained by drawing polygons of arrows in all 
possible ways such that each orbit is reached by 
no arrow or by the head of one and the tail of 
another. The permutation arrows may intersect. 

To find the coefficient of a given exchange 
integral in a matrix element, (I/F/PII), draw 
the vector-bond diagram for structure II, change 
it as indicated by the permutation diagram for P, 
and form the superposition pattern of I and PII. 
The coefficient is then given, except for the factor 
(—1)”, by the above rules for the Coulomb co- 
efficient; that is, it is (—1)?(—1)"2-@-®, 

This procedure may be simplified considerably 
for single exchange integrals in matrix elements 
between structures of a canonical set. We first 
form the superposition pattern of the two struc- 
tures. If the permutation involves the interchange 
of orbits in different islands, the number of 
islands is decreased by one; if the orbits are an 
odd number of bonds apart in the same island, it 
is unchanged; if they are an even number apart, 
it is increased by one. Hence (taking into con- 
sideration the factors (— 1)” and (—1)’) it is seen 
that the coefficient of a single exchange integral 
for two canonical structures is 2~“~—°f(p), where 
n is the number of bonds, 7 is the number of 
islands in the superposition pattern of the two 
structures, p the number of bonds in the super- 
position pattern along the path between the two 
orbits which are interchanged, and f(p) a factor 
with the values —3 for p=0, +1 for p=1, 3, 
5---, and —2 for p=2, 4, 6---. It is not neces- 
sary to use arrows in forming the superposition 
pattern for canonical structures. 

Similar rules may be formulated for more com- 
plex permutations. 

An example: Six electrons with only spin de- 
generacy. Placing the six orbits in a ring in the 
order abcdef, the five structures forming a canon- 
ical set are those given in Fig. 2. It is seen that 


JS 1 Fidv=Q+ab+cd+ef —4ac—}ad—4ae— af 
— bc — thd — be — bf — 3ce— cf — 3de—3df 
a oer, 


Groups, Second Edition, Cambridge University Press, 
1911. 
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e 


een 


ra i 
aan ee a ee 


= 
Topsy 
‘ae Ny Z a 


Fic. 2. The five canonical structures for n=3, and some of 
their superposition patterns. 


in which Q is the Coulomb integral (ab- - -f/F/ab 
---f), ab the single exchange integral (ab---f/ 
F/ba---f), etc. As in this case, it is true in gen- 
eral that in diagonal matrix elements the coeffi- 
cient-of bonding single exchange integrals is +1, 
of nonbonding —4. From the superposition pat- 
tern of I and II we obtain 


J IFIldv = }(Q+ab+bc+cd+de+ef+fatad 
+be+cf—2ac —2bd — 2ce— 2df — 2ea — 2fb 
+... -), 


and similar expressions for the other matrix 
elements can be quickly written.* The determina- 
tion of the coefficients of each of the 720 integrals 
would be, of course, quite a task, even by this 


simple method, but any one of them is easily 


found. Thus from Fig. 3 it is seen that in I FIIdv 


oy a i SS 


m3, 692, 763, (UJ. 


Fic. 3. The permutation diagram for P abcdef = bcfdae, and 
the superposition pattern for I PII. 


‘A. Sherman and H. Eyring (J. Am. Chem. Soc. 54, 
2661 (1932)) have published matrix elements for this 
six-electron system, giving the Coulomb and single 
exchange integrals. Their coefficients do not show the 
regularities which our treatment leads to, since their 
five functions do not form a canonical set. 


























the coefficient of (abcdef/F/bcfdae) is —}. The 
complete matrix for the forty-two independent 
structures given by ten electrons with only spin 
degeneracy, retaining only Coulomb and single 
exchange integrals, can be set up in a few hours. 


EXTENSION TO MORE GENERAL SYSTEMS 


The foregoing treatment can be extended to 
states of arbitrary multiplicity by the artifice of 
including in the system additional electrons and 
phantom orbits. To represent N electrons and NV 
orbits in states with resultant spin S (multiplicity 
2S+1), form a ring of the real and phantom 
orbits; a canonical set of electronic structures is 
then obtained by drawing in all possible ways 
such that no bonds intersect a phantom bond be- 
tween each of 2.S phantom orbits and a real orbit 
and a real bond between each pair of remaining 
real orbits. Rules regarding reversal of bond ar- 
rows (different numbers of arrow heads on phan- 
tom orbits correspond to different Ms values) 
and the calculation of exchange integral coeffi- 
cients are easily formulated. 

A canonical set of structures for a system with 
more orbits than electrons is obtained by arrang- 
ing all the orbits (including phantom orbits for 
S>0) in a ring and then drawing non-intersecting 
bonds to a number determined by the number of 
electrons and the multiplicity. If two electrons 
occupy the same orbit, forming an unshared pair, 
a loop is drawn with its ends at the orbit. 

The bond diagrams provide an obvious simple 
method of determining the allowed spectral terms 
for equivalent electrons with Russell-Saunders 
coupling, which may be convenient for the reason 
that it separates states of different multiplicity at 
the start. 

The methods developed in this paper have been 
applied in a discussion of the structure of aro- 
matic substances, free radicals, etc., to be pub- 
lished soon. 
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The Thermal Decomposition of Ethyl Mercaptan and Ethyl Sulphide 


W. M. MALIsorF AND E. M. Marks, Allantic Refining Company, Philadelphia, Pennsylvania 


N the first issue of this journal (Volume 1, page 
77), in their paper with the above title, Trenner 
and Taylor make several references to our own 
work on a similar subject. We believe they have 
misrepresented us on a sufficiently large number 
of points to give a completely misleading impres- 
sion of its nature to any one reading their article. 
For example, according to Trenner and Taylor: 
“It has frequently been suggested that the ali- 
phatic mercaptans behave thermally in a manner 
similar to the amines, and might therefore be sus- 
pected of decomposing in a unimolecular manner. 
In fact, that is the supposed finding of Malisoff 
and Marks recently for propyl, butyl, amyl and 
heptyl mercaptans, both iso and normal.” 

Here is what we say in our summary: “‘It is 
impossible with data of this order of accuracy to 
decide finally between a first and a second-order 
reaction. The trend is much closer to a first-order 
reaction and is confirmed by an approximate in- 
dependence of the rate of decomposition of the 
concentration used and the simultaneous poor 
agreement with tests for bimolecularity.”’ And in 
the body of our article we say ‘“‘This leaves much 
room for vacillation in deciding between first and 
second-order reactions,’ and we express our 
doubts in some detail. 

According to Trenner and Taylor: aa 
hardly surprising that Malisoff and Marks 
should find no effect from hydrogen and carbon 
dioxide.” 

What we said in our summary: ‘“The decompo- 
sition of amyl mercaptan is to some extent de- 
pressed by the admixture of gases like hydrogen 
and carbon dioxide.’’ And in the article itself: 
“It was found that gases like hydrogen and car- 
bon dioxide have no marked effect on the order 
of magnitude of the amount decomposed.”’ 

Our conclusions are “hardly justifiable’’ be- 
cause ‘‘the intermediate, as has been shown, 
reacts analytically as a mercaptan.’’ Despite 
their impression to the contrary Trenner and 
Taylor have not shown the presence of an inter- 
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mediate; the compound they obtained might 
have been one, but it might just as well have 
been a secondary product. And concerning this 
compound, which is supposed to make the results 
of our analyses worthless, they themselves say 
that ‘‘as has been shown, a maximum concentra- 
tion of the intermediate is found in the neighbor- 
hood of the point of inflection’’—and in this 
neighborhood, working at pressures of the order 
of atmospheric, they obtained ‘‘small yields,” be- 
cause “its rate of formation is only slightly 
greater than its rate of decomposition.’’ To secure 
enough product for an analysis to be made they 
worked at lower temperature and at 150 at- 
mospheres. The mercaptans in our experiments 
were at pressures of the order of 10 millimeters. 

The difference between the temperature coeffi- 
cients found by us and by Trenner and Taylor 
may be ascribable to other things than our “‘un- 
justifiable averaging of numerous erratic results.” 
There is, for one thing, the difference in condi- 
tions of experiment, and, for another, the fact 
that Trenner and Taylor had first to assume that 
certain points on their curves were ‘‘comparable” 
and then to determine these points by inspection 
—a process of a rather subjective nature, as can 
be seen from their Table VI, which contains 
points of inflection read from Fig. 3. Our values 
were determined directly from the amount de- 
composed found by analysis. 

Trenner and Taylor talk of our ‘‘conclusions 
that mercaptan decomposition involves a straight 
split into olefine and hydrogen sulfide in a manner 
satisfactorily accounted for by a single uni- 
molecular rate.’”’ The ‘‘straight split’ was a sug- 
gestion on our part, but was so far from being a 
conclusion that we spend about a quarter page in 
discussing a possible chain mechanism. 

It is chiefly to this sort of reference by Trenner 
and Taylor that we object. We make a suggestion 
concerning the mechanism, a suggestion justified 
by the experimental results, with full knowledge 
of its nature as an assumption, and we state that 
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the evidence favors this rather than the other 
hypothesis so far advanced. We do not consider it 
a conclusion, in fact we advance still another 
hypothesis. But Trenner and Taylor speak of our 
“conclusions” and ‘“‘finding,’’ and having thus 
set up a straw man proceed to demolish him. 
This seems to us decided misrepresentation, 
especially as Trenner and Taylor proceed from 
their supposed refutation to cast doubt upon the 
validity of our experimental results. 

As for the other details of criticism on the page 
which Trenner and Taylor devote to our results, 
it is well to realize that this article of theirs deals 
with the decomposition of ethyl mercaptan in the 
absence of solvents, by a static method, at pres- 
sures mostly of the order of several hundred mm, 


DECOMPOSITION OF ETHYL MERCAPTAN 





285 


with times of contact of the order of a half hour, 
at temperatures from 386°—420°. Our experi- 
ments were done with different mercaptans, in 
hydrocarbon diluents where the total mercaptan 
pressure in the gas was of the order of 10 mm, 
with times of contact from a few seconds to three 
minutes, at temperatures from 300° to 475°. We 
analyzed the products, they simply followed the 
pressure changes. 

It is well to repeat that the title of our paper 
contains the words “in hydrocarbon solvents” 
and all our conclusions refer to decomposition 
under those conditions only. It is quite possible 
for some of them to diverge from results of 
Trenner and Taylor without our being entirely 
incorrect. 
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The Thermal Decomposition of Ethyl Mercaptan and Ethyl Sulphide. A Reply 


NELSON R. TRENNER AND H. Austin Taytor, Department of Chemistry, New York University 


N reply to the communication of Malisoff and 

Marks, it is our belief that a detailed perusal of 
their paper! would lead anyone to believe that 
they had found the decomposition of mercaptans 
to be (1) entirely homogeneous, (2) a simple spiit 
into an olefine and hydrogen sulphide, and (3) 
of order between one and two, probably closer to 
the former. Disregarding entirely our suggested 
mechanism our actual experimental results can- 
not substantiate these three points and therefore 
some criticism is justified. 

We believe that, accepting the above three 
points, the logically scientific conclusion would 
be that the mercaptan decomposition was a uni- 
molecular reaction proceeding in a concentration 
range where a deviation towards second order 
had already set in. We believe further that this 
was the intent of their paper, although in the 
summary they do not so commit themselves. 
Why else should they refer to ‘“The unimolecular 
nature of such decompositions... ,” ‘‘For com- 
parison with other unimolecular reactions... ,” 
and finally, ‘Should it be necessary to invoke 
chain reactions to explain the speed of such uni- 
molecular reactions .-.. ,” this latter in spite of 
the complete absence of any evidence of a chain 
mechanism? 

We demonstrated in our paper that an im- 
portant feature of mercaptan decomposition lies 


1 Malisoff and Marks, Ind. Eng. Chem. 23, 1114 (1931). 





in the existence of an induction period, markedly 
affected by surface and foreign gases when the 
latter were present in amounts comparable with 
the amount of mercaptan, though not necessarily 
so in a system already excessively diluted with 
other foreign gases namely solvent vapors. This 
induction period increased rapidly with a de- 
crease in initial pressure and should therefore 
have been quite obvious at the low partial pres- 
sure of 10 mm used by Malisoff and Marks. 

Regarding their temperature coefficient which 
they claim was not an averaged value a study of 
their Fig. 6 will show an energy of activation of 
23,000 calories for ten percent decomposition 
with obviously higher values for fifteen and 
twenty percent. The value claimed by them is 
33,500 calories. A haphazard choice of two other 
experiments. 7528 and 7463 yields a value of 
43,700 calories. Finally it may be added that a 
decision on mechanism is not necessary to the 
calculation of a temperature coefficient and since 
we had found experimentally a comparable point 
on all our reaction rate curves, our calculation 
was entirely justified. 

In conclusion it would appear from our find- 
ings that the effect of ten gases was not such as to 
alter materially the latter course of the reaction 
and the assumption therefore that mercaptan 
decomposition in hydrocarbon solvents would be 
different from that in the pure state is unjustified. 
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The Thermal Decomposition of Ethylene Iodide. Errata 


L. B. ARNOLD, JR. AND G. B. Kistiakowsky, Chemical Laboratory, Harvard University 


N the abstract on page 166, Vol. I, February 
issue of this Journal, ‘the two reactions are 
31,400 cal. and 37,500 cal., respectively.” 
should read, ‘‘the two reactions are 30,200 cal. 


and 36,600 cal., respectively.” 

On page 169, ‘‘an activation energy of 31,400 
cal.” should read, “an activation energy of 
30,200 cal.” 


Zeeman Effect in Solids. Errata 


F. H. SpeppinG, Department of Chemistry, University of California 


ECAUSE of an accident in the mails, the 
author did not personally read the proof of 
his article in the February issue before it went to 
press. Therefore the following errors, which 
ordinarily would have been corrected in the 
proof, appeared in the publication. 
Page 145, col. 2, line 12, read 0.9A° instead of 
9.A°. 
Page 148, col. 2, second curve from top, the 
field of 18,500 gauss was parallel to the ‘‘a” axis 
of the crystal. In all other curves on the page it 
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was parallel to the ‘‘c”’ axis. 

Page 149, col. 1, third line from bottom, read 
“‘’’ factors instead of three factors. 

Pages 151 and 152, read “‘a”’ axis instead of 
“d”’ axis. 

Page 153, above curves in col. 1, insert light 
ray parallel to ‘‘b’”’ axis. Above column 2 light 
ray parallel to ‘‘c’’ axis. 

Page 153. The lower spectrum is reversed and 
should be turned end for end to correspond with 


the spectrum with zero field. 
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Editorial 


In response to requests from our contributors, and after 
a canvass of the opinions of the editorial board, the Journal 
of Chemical Physics wishes to announce that a ‘‘Letters to 
the Editor” section will be carried by the Journal in the 
future. This section will accept reports of new work, 
provided these are terse and contain few figures and espe- 
cially few halftone cuts. The Editorial Board will not hold 
itself responsible for opinions expressed by the cor- 
respondents. 


Contributions to this section must reach the office of the 
Managing Editor not later than the 15th of the month 
preceding that of the issue in which the letter is to appear. 
No proof will be sent to the authors. The usual page charge 
($3.00 per page) will not be made and no reprints will be 
furnished free. 

Haro_p C. UREy 
Managing Editor 


The Isotopic Fractionation of Water* 


When water is subjected to electrolysis an isotopic frac- 
tionation occurs (Washburn and Urey, Proc. Nat. Acad. 
Sci. 18, 496 (1932)). It has been found that the heavier 
isotope of hydrogen and the heavier isotopes of oxygen are 
concentrated in the residual water. The specific gravity of 
the residual water rises continuously as the electrolysis 
proceeds. The rise in specific gravity is accompanied by a 
rise in the freezing point and in the boiling point and by a 
decrease in the refractive index. Thus for water of specific 
gravity 1.0014, the changes are; F.P., +0.050°C; B.P., 
+0.02°C; np*, —(60+2) x 107%. 

No indication of approach to an electrolysis equilibrium 
has been found and there is every reason to hope that it will 
be possible to obtain the various isotopes of hydrogen and 
oxygen in a pure state, certainly in highly concentrated 
form. 

If the oxygen from the electrolysis of normal water is 
combined with normal hydrogen, the water produced has a 
lower specific gravity than normal water. If this water be 
again partially electrolyzed and the oxygen combined again 
vith normal hydrogen, a further drop in the specific grav- 





ity of the water occurs, the total drop thus far obtained 
being 18 parts per million. It will consequently be possible 
to prepare isotopically pure water of the composition 
H'0"'*H! and hence to determine accurately the atomic 
weight of normal oxygen on the O'*= 16 scale. 

With the different isotopes of hydrogen and oxygen avail- 
able in pure form, a new field of chemistry and possibly 
also of biology will be opened up, since the different iso- 
topes of hydrogen, at least, may be expected to exhibit 
pronounced differences in chemical behavior. 

A survey is being made of water from different natural 
sources to find out whether differences in isotopic composi- 
tion occur in nature. 

EDWARD W. WASHBURN 
EpGAR R. SMITH 
MIKKEL FRANDSEN 
Washington, D. C. 
March 6, 1933 


* Published by permission of the Acting Director of the 
Bureau of Standards, Department of Commerce. 


VOLUME 1 








ao. = a= - 


— + 





